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Preface to the Second Russian Edition 


In the mid-twentieth century the theory of partial differential equations was 
considered the summit of mathematics, both because of the difficulty and 
significance of the problems it solved and because it came into existence later 
than most areas of mathematics. 

Nowadays many are inclined to look disparagingly at this remarkable area 
of mathematics as an old-fashioned art of juggling inequalities or as a testing 
ground for applications of functional analysis. Courses in this subject have 
even disappeared from the obligatory program of many universities (for ex- 
ample, in Paris). Moreover, such remarkable textbooks as the classical three- 
volume work of Goursat have been removed as superfluous from the library of 
the University of Paris-7 (and only through my own intervention was it possi- 
ble to save them, along with the lectures of Klein, Picard, Hermite, Darboux, 
Jordan, ...). 

The cause of this degeneration of an important general mathematical the- 
ory into an endless stream of papers bearing titles like “On a property of 
a solution of a boundary-value problem for an equation” is most likely the 
attempt to create a unified, all-encompassing, superabstract “theory of every- 
thing.” 

The principal source of partial differential equations is found in the 
continuous-medium models of mathematical and theoretical physics. Attempts 
to extend the remarkable achievements of mathematical physics to systems 
that match its models only formally lead to complicated theories that are 
difficult to visualize as a whole, just as attempts to extend the geometry of 
second-order surfaces and the algebra of quadratic forms to objects of higher 
degrees quickly leads to the detritus of algebraic geometry with its discourag- 
ing hierarchy of complicated degeneracies and answers that can be computed 
only theoretically. 

The situation is even worse in the theory of partial differential equations: 
here the difficulties of commutative algebraic geometry are inextricably bound 
up with noncommutative differential algebra, in addition to which the topo- 
logical and analytic problems that arise are profoundly nontrivial. 
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At the same time, general physical principles and also general concepts 
such as energy, the variational principle, Huygens’ principle, the Lagrangian, 
the Legendre transformation, the Hamiltonian, eigenvalues and eigenfunc- 
tions, wave-particle duality, dispersion relations, and fundamental solutions 
interact elegantly in numerous highly important problems of mathematical 
physics. The study of these problems motivated the development of large 
areas of mathematics such as the theory of Fourier series and integrals, 
functional analysis, algebraic geometry, symplectic and contact topology, the 
theory of asymptotics of integrals, microlocal analysis, the index theory of 
(pseudo-)differential operators, and so forth. 

Familiarity with these fundamental mathematical ideas is, in my view, 
absolutely essential for every working mathematician. The exclusion of them 
from the university mathematical curriculum, which has occurred and contin- 
ues to occur in many Western universities under the influence of the axiomati- 
cist /scholastics (who know nothing about applications and have no desire to 
know anything except the “abstract nonsense” of the algebraists) seems to me 
to be an extremely dangerous consequence of Bourbakization of both math- 
ematics and its teaching. The effort to destroy this unnecessary scholastic 
pseudoscience is a natural and proper reaction of society (including scientific 
society) to the irresponsible and self-destructive aggressiveness of the “super- 
pure” mathematicians educated in the spirit of Hardy and Bourbaki. 

The author of this very short course of lectures has attempted to make 
students of mathematics with minimal knowledge (linear algebra and the ba- 
sics of analysis, including ordinary differential equations) acquainted with a 
kaleidoscope of fundamental ideas of mathematics and physics. Instead of the 
principle of maximal generality that is usual in mathematical books the author 
has attempted to adhere to the principle of minimal generality, according to 
which every idea should first be clearly understood in the simplest situation; 
only then can the method developed be extended to more complicated cases. 

Although it is usually simpler to prove a general fact than to prove nu- 
merous special cases of it, for a student the content of a mathematical theory 
is never larger than the set of examples that are thoroughly understood. That 
is why it is examples and ideas, rather than general theorems and axioms, 
that form the basis of this book. The examination problems at the end of the 
course form an essential part of it. 

Particular attention has been devoted to the interaction of the subject 
with other areas of mathematics: the geometry of manifolds, symplectic and 
contact geometry, complex analysis, calculus of variations, and topology. The 
author has aimed at a student who is eager to learn, but hopes that through 
this book even professional mathematicians in other specialties can become 
acquainted with the basic and therefore simple ideas of mathematical physics 
and the theory of partial differential equations. 

The present course of lectures was delivered to third-year students in the 
Mathematical College of the Independent University of Moscow during the 
fall semester of the 1994/1995 academic year, Lectures 4 and 5 having beeu 
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delivered by Yu.S. Il’yashenko and Lecture 8 by A.G. Khovanskiï. All the 
lectures were written up by V.M. Imaïkin, and the assembled lectures were 
then revised by the author. The author is deeply grateful to all of them. 

The first edition of this course appeared in 1995, published by the press of 
the Mathematical College of the Independent University of Moscow. A number 
of additions and corrections have been made in the present edition. 
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Lecture Í 


The General Theory for One First-Order 
Equation 


In contrast to ordinary differential equations, there is no unified theory of 
partial differential equations. Some equations have their own theories, while 
others have no theory at all. The reason for this complexity is a more com- 
plicated geometry. In the case of an ordinary differential equation a locally 
integrable vector field (that is, one having integral curves) is defined on a 
manifold. For a partial differential equation a subspace of the tangent space 
of dimension greater than 1 is defined at each point of the manifold. As is 
known, even a field of two-dimensional planes in three-dimensional space is in 
general not integrable. 


Example. In a space with coordinates x, y, and z we consider the field of 
planes given by the equation dz = ydz. (This gives a linear equation for the 
coordinates of the tangent vector at each point, and that equation determines 
a plane.) 


Problem 1. Draw this field of planes and prove that it has no integral surface, 
that is, no surface whose tangent plane at every point coincides with the plane 
of the field. 


Thus integrable fields of planes are an exceptional phenomenon. 

An integral submanifold of a field of tangent subspaces on a manifold is a 
submanifold whose tangent plane at each point is contained in the subspace 
of the field. If an integral submanifold can be drawn, its dimension usually 
does not coincide with that of the planes of the field. 

In this lecture we shall consider a case in which there is a complete theory, 
namely the case of one first-order equation. From the physical point of view 
this case is the duality that occurs in describing a phenomenon using waves or 
particles. The field satisfies a certain first-order partial differential equation, 
the evolution of the particles is described by ordinary differential equations, 
and there is a method of reducing the partial differential equation to a system 
of ordinary differential equations; in that way one can reduce the study of 
wave propagation to the study of the evolution of particles. 
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We shall write everything in a local coordinate system: x = (£1,..., £n) 
are the coordinates (independent variables), y = u(x) is an unknown function 
of the coordinates. The letter y by itself denotes the coordinate on the axis of 
values, and we denote the partial derivatives by the letter p: p; = su Sigs 


The general first-order partial differential equation has the form 


EE che Dis eats a) ea Ors 


Examples. 


Ou 
— —0; 1.1 


Ou\? (du ? 
— — = 1.2 
(2) f (=) (12) 
(the eikonal equation in geometric optics); 


Ut + Uug = 0 (1.3) 


(Euler’s equation). 


Consider a convex closed curve in the plane with coordinates +1, £2. Out- 
side the region bounded by the curve we consider the function u whose value 
at each point is the distance from that point to the curve. The function u is 
smooth. 


Theorem 1. The function u satisfies the equation (1.2). 


PROOF. Equation (1.2) says that the square-norm of the gradient of u 
equals 1. We recall the geometric meaning of the gradient. It is a vector 
pointing in the direction of maximal rate of increase of the function, and its 
length is that maximal rate of increase. The assertion of the theorem is now 
obvious. O 


Problem 2. a) Prove that any solution of the equation (1.2) is locally the 
sum of a constant and the distance to some curve. 

b) Understand where the wave-particle duality occurs in this situation. 
(In case of difficulty see below p. 14, Fig. 2.2.) 


Consider the field u(t, x) of velocities of particles moving freely along a line 
(Fig. 1.1). The law of free motion of a particle has the form z = y(t) = zo +vt, 
where v is the velocity of the particle. The function y satisfies Newton’s 
equation dg = 0. We now give a description of the motion in terms of the 
velocity field u: by definition ae = u(t, p(t)). We differentiate with respect 
to t, obtaining the Euler equation: 
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ult ae) 


Fig. 1.1. A particle on a line . 


Conversely, Newton’s equation can be derived from Euler’s, that is, these 
descriptions of the motion using Euler’s equation for a field and Newton’s 
equation for particles are equivalent. We shall also construct a procedure 
for the general case that makes it possible to reduce equations for waves 
to equations for the evolution of particles. First, however, we consider some 
simpler examples of linear equations. 


1. Let v = v(x) be a vector field on a manifold or in a region of Euclidean 
space. Consider the equation Ly(u) = 0, where the operator L, denotes the 
derivative in the direction of the vector field (the Lie derivative). 

In coordinates this equation has the form v ge veep Un $e = 0; it is 
called a homogeneous linear first-order partial differential equation. 

For the function u to be a solution of this equation it is necessary and 
sufficient that u be constant along the phase curves of the field v. Thus the 


solutions of our equation are the first integrals of the field. 


n 
For example, consider the field v = 5° Ti g- in Fig. 1.2. Let us solve the 
i=l 
equation L,(u) = 0 for this field v. The phase curves are the rays x = ezo 
emanating from the origin. The solution must be constant along each such 
ray. If we require continuity at the origin, we find that the only solutions are 
constants. The constants form a one-dimensional vector space. (The solutions 
of a linear equation necessarily form a vector space.) 


Sa 


Fig. 1.2. An Eulerian field ae 


In contrast to this example, the solutions of a linear partial differential 
equation in general form an infinite-dimensional space. For example, for the 
Ou 


equation a = 0 the solution space coincides with the space of functions of 


n — 1 variables: 


Uu (L508 ng Bei) 


It turns out that the same is true for an equation in general position in a 
neighborhood of a regular point. 
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The Cauchy Problem. Consider a smooth hypersurface T”! in x-space. 
The Cauchy problem is the following: find a solution of the equation L,(u) = 0 
that coincides with a given function on this hypersurface (Fig. 1.3). 


u pr-l = UO 


Fig. 1.3. The Cauchy problem 


r”-! 


A point of the hypersurface is called noncharacteristic if the field v is 
transversal to the surface at that point. 


Theorem 2. The Cauchy problem has a unique solution in a neighborhood of 
each noncharacteristic point. 


PROOF. Using a smooth change of variables we can rectify the vector field 
and convert J’ into the hyperplane zı = 0. Then in a small neighborhood of 
a noncharacteristic point we obtain the following problem: 


Ou 
TR =, Dil ait SUG Toye Ln), 


which has a unique solution. O 


2. Consider the Cauchy problem for a more general, inhomogeneous linear 
equation: 
Lalu)= f, Ulas = ug. 
The solutions of such a problem form an affine space. (The general so- 
lution of the inhomogeneous equation is the sum of the general solution of 
the homogeneous equation and a particular solution of the inhomogeneous 


equation.) 
By a smooth change of variable the problem can be brought into the form 
ðu 
{< EE TA E OE E 


This problem has a unique solution: 


Tı 


ulzi...) = uol) + f E)E. 


0 
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3. An equation that is linear with respect to the derivatives is quasilinear. 
In coordinates a first-order quasilinear equation has the form 


mlu e t t anle) = fle); (*) 

We remark that in the first two cases the field v is invariantly (indepen- 
dently of the coordinates) connected with the differential operator. How can 
a geometric object be invariantly connected with a quasilinear equation? 

Consider the space with coordinates (£1,...,Zn, Y), the space of 0-jets of 
functions of (£1,..., £n), which we denote JO(R”, R) or, more briefly, J°. 

We recall that the space of k-jets of functions of (z1,..., £n) is the space 
of Taylor polynomials of degree k. 

We note that the argument of (21,...,2n,Y,P1,---,Pn) in a first-order 
equation is the 1-jet of the function. Thus a first-order equation can be in- 
terpreted as a hypersurface in the space J!(R”, R) of 1-jets of functions. The 
space of 1-jets of real-valued functions of n variables can be identified with 
a (2n + 1)-dimensional space: J!(R”, R) ~ R*"*!. For example, for functions 
on the plane we obtain a five-dimensional space of 1-jets. 

The solution of the equation (*) can be constructed using its characteris- 
tics (curves of a special form in J°). The word “characteristic” in mathematics 
always means “invariantly connected.” For example, the characteristic poly- 
nomial of a matrix is invariantly connected with an operator and independent 
of the basis in which the matrix is formed. The characteristic subgroups of 
a group are those that are invariant with respect to the automorphisms of 
the group. The characteristic classes in topology are invariant with respect to 
suitable mappings. 

The vector field v (in the space of independent variables) is called the 
characteristic field of the linear equation L,(u) = f. 


Definition. The characteristic field of the quasilinear equation (*) is the vec- 
tor field A in J? with components (Oise is Qnad )« 


Claim. The direction of this field is characteristic. 


Indeed, let u be a solution. Its graph is a certain hypersurface in J®. 
This hypersurface is tangent to the field A, as the equation itself asserts. The 
converse is also true: if the graph of a function is tangent to the field A at 
each point, then the function is a solution. 

The method of solving a quasi-linear equation becomes clear from this. 
We draw the phase curves of the characteristic field in J°. They are called 
characteristics. If a characteristic has a point in common with the graph of a 
Solution, it lies entirely on that graph. Thus the graph is composed of char- 
acteristics. 

The Cauchy problem for a quasilinear equation is stated in analogy with 
the preceding cases. To be specific, suppose given a smooth hypersurface [7 ! 
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in x-space and an initial function ug defined on the hypersurface. The graph of 
this function is a surface I in J°, which we regard as the initial submanifold 
of Fig. 1.4. 


Fig. 1.4. Characteristics 
Tn of a quasilinear equation 
passing through the initial 


pr~! manifold ["~} 


If the characteristics are not tangent to the hypersurface I, the graph of 
the solution is composed locally of them. 

In this case two conditions are needed for a point to be noncharacteristic: 
the field A must not be tangent to ["~! and, in order for an actual graph to 
result, the vector of the field must not be vertical, that is, the component a 
must not be zero. 

Points where a = 0 are singular; at these points the differential equation 
vanishes, becoming an algebraic equation. 


Example. For the Euler equation uz; + wuz = 0 the equation of the charac- 
teristics is equivalent to Newton’s equation: t = 1, t = u, ù = 0. 


Let us now pass to the general first-order equation. 

Consider the space of 1-jets J!(R”,R). Instead of R” one can consider an 
n-dimensional manifold B”; in that case we obtain the space J!(B”, R). Let 
(x,y, p) be local coordinates in this space. 

A first-order partial differential equation is a smooth surface in Jt: 
Pee a 

For example, when n = 1 we obtain an implicit ordinary differential equa- 
tion (not solved with respect to the derivative). 

It turns out that there is a remarkable geometric structure in our space 
J', an invariantly defined distribution of 2n-dimensional hypersurfaces. For 
example, when n = 1 we obtain a field of planes in three-dimensional space. 
The structure arises only because the space is a space of 1-jets. An analogous 
structure appears in spaces of jets of higher order, where it is called a Cartan 
distribution. 

Each function in the space of k-jets has a k-graph. For a 0-jet this is the 
usual graph — the set of 0-jets of the function: Tu = {jeu zE R”} = 
{(x,y): y = u(x) }. In the case of 1-jets a point of the 1-graph consists of the 
argument, the value of the function, and the values of the first-order partial 


derivatives: {jiu : z € R”} = {(z,y,p): y = u(x), p = S*} (see Fig. 1.5 
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for n = 1). We remark that the 1-graph is a section of the bundle over the 


domain of definition. 
E -Menn p is the slope 
l 


Fig. 1.5. A point 
of the space of 1-jets 


Remark. The surface of the 1-graph is diffeomorphic to the domain of defi- 
nition of the function, x is the n-dimensional coordinate on this surface. The 
smoothness of this surface is 1 less than the smoothness of the function, but 
smoothness is preserved for an infinitely differentiable or analytic function. 


Consider the tangent plane to the 1-graph. This is an n-dimensional plane 
in a (2n + 1)-dimensional space. 


Theorem 3. All tangent planes of all 1-graphs at a given point lie in the same 
hyperplane. 


PROOF. Along any tangent plane we have dy = )> = d; = $ pidt or 
dy = pdz. Since p is fixed at a given point of the space of 1-jets, we obtain 
an equation for the components of the tangent vector that determines the 
hyperplane. Thus the tangent plane to any 1-graph lies in this hyperplane. O 


For example, when n = 1 the equation dy = pdz defines a vertical plane 
in space with the coordinates zx, y, p. The tangents to the l-graphs are the 
non-vertical lines lying in this plane (Fig. 1.6). 


Fig. 1.6. The contact plane y 
in the space of 1-jets 


In this case one can see that the hypersurface itself is the closure of the 
union of the tangents to all 1-graphs passing through the given point. 


Problem 3. Prove that this is true for any dimension. 
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Corollary. The field of hyperplanes dy = pdx constructed is invariantly de- 
termined, that is, in other coordinates it is given by the equation dy = pdz. 


Definition. The field of hyperplanes in J! is a Cartan distribution or a stan- 
dard contact structure. 


Problem 4. What are the dimensions of the integral manifolds for the field 
of contact planes? (A manifold is called integral if the tangent plane at each 
of its points is a subspace of the contact plane.) 

A 1-graph is always an integral manifold, that is, a 1-dimensional integral 
manifold exists. But are there integral manifolds of higher dimension? 


ANSWER. A field of contact planes has no integral manifolds of dimension 
greater than half the dimension of the contact plane. 


Definition. An integral submanifold of a field of contact planes whose di- 
mension is maximal (that is, equal to half the dimension of the contact plane) 
is a Legendre manifold. 


For example, 1-graphs are Legendre manifolds. 

Let us now return to the differential equation. 

An equation is a 2n-dimensional submanifold T?” in J?. 

At each regular point of this surface one can distinguish a characteristic 
direction determined by the surface and the contact structure. We shall con- 
struct the characteristics (integral curves of this direction field), and then use 
them to form the integral manifolds. 

At a point of the surface T?” we consider the intersection of the tangent 
and contact planes. These planes either coincide or have a (2n—1)-dimensional 
intersection. In the first case the point is singular, in the second case it is 
regular. 

We remark that for a surface I’ in general position the singular points 
are isolated. Indeed there are 2n coordinates on J’. Consider the normal to 
the tangent plane and the normal to the contact plane. A point is singular 
if these normals have the same direction. This means that 2n functions of 
2n variables vanish simultaneously. In general position this can happen only 
at isolated points. 

Thus, at regular points there are (2n — 1)-dimensional intersections of the 
tangent and contact planes. These are lines when n = 1, but not for n > 1. 
How can we distinguish a one-dimensional direction? 

In local coordinates the contact field is given by the zeros of the 1-form 
a = dy — pdz, and this form can be multiplied by any nonzero function 
without changing the field (contact structure). 

The 2-form w? = da, which is the exterior differential of the form a, is 
no longer defined invariantly by the contact structure. However, the following 
proposition holds. 


Proposition 1. The form 7 is defined tnvariantly up to multiplication 
by a nonzero number at each point. 
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PROOF. Let à = fa. Then da = df Aa+ fda. 
dé|,_, =fda|,_o. 


that is, ©? differs from w? by multiplication by a number at each point (the 
conformal type of the form w? is said to be invariantly defined). We remark 
that & = 0 wherevera=0. O 


Proposition 2. w| o is a symplectic structure. 


We recall that a symplectic structure is a nondegenerate skew-symmetric 
bilinear form in an even-dimensional space. 
Nondegeneracy of the form w means that YVE #0 3n w(E,n) £0. 


PROOF. In local coordinates our form is written as: da = — >> dp; A da;, 
where p; and x; are coordinates in the plane œ = 0. O 


Exercise. Write out the matrix of the form ` dp; A dz; and verify that this 
form is nondegenerate. 


This form is called the skew-scalar product. Let us clarify its geometric 
meaning. 

Let n = 1. Then w = dz A^Adp. The value of this form on a pair of vectors is 
the oriented area of the parallelogram spanned by the vectors, as in Fig. 1.7. 
In the case of higher dimension w(£,ņ) is the sum of the oriented areas of 
the projections of the parallelogram with sides € and 7 on the plane with 
coordinates (2;, pi). 


Fig. 1.7. A symplectic 
structure 


We recall that in Euclidean space there is a concept of orthogonal comple- 
ment. In n-dimensional space the orthogonal complement to a k-dimensional 
subspace is an (n — k)-dimensional subspace. The proof of this fact uses only 
the bilinearity and nondegeneracy of the scalar product, not its symmetry, 
so that the same is true in the case of the skew-scalar product. 

Thus the skew-orthogonal complement to a (2n — 1)-dimensional plane in 
2n-dimensional space is a line. But, in contrast to the Euclidean case, the line 
lies in this plane! 


10 1. The General Theory for One First-Order Equation 


Lemma. The skew-orthogonal complement to a hyperplane in symplectic 
space is a line lying in the hyperplane. 


PROOF. Let the line be spanned by the vector £. Its skew-orthogonal comple- 
ment is the hyperplane {7 : w(£,7) = 0}. The vector £ lies in this hyperplane 
since w(€, £) = wl.) =0. 0 


Definition. A characteristic direction in a contact plane is the skew-orthogo- 
nal complement to the intersection of the contact plane and the plane tangent 
to I’ at a regular point. 


This skew-orthogonal complement is a line. Thus there is an invariant con- 
tact structure and an invariant relation of skew-orthogonality in each contact 
plane. Consequently a characteristic direction is invariantly distinguished at 
each regular point (see Fig. 1.8). 


z is a regular point 
contact plane 
tangent plane 


the (2n — 1)-dimensional intersection 
contains the one-dimensional 
skew-orthogonal complement 
equation 


singular point 


Fig. 1.8. A characteristic direction for a general first-order equation 


The integral curves of this direction field are called characteristics. 


Problem 5. Compute the characteristic direction field in coordinates z, y, p, 
that is, determine it by a system of differential equations t = ?, y= ?, p = ?. 
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Lecture 2 


The General Theory for One First-Order 
Equation (Continued) 


We are considering a general first-order partial differential equation 
F(ayp)-=]-0, where = (hse Ons P = Digests Pale De = pe, and 
y = u(x) is an unknown function. The equation determines a 2n-dimensional 
hypersurface V?” in the space J! of 1-jets of functions of (£1,.-., £n). Each 
function has a l-graph in J}; it is a solution of the equation if its 1-graph is 
a submanifold in V2". 

At each point of J! there is a contact plane K?” defined in local coordinates 
by the equation dy = pdz; geometrically it is the closure of the union of the 
tangent planes to all the 1-graphs passing through the point. 

At points of the surface V?” the tangent plane to V2” intersects the contact 
plane. If the intersection is (2n — 1)-dimensional, the point is regular; other- 
wise it is singular. In general position singular points are isolated. At regular 
points z we obtain a distribution of (2n—1)-dimensional planes (T,V?7")NK2”, 
which are subspaces of the contact planes K2”. 

Each contact plane is a symplectic space, the symplectic structure being 
defined by the differential 2-form w? = da! gen» Where a = dy — pdr. 


n 
In coordinates we have w? = dx A dp := Y da; A dp;; (x, p) constitute 
i=l 
coordinates in the contact plane. 


Lemma. The form w? is nondegenerate. 


PROOF. Its matrix in the coordinates z1, p1, £2, p2,... has the form 
0 1 0 
—1 0 
0 1 g 
-1 0 
0 


Hence the form w? defines a skew-scalar product in K?", and the skew- 
orthogonal complement to each subspace in K?” is defined and has the com- 
plementary dimension. In particular, the skew-orthogonal complement to the 
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(2n — 1)-dimensional intersection (T, V?”) N K?” is one-dimensional. This di- 
rection is characteristic for the equation. 


Problem 1. Prove that the characteristic line lies in (T,V°") N K2”. 


Let us compute the components of the characteristic vector explicitly in 
terms of the equation. Let (x, y, p) be the coordinates in T, J’. The character- 
istic vector must be tangent to V2": by differentiating the equation, we obtain 
the first condition 

Fzt + Fyy + Fpp =0. (2.1) 
Moreover, the vector must lie in the contact plane; hence we obtain the second 
condition 

y= pi. (2.2) 

Eliminating y from these two equations, we obtain the equation of the 
(2n — 1)-dimensional intersection K?"~' of the tangent plane and the con- 
tact plane: 

(Fz + Fyp)t+ Fpp =0. (2.3) 
This equation is written in (x, p)-coordinates, which we can regard as coordi- 
nates in the contact plane, since they project one-to-one onto the (x, p)-hy- 
perplane in J! (see Fig. 2.1). 


K2n-) 


K?” 


Ki2n-l 


Fig. 2.1. x and p are coordinates in the contact plane 


It remains for us to find the skew-orthogonal complement to K2"~!. The- 
oretically this problem could be solved algorithmically by writing a suitable 
system of linear equations, but we shall take advantage of a useful observation. 

As an example let us consider the case n = 1. Then the value of the skew- 
scalar product dz Adp on a pair of tangent vectors (4, p), (x', p’) is simply the 
determinant of their coordinates: 


dz app) (a ep )) = a H = tp — r'p. 


It turns out that the same is true in general if x and p are interpreted as 
vectors and their product as the scalar product. 
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Problem 2. Prove that 
w’ ((t, p), (£, p) = tp! — x'p. (2.4) 


But let us examine (2.3) carefully. The equation itself has the form (2.4), 
that is, it expresses the fact that the skew-scalar product equals 0. Thus we can 
write out the skew-scalar complement to K2"~! directly from this equation: 
p' =F,+pF,, £ = —F,. (To be sure, it is the oppositely directed vector that 
is traditionally called the characteristic vector.) As a result, we have proved 
the following theorem. 


Theorem 1. The direction field of the vector field given by x’ = Fp, P = 
—(F, + pF,), y = pF, is characteristic (the last component being obtained 


from (2.2). 


We remark that the field is determined by the function F itself, not just 
by its zeroth level surface, so that it is defined in the entire space Jt. 


Example 1. Suppose F is independent of y. The corresponding equation is 
called a Hamilton-Jacobi equation. The traditional notation for it is F = 
H(z, p). A Hamilton-Jacobi equation is an equation of the form H(z, gu) = 0. 
The equations of its characteristics have the form 

, . OH , OH jC 

OC apy EE oa a ae te (2.5) 

The first two equations are the usual canonical Hamilton-Jacobi equations. 

(And this determines the choice of sign for the characteristic field indicated 
above: for the Hamiltonian of a free particle of unit mass H = p?/2; with this 
choice of sign we obtain x’ = p, that is, the momentum is numerically equal 
to the velocity.) 


By projecting the system (2.5) onto the subspace with coordinates (x, p), 
we obtain a separate system of equations in this space. Thus we have carried 
out a factorization or splitting of the original equation (2.5). This vector 
field in the space (x,p) can also be invariantly defined, without resorting to 
coordinates. The surface of the original equation in J! was cylindrical, not 
depending on y. Such a surface projects very well onto the (x, p)-space. The 
Hamiltonian field is defined not only on the zeroth level surface of H, but 
everywhere. The function H usually has the physical meaning of energy. 

Suppose, for example, H = (p? — 1)/2. Then the equation has the form 
p?=lor Gaye = 1, that is, it is the eikonal equation of geometric optics. 

The system of equations of characteristics z’ = p, p’ = 0 describes the 
motion of particles along straight rays at constant speed (see Fig. 2.2). 

If a Riemannian metric is defined on the manifold, the characteristics of 
the eikonal equation turn out to be the geodesics of this metric. 
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Fig. 2.2. The characteristics 
of the eikonal equation 


Example 2. Consider the Euler equation us + wuz = 0. When n = 2, the 
space of 1-jets J! is five-dimensional, and F(t, £,Y, Pt, Pr) = Pt + yPr- Let us 
write the system of equations of the characteristics: 


fan Po 
Pr = Pr » tT =y, 
The first two equations are the so-called conjugate equations of evolution 
of the derivatives, which we shall not study at the moment. Let us take up 
the remaining equations. The variable t represents time, and we write the 
remaining equations as dr/dt = y, dy/dt = 0, whence d?x/dt? = 0, that is, 
the equation of the characteristics has turned out to be Newton’s equation 
(see the previous lecture). 


We now apply the characteristics to solve the equations in the general case. 


Theorem 2. Let I” be the 1-graph of a solution, and suppose a characteristic 
passes through a point of this graph. Then the characteristic lies entirely in 
I”, that is, the graph of a solution can be fibered into characteristics. 


PROOF. Consider a regular point z, and let £ be the vector of the character- 
istic direction at the point z. The following are obvious: 


1. 7,1" C T,V?”, since I is a submanifold in V2”. 

2. T,I'” C K?” by definition of K?” , and therefore T,” C (T,V2")N(K2"). 
We remark that the space T, I” itself is n-dimensional. 

3. All vectors in T, I” are pairwise skew-orthogonal. Indeed, al pa = 0, since 
T” is a l-graph; therefore da| „„ = 0, that is, w?| „„ = 0. 


rr 


Exercise. Verify that in coordinate notation the vanishing of the form w? | pi 
follows from the symmetry of mixed partial derivatives. 


Definition. A subspace of symplectic space all of whose vectors are pairwise 
skew-orthogonal is isotropic (“crazy” according to Lie). 


For example, every line is isotropic on the symplectic plane. 


Lemma. In 2n-dimensional symplectic space the dimension of an isotropic 
subspace is at most n. 
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(Isotropic subspaces of maximal dimension exist and are called Lagrangian 
subspaces. An example of a Lagrangian subspace is a line in the symplectic 
plane.) 


PROOF. Suppose an isotropic subspace is m-dimensional. Being isotropic, it 
is contained in its skew-orthogonal complement, and hence the dimension of 
the skew-orthogonal complement is at least m. We then have 2n — m > m, 
so that m < n. The lemma is now proved. O 


We now return to the characteristic vector € at the point z. 
Claim. The vector € lies in T,I"™. 


Suppose the contrary. Then, being characteristic, € is skew-orthogonal to 
T,V2" N K2". By virtue of item 2 above, € is skew-orthogonal to T, I”. But 
then the linear span of € and T, I" is (n + 1)-dimensional and is isotropic by 
virtue of item 3, all of which contradicts the lemma. 


Thus, at each point the characteristic direction is tangent to the 1-graph 
of the solution, and so the entire characteristic lies on the graph. The theorem 
is now proved. U 


A recipe for constructing the 1-graphs of solutions now follows. (We remark 
that the existence of even one solution has not yet been proved.) One must 
take an (n —1)-dimensional submanifold of the surface V?” that is not tangent 
to a characteristic and draw the characteristics through its point, thereby ob- 
taining the 1-graph of a solution locally. To obtain a real graph it is necessary 
that the plane spanned by the tangent plane to an isotropic submanifold and 
the characteristic direction project one-to-one onto the x-space. 

There is a standard method for constructing initial isotropic manifolds. 
Consider an (n—1)-dimensional surface y”~! in x-space, and suppose an initial 
function uo is defined on it. Substitute these data into the equation and regard 
it as an equation with respect to p. We have already defined the derivatives 
along the directions tangent to y”~!. The derivative in a transversal direction 
can be found from the equation using the implicit function theorem if F is 
not tangent to the initial hypersurface y”~! at the point in question. Thus we 
obtain an (n — 1)-dimensional isotropic initial submanifold; data of this type 
are called Cauchy data. For details see [1, §8, pars. J and L]. 


Theorem 3. Consider an initial (n — 1)-dimensional submanifold V"-! c 
V2" that is not tangent to a characteristic. Let the n-dimensional plane 
spanned by its tangent plane and the characteristic direction project isomor- 
phically to x-space. Draw the characteristics through its points. Then we obtain 
the 1-graph of a solution locally (see Fig. 2.3). 
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Fig. 2.3. Construction of the l-graph 
of a solution 


yen 
yr-l 


PROOF. The submanifold V”~! is isotropic: da rye = 0. Drawing the 
characteristics through all of its points in a neighborhood of a regular point 
z, we obtain an n-dimensional submanifold T” in V2". We shall prove that 
da pn = 0. Let v be the vector field of characteristic directions. By a homotopy 
formula, Lya = ida + d(iya), where i, da(€é) := da(v,&) and ipa = a(v). 
It is obvious that iœ = 0 and i, da(é) = da(v,€) = 0 for every vector £ 
that is tangent to V2", by definition of a characteristic direction. Thus, the 
derivative of the form in the direction of the field v on V?” equals 0. Therefore 
the phase flow of the characteristic field maps the form into itself. This means 
that the value of the form on any vector equals the value of the same form on 
the vector translated back to the initial manifold. 

It follows that da is also mapped into itself by the field v. At the origin 
the space T, I” is isotropic. (V"~! is isotropic and v is characteristic.) This 
property carries over to any point z’ of the manifold IT” by the phase flow of 
the field v: the space T; I” is isotropic. Thus al pa = 0. Hence I” is an integral 
submanifold of the contact structure. Thus I” is a Legendre submanifold. And 
by construction J” is a submanifold of the equation. It remains to be proved 
that it is locally a 1-graph. That follows from the following general theorem. 


Theorem 4. Suppose an n-dimensional surface I” in the space of 1-jets is a 
Legendre manifold (an integral manifold for the contact structure) and locally 
projects in a one-to-one manner onto the x-space. Then it is locally a 1-graph 
(see Fig. 2.4). 


Fig. 2.4. The integral 
surface as a l-graph 
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PROOF. We write the condition for a Legendre manifold a| pa = 0 (the surface 
is an integral surface for the contact structure): 


n Of n E _ of 
(2 Ox; i? giite =o] => (s = ax) 


since x; are local coordinates on J”. But this means that I” is a 1-graph. O 


Example. Let us solve the Cauchy problem for the Euler equation u;+uuz = 
0. The initial curve in the (z,t)-plane is usually chosen to be the line ¢ = 0. 
We prescribe an initial function y = uo(zx), as in Fig. 2.5. 


Fig. 2.5. The initial condition 
for the Euler equation 


The equations of the characteristics are t = 1, y’ = 0, x’ = y, and the 
transversality conditions are satisfied. We remark that for a quasi-linear equa- 
tion everything can be considered in the space J°, since the characteristics of 
a quasilinear equation in J’ are the projections of the “true” characteristics 
from J!, as shown in Fig. 2.6. 


J’, dim =5 


1-graph 


Fig. 2.6. The characteristics 
of a quasilinear equation in J° are 
projections of the characteristics of J! 


“forgetting p” 


J°, dim = 3 


0-graph 


Thus we shall study the 0-graph in J°. On the (x, y)-plane we draw a 
Series of sections of the graphs for t = const, as in Fig. 2.7. The values of the 
Solution at subsequent times are obtained by translation of the initial values 
along the characteristics. 
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Fig. 2.7. The graphs 
of the solution of the Euler 
equation at successive times 


We see that from some time on the curve ceases to be a graph. The integral 
surface of the equation ceases to project in a one-to-one manner onto the 
(x,t)-plane. The curve of critical values for the projection has a cusp, as in 
Fig. 2.8. 


y graph non-graph 
x 
place where 
the tangent 
is vertical t 
(a) (b) 


Fig. 2.8. (a) The integral surface ceasing to be a graph; (b) the curve of critical 
values of the projection 


Physically the Euler equation describes the evolution of the velocity field 
of noninteracting particles. This model is good, for example, for describing 
star streams; the loss of uniqueness of the solution can be interpreted as the 
free passage of different streams through each other. 

On the other hand, for large densities the particles begin to collide, and 
after a certain amount of time the Euler equation ceases to hold. It gets 
replaced by another equation that takes account of interaction, for example, 
Bürgers’ equation uz + wuz = EUzz. For small £ its solutions approximate the 
solutions of the Euler equation up to the critical time, while for large values of 
time they are of shock wave type: the graph of u is nearly vertical in a small 
neighborhood (of order £) of a moving point, as shown in Fig. 2.9. 

To the right and left of this point the solutions are also close to solutions 
of Euler’s equation. It is remarkable that Burgers’ equation can be solved 
explicitly — it reduces to the heat equation, which we shall study below. 
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Fig. 2.9. Solutions of shock wave type 


Remark. The 1-graph of a solution of a first-order equation consists of the 
characteristics of the hypersurface that defines the equation in the space of 
1-jets, as in Fig. 2.10. 


F(z,y,p) = 0 y = u(x) 


p = ðu/ðx 


Fig. 2.10. The l-graph 
of a solution in the space 
of l-jets 


Thus we have reduced the equation of the field to the equations of motion of 
particles. But, as Jacobi remarked, one can apply the theorem in the opposite 
direction: the equations of motion of particles can be reduced to the equation 
of the field. This is useful, since the systems of equations for particle motion 
cannot as a rule be solved explicitly. 

If we find one solution of the equation of the field (a partial differential 
equation), we shall know a Legendre submanifold composed of characteris- 
tics. If we find another solution of the field equation, we shall have imposed 
another condition on the characteristics, and the integral manifold of inter- 
section decreases in dimension by 1. If we succeed in finding sufficiently many 
solutions, we shall have reduced the dimension of the integral manifold. If a 
sufficiently large (finite-parameter) family of solutions of a partial differential 
equation are known, all the characteristics can be obtained by intersecting 
the corresponding Legendre manifolds until we isolate the (one-dimensional) 
characteristic itself. 

This is Jacobi’s method of integrating ordinary differential equations by 
representing them as the equations of the characteristics of some partial dif- 
ferential equation. 

Jacobi applied this method to the solution of Hamiltonian systems — 
the systems of equations of characteristics for the Hamilton-Jacobi equation 
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(in which F does not explicitly depend on y). One result was, for example, 
the famous theorem of Liouville on the integration of Hamiltonian systems 
for which a complete set of integrals in involution is known. 

Astonishingly, this method could be “lifted” to integrate partial differen- 
tial equations as well! To do this it is necessary to regard them as infinite- 
dimensional Hamiltonian systems. In this way the well-known sine-Gordon 
equation and the Korteweg-de Vries equation were integrated (see (2, Ap- 
pendix 16]). 

But it follows from our remark that the Jacobi method can be applied to 
a larger class of systems than Hamiltonian systems, namely to those that are 
the equations of the characteristics of some partial differential equation. To be 
sure, it is difficult to find out whether a specific vector field is characteristic 
for some equation. The corresponding theory seems to be undeveloped as yet; 
at least I do not know any applications of it to physically interesting problems. 
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Lecture 3 


Huygens’ Principle in the Theory of Wave 
Propagation 


We consider space-time M”+1 = B” x R, where B” is “physical space,” and 
R is the “time axis.” (Most of what follows carries over word for word to the 
case of the bundle M”"+! — R over the time axis, and a great deal of it also 
carries over to the case when the role of time is played by a fibration into 
“isochrones” or even the zeros of a nonclosed differential form “dt” on M.) 

In geometric optics (and also in the calculus of variations and optimal 
control theory) a cone of possible velocities of motion is prescribed at each 
point of a manifold M. 


Example. Let B be a Riemannian manifold. The graphs of all possible “mo- 
tions” q : R — B with unit velocity at each point is tangent to the quadratic 
cone (which defines the Riemannian metric) ||dg|| = |dt|. Thus, M is endowed 
with a field of Lorentz quadratic cones. 


The tangent planes to the cone of possible directions at the point m that 
forms its vertex belong to the projectivization of the tangent space to M at 
that point: PTX M = RP”. 

It is often useful to endow these tangent planes with co-orientations (in- 
dicating the direction of propagation of a disturbance). The oriented tangent 
planes of the cone of possible velocities belong to the spherization ST% M of 
the cotangent space of space-time M. This spherization is diffeomorphic to 
the sphere S”. 

The non-co-oriented (resp. co-oriented) hyperplanes in the tangent space of 
a manifold are called the contact elements (resp. co-oriented contact elements) 
of the manifold. They form the manifold of the bundle of contact elements 
(resp. co-oriented contact elements) PT*M — M (resp. ST*M — M) with 
fiber RP” (resp. S”), where n = dim M — 1. 

Thus the field of cones of possible velocities defines a hypersurface in the 
manifold of (possibly co-oriented) contact elements of space-time. In geometric 
optics this hypersurface is called the Fresnel hypersurface. It is a fundamental 
geometric object of geometric optics, calculus of variations and optimal control 
theory. 
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Remark. In the theory of wave propagation the vectors of the space dual to 
the velocity space are often called retardations. 


The hypersurface constructed above can be described as the field of cones 
in the space T*M of retardations on space-time. 

On the manifolds of contact elements of any manifold (co-oriented or not) 
there is a remarkable geometric structure — the contact structure. That is the 
field of hyperplanes in the tangent spaces of the manifold of contact elements, 
which is invariantly determined by the bundle of the manifold of contact 
elements over the original manifold M. 

Such a field exists and is unique; it is called the tautological field and is 
defined by the following construction. Each point of the space of the bundle of 
contact elements on M is a hypersurface in the tangent space of the original 
manifold M. The preimage of this plane under projection of the space of the 
bundle on M is the hyperplane of the tautological field at the initial point. 


Remark. The tautological contact structure gives a condition on the velocity 
of motion of a contact element. This condition is called the skate condition. 
The meaning of the condition is that a skate (the contact element on the plane 
of a rink) can rotate freely in place and can move in the direction determined 
by itself, but resists any attempt to move it in transversal directions. 


Let us consider the integral manifolds of the tautological contact structure. 


Example 1. For any hypersurface in M the set of contact elements tangent to 
it forms the integral manifold of the tautological contact structure in PT* M 
(or in ST* M, if they are co-oriented). The dimension of this integral manifold 
is slightly less (by 1/2) than half the dimension of the manifold of contact 
elements. 


Example 2. The set of contact elements at a point of the manifold M 
forms the integral manifold of the tautological contact structure in PT* M 
(or ST*M). The dimension of this integral manifold is the same as the di- 
mension of the integral manifold of Example 1. 


Example 3. For any submanifold (of any dimension) in M the set of contact 
elements on M tangent to it forms an integral submanifold of thestautological 
contact structure in PT*M (or ST* M) of the same dimension as in Examples 
1 and 2. 


Problem 1. Prove that the manifold of contact elements on M”t! with 
its tautological structure has no integral submanifolds of dimension greater 
than n. 


Problem 2. Do there exist smooth integral submanifolds of dimension n for 
the tautological contact structure on the manifold of contact elements in M”+! 
not obtained by the construction of Example 3? 
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HINT. Consider the contact elements that are tangent to the semicubical 
parabola x? = y? on the plane (n = 1). 


Definition. Integral manifolds of maximal dimension (which is n for a con- 
tact manifold of dimension 2n + 1) are called Legendre manifolds. 


Theorem 1 (the theory of support functions). The manifold of 1-jets of 
functions from the sphere S”~* into R” (with its natural contact structure) is 
contact diffeomorphic to the manifold of co-oriented contact elements in R” 
(with its tautological contact structure). 


PROOF. We denote by q € S”! C R” a point of the unit sphere in Euclidean 
space. The tangent vectors to the sphere at the point q can be regarded as the 
vectors in the space orthogonal to q. They can also be regarded as cotangent 
vectors, since the Euclidean structure identifies a vector p with the linear 
functional (p,-). Thus the 1-jet of a function f on S”~-! is given by a triple 
(ge S}, p=df|, € R”, z= f(q)). 

The required diffeomorphism is defined by the formula 


Q=qtpz, P=q, 


where Q is the point of application of the contact element of R” co-oriented 
by the vector P normal to it. O 


Problem. Prove that: 1) this is indeed a diffeomorphism J1(S"~1,R) ~% 
ST*R” (x~ S} x R”); 2) it maps the natural contact structure in J? into 
the tautological contact structure ST*. 


HINT. Consider what the Legendre manifolds map to. 


Remark. The distance from the origin to the tangent plane orthogonal to q 
of a convex hypersurface in R” is called the value of the support function of 
the given hypersurface at the point q € S"~!. The support function is defined 
on the sphere of (outward) normals and defines a hypersurface. 

A function defined on the sphere determines a hypersurface in the space 
(the envelope of the family of hypersurfaces corresponding to the function). 
This hypersurface, however, may have singularities. 


Problem. Study the curves defined by the functions x? + 2y? + t (where t is 
a parameter) on the circle x? + y? = 1. 


Let us return to the Fresnel hypersurface of retardations of fronts lying in 
the manifold ST*M”t! of co-oriented contact elements of space-time. 


Definition. The rays (on this hypersurface) are its characteristics, that is, 
the integral curves of the field of characteristic directions of the hypersurface. 
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We recall that a characteristic direction of a hypersurface at a point of a 
manifold with a contact structure (defined as the field of zeros of a differential 
1-form @) is the skew-orthogonal complement (in the sense of the symplectic 
structure da on the plane œ = 0 at each point) to the intersection of the 
tangent hypersurface of the hypersurface with the hyperplane of the contact 
structure. 


Example. Let (q1,...,@n) be the coordinates in “physical space” B”, let 
qo = t be time, and let (po, P1,- --, Pn) be the corresponding components of the 
momentum. The tautological contact structure is given locally by the 1-form 
a = po dgo + pı dqi +--+ + Pn dgn, where [po : pi : -++ : Pn] are homogeneous 
coordinates in RP” and qo = t. In affine coordinates, for which pp = —1, we 
have a = —dt + pı dq, +--: + Pndqn. The retardation hypersurface is given 
by the equation 
Po = H (pis. +5 Pai qis: --30n;t) , 


where H is a homogeneous function of degree 1 with respect to the variables 
p;. The equation of the characteristics has the form 

dpi OH dq; ƏH 
dt ðq’ dt Op, 


Thus the rays are defined by the Hamilton equations, the Hamiltonian 
being given by the retardation hypersurface. 


Let us now consider the geometry of wave propagation in a medium in 
which the local velocity of propagation of disturbances is defined by this re- 
tardation hypersurface in the space of contact elements of space-time. 

A typical example is the family of equidistants of a submanifold of a Rie- 
mannian manifold. The retardation hypersurface defines a second-order hy- 
persurface in each projectivized cotangent space of space-time. The family 
of equidistants (hypersurfaces lying at distance t from a given initial hyper- 
surface in physical space) can be regarded as a single hypersurface in space- 
time (whose sections by different isochrones t = const give the entire family 
of equidistants). Similarly one can define the large front that describes the 
propagation of fronts of perturbations by means of a single hypersurface in 
space-time with a more general retardation hypersurface in the manifold of 
contact elements of space-time. 

Let us consider the contact elements of space-time that are tangent to the 
large front. They all belong to the retardation hypersurface (that is the content 
of the local law of propagation of disturbances). The following statement is 
obvious. 


Proposition. The contact elements tangent to a large front form a Legendre 
submanifold of the tautological contact structure of the space of contact ele- 
ments of space-time that projects to the large front and lies in the retardation 
hypersurface. 
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Let us now consider the initial condition given by the instantaneous front 
of a disturbance at time t = 0. This front determines a submanifold of codi- 
mension 1 on the large front, and this submanifold is also integral for the 
tautological contact structure of the manifold of contact elements of space- 
“time. It is not a Legendre manifold, however, since its dimension is one less 
than that of the Legendre manifolds. 

This initial integral submanifold consists of contact elements of space-time 
that are tangent to the initial front and in addition belong to the retardation 
hypersurface. It is this last condition that renders it possible (under suitable 
nondegeneracy conditions) to choose the contact element of space-time that 
is tangent to the cone of possible velocities out of the one-parameter family 
of contact elements of space-time containing the tangent space to the initial 
front at the given point. 

The basis of the entire theory of propagation of disturbances is provided 
by the following simple and general fact of contact geometry (discovered in 
essence by Huygens and hence deserving the name “Huygens’ principle” ). 


Theorem 2. A Legendre submanifold of a hypersurface in contact space con- 
tains, along with each point, the entire characteristic of that surface passing 
through the point. 


It is assumed here that the tangent plane to the hypersurface nowhere 
coincides with the contact plane a = 0. 


PROOF. Along a Legendre manifold we have œ = 0, so that da = 0 in the 
tangent plane to the manifold. If a characteristic direction vector € did not lie 
in this plane, it would be skew-orthogonal to it in the sense of the form da, and, 
adjoining € to the tangent plane of the Legendre manifold, we would obtain 
a plane skew-orthogonal to itself in the symplectic vector space a = 0 having 
dimension greater than half the dimension of the space. This is impossible, and 
therefore the characteristic directions are tangent to the Legendre manifold, 
as asserted. UJ 


Corollary. The Legendre manifold corresponding to a large front is obtained 
from the initial condition by the following construction. The initial front is 
lifted to the integral manifold of the space of contact elements of space-time 
lying in the retardation hypersurface. We then consider the characteristics 
of this hypersurface passing through points of the integral manifold just con- 
structed. They form a Legendre manifold. The projection of this manifold on 
Space-time is a large front. The sections of the large front by the isochrones 
t = const are the instantaneous fronts. 


This constitutes a description of the propagation of disturbances using 
waves (fronts) and rays (characteristics). 

Since the characteristics are determined by ordinary differential equations, 
each of them is determined by any one of its points. They will be the same 
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under completely different initial conditions for which only the tangent hyper- 
planes of the initial fronts coincide at one point at the initial time. Physically 
this means that infinitely small pieces of a wave front propagate (along the 
characteristics) independently of one another. That makes it possible to in- 
voke particles (whose motion is determined by ordinary differential equations, 
namely the Hamilton equations for the characteristics of the retardation hy- 
persurfaces) to describe the propagation of the front (waves, which are deter- 
mined by partial differential equations). 

In particular, the initial condition can even be taken to be a point front 
(corresponding to the Legendre manifold forming a fiber of the bundle of 
contact elements). The theorem is applicable in this case also (although it 
appears to be degenerate from the point of view of the corollary just given). 

A “spherical” wave corresponds to a point initial front. Comparing the 
propagation of disturbances from an arbitrary initial front with those from 
one of its component point sources, we see that the Legendre manifolds cor- 
responding to it contain a common characteristic. It follows that the instan- 
taneous front is subsequently tangent to the spherical front of the original 
point source at a point on the ray emanating from the source in the direction 
determined by the direction of the initial front. 

It follows that at time t the instantaneous front is the envelope of the family 
of spherical fronts of point sources belonging to the original front. This was 
the original formulation of Huygens’ principle (known under different names 
in different branches of mathematics: Pontryagin’s principle, the canonical 
Hamilton equations, and so forth). 

Incidentally, the equations of characteristics are precisely the infinitesimal 
version of Huygens’ assertion about the envelope, corresponding to infinitely 
small values of t and consequently to infinitely small spherical fronts. 

Naturally it was assumed in the geometric theory described above that 
the necessary determinants are nonzero, so that the implicit function theorem 
could be applied. For example, the field of Lorentz cones in space-time cannot 
be tangent to the isochrones t = const at any point, and so on. It is not 
difficult to verify that these conditions hold automatically in such examples 
as the problems of Riemannian geometry, even with a time-dependent metric. 
In other problems, for example, in optimal control theory, the situation often 
turns out to be more complicated and the need to investigate singularities 
arising when the determinant vanishes constitutes a major difficulty. 

The fact that the characteristics passing through points of the initial inte- 
gral manifold form a Legendre manifold needs the proof. We shall not go into 
detail on this point, since it is exactly the same as the proof of Theorem 3 
given in the preceding lecture. 


Lecture 4 
The Vibrating String (d’Alembert’s Method) 


We consider a string, which can be pictured as the actual string of a musical 
instrument, attached to a sounding board, as in Fig. 4.1. 


u(x, t) 


Let u(x,t) be the displacement of the string at the point x at time ¢ in 
the orthogonal direction from its equilibrium position. It is known that if the 
string is homogeneous and the displacements are small compared with the 
length of the string, the function u(x,t) satisfies the following linear second- 
order partial differential equation: 


Fig. 4.1. The vibrating string 


4 


Utt = a? Urz ; (4.1) 


We shall not discuss the derivation of this equation in the present lecture. 
Here a is a constant having the physical dimension of velocity, as will be 
shown below. 


4.1. The General Solution 


Problem 1. Show that the change of variables 
€=x-at, n=x-+at 
reduces the equation (4.1) to the form 
Uen Oss (4.2) 


(When a = 1, this change of variables is the composition of a rotation and a 
dilation in the plane of the independent variables.) 


28 4. The Vibrating String (d’Alembert’s Method) 


The equation (4.2) is easy to solve. We write it in the form (ug), = 0, 
which means that the function uç is independent of 7, that is, constant along 
vertical lines on the (€,7)-plane, as shown in Fig. 4.2. 


g 


Fig. 4.2. The characteristics of the vibrating string equation 


Hence ue = f(E), that is, the differential equation has been reduced to a 
simple equation of the form F’ = f, so that u = f f(€) d€ = F(€) + G(n); we 
remark that the constant of integration in general depends on the line along 
which the integration occurs, that is, it depends on 7. 

If we seek a solution in the class of functions that are continuous at every 
stage of the computation, we obtain a general solution of the form f(€)+g(n), 
where the function f is smooth and g is continuous. We have thus arrived at a 
“solution” that in general might have no derivatives! Moreover, if we solve in 
the opposite order, we obtain an asymmetry in regard to smoothness. How- 
ever, equality of mixed partial derivatives holds only in the case of sufficiently 
smooth functions. This observation serves as motivation for the introduction 
of generalized functions (distributions), a class in which one can consider a 
solution of the form f(€) + g(7) with nonsmooth f and g. 

In the present lecture we assume that there exists a solution in the class 
of sufficiently smooth functions, so that f and g are sufficiently smooth. 

Thus, u(x,t) = f(x — at) + g(x + at) is the general solution of the equa- 
tion (4.2). 


4.2. Boundary-Value Problems and the Cauchy Problem 


The Cauchy problem reveals both the similarities and the differences between 
the theories of ordinary and partial differential equations. In the theory of 
ordinary differential equations the phase space is finite-dimensional, whereas 
in this book we shall be dealing with an infinite-dimensional phase space. 

The Cauchy problem for the vibrating string consists of the equation (4.1) 
with initial conditions 


ulo = (2) > ulizo = ¥(2) (4.3) 
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The value of the second derivative at the initial instant of time need not be 
prescribed, since it can be determined from the equation. 

We assume that x € R, that is, the string is infinite. This model gives a 
good approximation to physical reality, provided we consider deviations that 
are small compared to the length of the string over brief time intervals. 

The first boundary-value problem includes the equation (4.1) for x € (0,1), 
the initial conditions (4.3), and the boundary conditions 

lil Ml Og (4.4) 
which express the fact that the finite string is clamped at both ends. 

This problem also describes the longitudinal vibrations of a rod with fixed 
ends. In this case u(x,t) is the displacement of the point x of the rod from its 
equilibrium position at time t, as in Fig. 4.3. 


Fig. 4.3. The vibrating rod P 


u(x,t) 


One can consider the vibrations of a rod with one or both ends free, ex- 
pressed by the condition 


=0 and/or uzl ,=0. (4.5) 


Uz | r=0 r= 
The second and third boundary-value problems are obtained as combinations 
of these conditions. 

The fourth boundary-value problem involves consideration of the periodic 
boundary conditions 

u(x,t) = u(x +l,t), (4.6) 

so that the function u is defined on a circle. 

We are now going to obtain the complete solutions of all these problems 
(some of them in the exercises). 


4.3. The Cauchy Problem for an Infinite String. 
d’Alembert’s Formula 
Thus, consider the Cauchy problem uy = a Uzr, Ul = (zx), tl = y(x). 


Theorem (d’Alembert’s formula). The solution of the Cauchy problem is 
given by the formula 


z+at 
Gat i 
u(t) = EEEE foe. a) 


x—at 
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OUTLINE OF THE PROOF. We know that the general solution of the equation 
has the form u(x,t) = f(x — at) + g(x + at). Consider the term f(x — at). At 
a fixed time t the graph of the function f(x — at) is the graph of the function 
f(x) shifted to the right by at when t > 0. This term is called the forward 
wave. The analogous term g(x + at) is called the reverse wave. We first put 
t = 0 in the formula for the general solution; then, differentiating with respect 
to t in this same formula, we obtain the following system: 


u| o = (£) = f(z) + g(2) , 
url, = V(x) = -a f'(x) + ag'(2) . 


Solving this system, we find f and g, substitute them into the formula for the 
general solution, and obtain (4.7). O 


Problem 2. Fill in the details of this proof. 


(4.8) 


“Movie” Problems 


Using d’Alembert’s formula, one can draw the successive “frames of a movie 
of the vibrations of the string,” knowing the graphs of the Cauchy conditions. 


Example. Let y = 0, and let the graph of y have the form shown in Fig. 4.4, 
where y # 0 on an interval of length 1. 


p 
ee a (re ,w=0 Fig. 4.4. The initial conditions 


Then by d’Alembert’s formula, which has the form u(x,t) = (p(x — at) + 
y(z+at)) /2 in the present case, we obtain the shape of the string at successive 
times, as in Fig. 4.5. 


t&l Z7 vD Fig. 4.5. The shape 
of the string at 
AEON ATN successive 


7 N Z ‘N d 
E EEE a O A o times 


Problem 3. Draw the movie for the initial conditions y = 0 and y a function 
of the form “~~ . 
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4.4. The Semi-Infinite String 


This problem includes the equation (4.1), the initial conditions (4.3), and the 
boundary conditions ul 6 = 0 (fixed endpoint) or ug ox = 0 (free endpoint). 

By combining the boundary conditions and the form of the initial data we 
obtain the following table of problems: 


boundary conditions 


free endpoint at x = 0 
ised endpolit S 1 


Problems 4—7. Draw the movies of the motion of the string for the condi- 
tions shown in the table. 


HINT. The problem for the semi-infinite string can be reduced to the problem 
for the infinite string so that the solution of the problem for the infinite string, 
when restricted to a half-line, yields the solution of the problem for the semi- 
infinite string. To do this we must extend the initial conditions to the entire 
line in such a way that the solution satisfies the boundary conditions at the 
point z = 0. 

Parity considerations are helpful in doing this. The initial condition can 
be extended to the entire line as an even function provided the relation 
Ur = = 0 holds at the initial instant of time. But will the solution be an even 
function of z at all times? First of all, this can be deduced from d’Alembert’s 
formula; second, one can make use of the following remarkable idea. 

The uniqueness of the solution of the Cauchy problem was proved in the 
derivation of d’Alembert’s formula. This equation is invariant under the trans- 
formation x +> —z. If the initial condition is even, that is, also invariant under 
this transformation, we then have two solutions u(x,t) and u(—z,t). Since the 
solution is unique, these two solutions coincide, so that u(x,t) = u(—z,t) and 
the solution is an even function. 

This is a general idea: If the problem supports some symmetry and the 
solution is unique, then the solution must also support that symmetry. 

Similarly one can use an odd extension when the condition fil 8 =) 
holds. 


4.5. The Finite String. Resonance 


The method of d’Alembert is not very convenient for solving the boundary- 
value problems in the case of a finite string. Later on we shall develop another 
very powerful method to handle this case. At present we illustrate the appli- 
cation of d’Alembert’s method to the problem of forced vibrations of a finite 
string. This problem consists of the equation (4.1), the initial conditions (4.3), 


and the following boundary conditions: u|__, = f(t), u|,_, = 0. 


32 4, The Vibrating String (d’Alembert’s Method) 


Problem 8. Find the general solution of this problem. 


HInT. The values brought from the boundary and the initial interval ¢ = 0 
along the characteristics x — at = const, x + at = const contribute to the 
solution at the point (x,t). The characteristics undergo bending as they are 
reflected from the boundary. As a result, the value of the solution is an alter- 
nating sum of the values at the nodes of the resulting broken lines, as shown 
in Fig. 4.6. The initial conditions can be arbitrary; for simplicity you may first 
take them to be zero. 


t u(z, t) 


Fig. 4.6. How to obtain the value 
of the solution at a given point 


Problem 9. Suppose the function f(t) is periodic with period T. Will there 
exist periodic solutions with periodic boundary conditions? If so, will the other 
solutions be periodic? 


ANSWER. There is no periodic solution, for example, in cases when resonance 
occurs (when the period of the external forces T is commensurable with the 
“periods of natural oscillation” 2//a). For almost all (in the Lebesgue sense) 
values of T there exists a periodic solution (at least if the function is smooth). 
This follows from number-theoretic considerations (the metric theory of Dio- 
phantine approximations). 


In conclusion we remark that all these results have required only elemen- 
tary mathematical analysis, remaining within the scope of the theory of differ- 
entiation and integration. The following sections require more powerful tools 
of analysis. 


4.6. The Fourier Method 


Let us return to the equation (4.1): utt = a? Urz. We can regard the right-hand 
side as a differential operator that maps one function space into another (or 
into itself, if we consider the space of infinitely smooth or analytic functions). 

We write the equation (and any problem we may consider along with it) 
more abstractly, for example: 


ut = Au, ul a u| _,=0 or u(x +27,t) = u(x,t). 
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We shall not specify any initial conditions, since we wish to learn how to 
solve the problem for arbitrary initial conditions. 

Thus A is a linear operator. A point of the phase space is a function u 
satisfying our boundary conditions. 

The operator A has certain supplementary properties in addition to lin- 
earity. Let us see what these properties yield, for simplicity starting with 
finite-dimensional models. 


Model 1. We consider the linear ordinary differential equation us = Au, 
where u is a vector in Euclidean space RY and A is a self-adjoint opera- 
tor: (Az, y) = (x, Ay), the parentheses denoting the inner product. Then the 
operator A has an orthogonal basis consisting of eigenvectors X,, and the cor- 
responding eigenvalues Àx are real. The functions e**' X, form a fundamental 
system of solutions. The general solution has the form u(t) = Y cge*** Xp. It is 
easy to find a solution with prescribed initial condition u(0) = ¢ (in this case 
these are a finite-dimensional vectors). We substitute t = 0: p = >> c,Xg, 
and we need to find c. We then take the inner product with X,, use the 
orthogonality of the basis, and obtain the result that c; = (Y, X;)/(X,, X;). 


The similarity to the wave equation here is that our second-order differ- 
ential operator d?/dx? is self-adjoint in a suitable space. 


Differential Operators in the Space of Trigonometric 
Polynomials 


One of Hilbert’s most remarkable ideas was to regard function spaces as 
Euclidean spaces. This idea is the foundation of functional analysis. 

But at present we shall restrict ourselves to finite-dimensional spaces of 
trigonometric polynomials. For any nonnegative integer N we consider the 
space 


k=N 
By =} ula) = >D cket, rest, aec}. 
k=-N 


As long as the summation extends only from —N to N, we remain within 
the limits of linear algebra. If it extends from —co to +00, we go into the area 
of Fourier series, which is a branch of functional analysis. 

On the space Eyn we introduce a Hermitian form: 


Problem 10. Verify that all the axioms of a space with a Hermitian inner 
product are satisfied. 


The norm of a polynomial is defined by the formula ||f||? = (f, f). 
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Problem 11. Prove that the polynomials e, = ekt ke {-N,...,N}, form 
an orthogonal basis. 


27 
SOLUTION. f ete"? dr = 2r if k = l. But if k—l = m Æ 0, then the 
0 


1 aimer 


2r 
À l 2 
integral equals J em? dr = $ e i 


- =Q. 

l 
Model 2. We consider a periodic boundary-value problem for the wave equa- 
tion with initial conditions in the space Eyn: 


Utt = Q Urr, ule Qn, t) = u(x,t), 
ul o=? uha = p,p EEN. 


In the next lecture we shall reduce this problem to a system of ordinary 
differential equations of the form uu = Au, that is, a system analogous to 
that of Model 1. 


Lecture 5 


The Fourier Method (for the Vibrating String) 


Consider the periodic boundary-value problem for the vibrating string: 
ie = Oe, Bie! SR WZ, (5.1) 
ieee OL) 7 = y(x). (5.2) 


The problem can be regarded as a Cauchy problem on the circle S1. On a 
connected compact manifold with boundary (for example, on an interval) it 
would be necessary to prescribe some boundary conditions. 


5.1. Solution of the Problem in the Space 
of Trigonometric Polynomials 


Consider the space of complex-valued trigonometric polynomials 


N 
En = [Dae » QKE c} : 
-N 


(One can also consider real-valued polynomials; the necessary and suffi- 
cient condition for the polynomial to be real-valued is a, = @—,.) 

Suppose the initial data y and w belong to Ey. Then the solution u(.,t) 
also belongs to Ey. There is a Hermitian product in Ey: 


Coe I fle)g@)de, IF = (F, f). 
Sl 


The following properties are easy to verify: 
1. The monomials (e] form an orthogonal basis of the space. This basis 
is not normalized: ||e*®]| = vV2r. We introduce the notation X; = e**. 
‘ 2 
2. The operator “L = a’ 4, mapping Eyn into itself is self-adjoint, having 
the X;, as eigenvectors: 


DX Sah Rice 
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5.2. A Digression 


Consider the Euclidean space R” and a self-adjoint nonsingular operator L : 
R” — R”. Let {Xx} be an orthogonal basis of eigenvectors and A, = —w? the 
eigenvalues. 

Consider the problem ë = La, x(0) = y, (0) = wv. 

The vector-valued functions sin(wzt) X, and cos(w,t) X, form a funda- 
mental system of solutions of this problem. Indeed, the dimension of the so- 
lution space of the problem is 2n (as can be proved by reducing the problem 
to a system of first-order equations t = p, p = La). These 2n functions are 
solutions and, in addition, are linearly independent (verify this!). 

Hence the general solution has the form 


r(t) = S (a cos(wkt) Xk + bk sin(wkt) Xk) ‘ 


It is also easy to find a solution with prescribed initial conditions: 


Se = (9, Xk) 
eat aay 
‘ ee i =. (Y, Xx) 
z(0) = ) bk kXk => bk = wr Xr, Xe) i 


5.3. Formulas for Solving the Problem of Section 5.1 


In this case we write u instead of x, we have Ey instead of R”, and L = Gn 
For simplicity we assume a = 1. We have X;, = e!**, wp = k. 
Then the general solution of the Cauchy problem has the following form: 


u(x,t) = N (ak cos(kt) + by sin(kt)) e'** , (5.3) 
ak = =< (ve) ’ (5.4) 
bk = az (We) ,k #0. (5.5) 


Here the summation on k extends from —N to N. 

Thus, the process of solving the Cauchy problem in the class Ey requires 
nothing more than linear algebra. But it is noteworthy that the solution, even 
in the general case, is given by formulas (5.3)-(5.5) (with summation from 
—oo to +00). 


5.4. The General Case 


Let y and y belong to C™(S!). (One can also consider the case of analytic 
initial data or, in the other direction, data with finite smoothness.) 


5.6. Convergence of Fourier Series 37 


Theorem. The solution of the Cauchy problem (5.1), (5.2) is given by for- 
mulas (5.3)-(5.5) with summation from —oo to +00. 


To prove this result one must verify that the series converges, that it can 
be differentiated termwise, and that the initial conditions are satisfied. 


5.5. Fourier Series 


We need to prove that the series (5.3) converges along with all its partial 
derivatives up to order two inclusive. 
Consider the space C®(S') with the Hermitian product 


(f.9) = J fgde, IFIP = (Ff). 
Si 


This space is not complete in this norm, and so it is called a pre-Hilbert 
space. (A complete space with a Hermitian product is called a Hilbert space.) 
In this space the system of functions {X,} = {e'**} forms a “basis” in a 
certain sense. 

We remark that the series (5.3) can be rewritten as a series of exponentials 
by using Euler’s formulas. The result is a series of terms of the form e/*(¢+=) 
and e!*(t-=) Such a series is easy to differentiate. 

To prove that the series (5.3) converges together with all the series obtained 
from it by termwise differentiation, it suffices to prove that the coefficients a, 
and b; tend to 0 faster than any power of k. 

To prove that the initial conditions are satisfied, it is necessary to show 
that the Fourier series of a function converges to that function. It is in this 
sense that the system {X;} is a basis: 


If pw So apXe, an =(9,Xk)/(Xk: Xe), then S > aXe gy. 


We remark that the part of the proof in which the mere convergence of 
the Fourier series of the initial data is proved (without regard to the limiting 
function) does not differ from the proof of the convergence of the serié (5.3). 
To do this, as noted above, it suffices to establish that the |a,| decrease rapidly. 


5.6. Convergence of Fourier Series 


Lemma 1. The Fourier coefficients of a function in C®(S') decrease faster 
than any power of the exponent. 


PROOF. The following estimate is obvious: 


lax] < (max lel) 27, 
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that is, the coefficients are at least bounded. We now integrate by parts: 


1 
an = — | o(z)e"** dx 
27 


== | oae) 


1 —ikx |27 —i 
~ Drik (roe Fha - Jee an) 
1 


/ -ikx 
= dz. 
saz | owe dx 


For the same reason as at the beginning of the proof, the last integral is 
bounded uniformly with respect to k. Continuing, one can again integrate by 
parts as many times as necessary. 

As a result we find that Vm 3Cm,y : |(y, e**) | < Cryp|k|~™. The lemma 
is now proved. UU 


Remark. For a function of class C™, we find that the Fourier coefficients 
decrease like 1/|k|™. 


Problem. Suppose the function y belongs to the class CĦ (St), that is, y is 
holomorphic in the strip |Imz| < 8 and periodic: y(z + 27) = (z). Prove 
that then the Fourier coefficients decrease exponentially: 


lag] < Ce PE, 


HINT. Shift the path of integration in the integral that defines the Fourier 
coefficients by +i@ (depending on the sign of k). 


Is the converse true? That is, if this last estimate holds for the Fourier 
coefficients, does the function admit a holomorphic extension to the strip? 


ANSWER. The sum of the series is holomorphic in the interior of the strip. 


Remark. If a function is continuous, its Fourier series converges to it in the 
metric of L?. The contemporary Swedish mathematician Carleson has shown 
that pointwise convergence holds on a set of full measure. 


Lemma 2. The Fourier series of a function in C? converges to the function. 


PROOF. The Fourier series of a function y in C? converges uniformly to some 
function w € C? (by Remark after Lemma 1). We must show that y = yw. 
Suppose not. Then d(x) := p(x) — y(x) # 0. Nevertheless 6 has a zero Fourier 
series: Vk (ô, X;,) = 0. 

Let us attempt to approximate 6 with a trigonometric polynomial: 


5 — Zax = Noll? + |E axl > ysl? > 0. o) 
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We see that it is impossible to get an arbitrarily close approximation. However, 
by the Weierstrass theorem, one can always approximate a continuous func- 
tion 6 uniformly to any degree of precision using trigonometric polynomials. 
That contradicts (x). The lemma is now proved. O 


Remark. One can consider the Fourier method for the problem on an in- 
terval with various boundary conditions. But we have actually developed a 
method of solving the completely general problem uş = Au, where u is a 
function on an arbitrary manifold. One need only find the eigenfunctions and 
eigenvalues for the operator A on the manifold, construct an orthogonal basis 
of eigenfunctions, and then follow the outline above. 


For the Laplacian A = A on, for example, a manifold with a Riemannian 
metric, a rich theory thereby results. It forms part of the spectral theory of 
differential operators. 


5.7. Gibbs’ Phenomenon 


The Fourier series of a discontinuous periodic function cannot converge uni- 
formly to it, but it can converge pointwise. For example, such is the case 
for discontinuous piecewise-smooth functions: the sum of the series coincides 
with the function on the intervals of smoothness. The sequence of graphs of 
the partial sums of the Fourier series of a piecewise-smooth function con- 
verges (uniformly!), but not to the graph of the original function. Rather it 
converges to a different curve. This curve is obtained from the graph of the 
given discontinuous function by adjoining vertical line segments above the 
points of discontinuity. It is interesting that these segments are longer than 
the segments joining the parts of the graph to the right and left of the point of 
discontinuity. Here the lengths of the additional vertical tails extending above 
and below the graph of the original function always constitute the same por- 
tion (each approximately 9%) of the magnitude of the discontinuity, as in 
Fig. 5.1. 


Fig. 5.1. Gibbs’ phenomenon 


To prove this amazing fact (known as Gibbs’ phenomenon) it suffices to 
carry out the computation of the maximum and minimum of the partial sum 
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of the Fourier series of a very simple discontinuous function of period 27 (say, 
the function equal to x for |z| < 7). The general case can be reduced to this 
one by removing the discontinuities through subtraction of a suitable linear 
combination of translates of this very simple discontinuous function. 

Gibbs’ phenomenon leads to some interesting artefacts in tomography 
(where the function of two variables is, for example, the X-ray density of 
a planar section of a human body, reconstructed by a computer using the 
summation of a (double) Fourier series). 

Since the density is a discontinuous function (for example, because of the 
presence of bones), Gibbs’ phenomenon arises. It manifests itself here in the 
form of additional lines to the lines of discontinuity: the double tangents and 
the inflection tangents of real boundaries of tissues of differing density (think 
about why this is so). 


1 The coefficients of this series are the Fourier coefficients of the so-called Radon 
transform of the original density. The Radon transform is the integral of the 
original density along a line, regarded as a function of the line; here it is regarded 
as a function of one variable, namely the distance of a line of a given pencil of 
parallel lines from a single one of these lines. 


Lecture 6 


The Theory of Oscillations. The Variational 
Principle 


It is an experimental fact well-known to physicists that the laws of nature can 
be described by variational principles. The absence of rational grounds for 
this fact has generated attempts to interpret it theologically, philosophically, 
and in other ways, as, for example, in the writings of Voltaire, Maupertuis 
and others in the collection Variational Principles of Mechanics. 

The variational principle asserts that “nature acts in the shortest way.” 
For example, according to Fermat’s principle, rays of light propagate along 
the shortest path. The principle of least action was stated by Hamilton in a 
form close to the modern one: 

ô / Ldt=0. 


The integral that occurs in this formula is called the action integral. The 
principle describes the motion of some mechanical system. The Lagrangian 
L(q,q,t) is a function of the state and velocity of the motion in the cor- 
responding configuration space. In this form the principle reads as follows: 
“A motion of a mechanical system q = q(t) is possible if and only if the 
variation of the action integral along the curve q = q(t) equals zero.” 

In other words, the true motion is a critical point of the action function de- 
fined on the infinite-dimensional space of smooth mappings of an interval into 
the configuration space, as shown in Fig. 6.1. Functions on infinite-dimensional 
spaces are usually called functionals. What form does the action functional 
take for a typical mechanical system? 

Consider the motion of a point q on a Riemannian manifold M”. The 
trajectory is defined by a mapping q = q(t), t € [to,t,]. We shall assume 
that the endpoints g(to) and g(t) are fixed. The tangent vector q(t) has the 
Riemannian square of the length. The quantity T = (1/2)q? is called the 
kinetic energy. We also consider a certain function U : M” — R called the 
potential energy. We set L = T — U; this functional is called the Lagrange 
function or the Lagrangian. 
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M” 
Fig. 6.1. The true motions 
tı , are critical points of the 
q = q(t) Pq] = J L (q(t), q(t), t) dt action functional 
0 
— 
to ti 


(There was a dispute in the nineteenth century over the notation for co- 
ordinates and momenta and also over the choice of sign in the Lagrangian. 
Naturally mathematicians were unhappy with the minus sign. As usual, the 
physicists triumphed, and the generally accepted notations are now q for co- 
ordinates, p for momenta, L = T — U for the Lagrangian, F = de for force. 
In any case, physical common sense is satisfied: when a rock is lifted from the 
surface of the earth, its energy increases. See Klein [2].) 


ti 
The action function has the form Sq] = f L(q(t), q(t), t) dt. 

to 
In coordinates the kinetic energy T is a quadratic form: T = 5 Saad) 
Let us introduce the momenta p; = 0T/04q;. ij 
Geometrically momentum is a linear functional acting on tangent vectors, 


that is, it is a cotangent vector. 


Theorem 1. The extremals (critical points) of the action functional satisfy 
the Euler-Lagrange equations 


a Sas qa 
(This is a system of n second-order ordinary differential equations for the n 
unknown functions q;(t).) 


pe er 


The theorem will be proved below; at this point we shall exhibit the con- 
nection the vibrations of a string have with the equation and examine several 
examples. 

Lagrange considered a “discrete string” made up of balls joined together 
with springs, as in Fig. 6.2. 


Fig. 6.2. The Lagrange model 
of the string 


It is not difficult to find the Lagrangian for such a system and derive the 
equations of motion, and then obtain the equation for the string by passing to 
the limit, understand its variational basis, and find the Lagrangian. We first 
consider some simple examples. 
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Example 1. Suppose a particle is moving freely in a Euclidean space, say the 
plane for simplicity. Then the potential energy U is zero, and the Lagrangian is 
L =T = (1/2)(q? +42). It is easy to verify that the extremals of the action are 
straight lines. Indeed, by Theorem 1 the extremals satisfy the Euler-Lagrange 
equations, which in the present case have the form pı = 0, po = 0. But pı = qy 
and p2 = q2, that is, the generalized momenta coincide with the velocities. 
Then g; = 0 and gz = 0. The solutions of these equations are lines, as shown 
in Fig. 6.3. 


Fig. 6.3. Lines are the extremals 
of the action for a free particle 


Example 2. We now introduce a potential energy U = U(q), that is, we let 
the particle move in a force field. We have the following system of equations: 


, OU j ðU 
PIS SSR SS a 
Oq1 Og 
or 
; OU E ðU 
Vi gome s 
Oq2 042 
Let us consider the case when U is a quadratic form. For example, when 
n = 1 we have U = aq’, and the equation g = —aq is the equation of pendulum 
motion. 


In the general case of a quadratic potential U we refer the quadratic form 
to principal axes via an orthogonal transformation (in the sense of the metric 
defined by the kinetic energy). We thereby obtain a system of independent 
oscillators. Returning to the vibrating string, one can say that it is an infinite 
system of independent pendulums. 


PROOF OF THEOREM 1. For simplicity, to avoid having to write subscripts, 
we shall assume that q € R. We have the Lagrangian L(q,q,t). Consider a 
variation of the motion dgq(t), d¢(to) = dg(t1) = 0, as in Fig. 6.4. 


Fig. 6.4. The variation of the motion 
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We compute the increment of the Lagrangian: 


: ; OL ðL.. 
L(q + ôq, ġ + 5q,t) — L(q,ġå,t) = oq ôq + gg OD +o(ô) . 


We now find the principal linear part of the increment of the action: 


tı 


ti ty 
OL ðL. (x) d OL OL ) 
pde=) | —“oge — 2 se ogy a di 
6 f Laatat IG, sat 5 (64) ) at rae dt wi) 1 Bq í 
i 


to to 0 
| daL ƏL 
= pe ede C Nde 
Jal a i 


The transition marked (*) is an integration by parts, taking account of 
the fact that the integrated term equals 0, since the variation of the trajectory 
at the endpoints of the interval is zero. 

If the trajectory is critical, this last expression is identically zero. Then 


dtôðġ g 


For if not, the variation ôq could be chosen so that the integral is not identically 
zero (this is an exercise for the reader). The theorem is now proved. O 


In the case of free particle motion the Lagrangian has the form L = T. It 
then follows immediately from the Euler-Lagrange equations that the critical 
trajectories are straight lines. The integral is minimized on them. For the 
problem of the shortest path, in which L = /2T, the extremals are again 
straight lines, and the minimum is reached on them. 


Important remark. This theory is independent of the coordinate system. 
For that reason, if a certain equation is an Euler-Lagrange equation, it has the 
form of an Euler-Lagrange equation in any other coordinate system. Hence, 
it suffices to perform the change of coordinates in the Lagrangian alone. 


Exercise. Consider a uniform motion q(t) of particles along lines. Then q(t) 
is the solution of the equation g = 0. Write the corresponding equation in 
polar coordinates. 


For more details on the methods of the calculus of variations, see {1, Ch. 3], 
which contains further references. 

Our theory provides a powerful method of studying systems in a neigh- 
borhood of a critical point in relation to the deviation from the equilibrium 
position of the potential energy, for example, a minimum. We replace the po- 
tential by its quadratic part. The error in the right-hand side of the equation 
(that is, in the force) will be of second order of smallness in comparison with 
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its distance from the equilibrium position. The linear portion does not change 
when this is done, so that the change just described is simply a linearization. 
We replace the quadratic form of the kinetic energy, which depends both on 
a parameter and on the position g, by a constant form taken at the minimum 
point. In the finite-dimensional case the kinetic and potential energy have the 
form 


Te ON x 1 
T = 54, å) , U = 5(Ba,q) - 


The first form is positive-definite, but the second need not be if the initial 
critical point of the potential energy is not a minimum. 

It is known that two such forms can be simultaneously referred to prin- 
cipal axes. This is done by bringing the second form to principal axes via 
an orthogonal transformation for which the first form defines the Euclidean 
structure of the space. To do this in practice it is necessary to solve the char- 
acteristic equation det(B — AA) = 0. This equation arose historically in the 
work of Lagrange on the secular perturbations of the orbits of planets near 
Keplerian orbits. For that reason the characteristic equation was called the 
secular equation. 

Thus in the theory of small oscillations, by definition, the kinetic energy 
T defines a Euclidean structure in the configuration space and the potential 
energy U is a quadratic form in that space. In Cartesian coordinates the 
Euler-Lagrange equations have the form g = -VU = —Bgq. 

There is a simple geometric procedure for finding the principal axes. Con- 
sider a level ellipsoid of the quadratic form (Bq, q). The point farthest from 
the origin in the sense of the given Euclidean structure determines the first 
eigenvector (principal axis). The other principal axes lie in the orthogonal 
complement of the first, and this hyperplane intersects the ellipsoid in an el- 
lipsoid of smaller dimension, in which the same procedure can then be applied, 
and so on, as in Fig. 6.5. It is surprising that this method also works in an 
infinite-dimensional space. 


> 


Fig. 6.5. Finding the principal axes 5 


Nr / 
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In principal axes the system splits into one-dimensional equations that 
can be solved explicitly. (We write out the solutions in the case of a positive- 
definite form B): 


Ük = Akk, Q(t) = ax cosuwzt + bksinwgt, where wk = Nee: 

These formulas describe the independent natural harmonic oscillations 
along orthogonal directions. If the frequencies are incommensurable, the mo- 
tion as a whole is non-periodic. 

It is convenient to write the solution corresponding to the natural oscil- 
lations in the complex form q(t) = Re(A,el**té,). Here Ax is the complex 
amplitude and £x is the eigenvector. All possible solutions are obtained by 
adding these solutions: every oscillation is a superposition of natural oscilla- 
tions. The original system thus turns out to be a system of n noninteracting 
oscillators. 

As the potential energy U increases all the natural oscillations increase 
(for the proof see |1]). From this we obtain the following simple but amazing 
facts about the ellipsoids: 


e if one ellipsoid lies inside another, each axis of the smaller ellipsoid is less 
than the corresponding axis of the larger; 

e consider a section of the ellipsoid by a plane passing through its center; 
the axis of the section has a length intermediate between those of the axes 
of the original ellipsoid, as in Fig. 6.6. 


a2 < ai 


T bə < by 
\/ a<b<e 


A<B 
axA<b<B<e 


Fig. 6.6. Relations between the axes of ellipsoids 


This effect is observed in practice, for example, when a crack in a bell is 
filled in; the cracked bell has a dull sound, which becomes more resonant when 
the crack is filled in. 

For more details on small oscillations see [1, Ch. 5] and also [3, § 25, par. 6]. 
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To generalize our theory to infinite-dimensional systems we need some 
groundwork. But first let us see what results are to be obtained. 

In the case of a one-dimensional string with endpoints fixed the configu- 
ration space is the space of functions u(x): u(0) = u(l) = 0. We shall assume 
that these functions are smooth. The kinetic energy has the form 


Lp (aur? 
u 
0 


The potential energy is the work done in bringing the string into the given 
shape u(x). We shall assume that this form differs only by a small amount 
from the stationary position u = 0. In computing the potential energy in the 
approximation provided by the theory of small oscillations we need to take 
account of terms of order u?, but we can neglect terms of higher order (u? and 
so forth). Since no work is expended in bending the string, all the work goes 
into lengthening the string from its original length l to a new length, that of 
the graph of u(z). 


Lemma. The potential energy of a tightly stretched string is proportional to 
its lengthening (in the approximation of small oscillations) with a coefficient 
proportional to the tension. 


PROOF. By Hooke’s law the tension is proportional to the lengthening. 
Lengthening is a quantity of second order of smallness relative to the length of 
the string. Consequently, the tension in the deformed string can be regarded 
as constant and independent of the shape of the string. (This approximation 
is valid within the limits of the approximation of small oscillations, that is, 
when the energy is computed taking account of quantities of second-order 
smallness relative to the deviation but not third-order. In this approximation 
it suffices to carry out computations retaining terms of first order of smallness 
relative to the deviation and neglecting terms of second-order smallness.) 

If the tension is constant throughout the deformation, the energy of an 
element of the string can be computed as the work of this constant force on 
an element of the path and therefore (in this approximation) proportional to 
the lengthening of this element (and the magnitude of the tension). 

Summing the potential energies of all the elements of the deformed string, 
we find that in the present approximation its potential energy equals the 
product of the tension and the lengthening of the entire string. 

Neglecting infinitesimals of second order relative to u, we finally obtain 
the following expression for the potential energy: 


b b 
Ou 2 F Ou 2 
u=F | i+ (Fe) =] aag f (FE) dz. g 
0 0 
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Remark. The same expression could have been derived by considering the 
model with beads and springs and computing the force acting on the it” bead 
from sides (2 — 1) and (7 + 1) (assuming that this bead moves in the direc- 
tion perpendicular to the unperturbed string). The projections of these two 
forces in the direction normal to the unperturbed string are in first approxi- 
mation proportional to the differences in the displacements of successive beads 
Q—Q—-1 and g,41—q, respectively (with opposite signs and common coefficient 
F equal to the tension in the string). 

The potential energy U = ` £(h41 —q,)* leads precisely to the following 
expression for the force acting on the ¿t bead: 


= Fl (Q41 =.) (a -— q—1)] = F(q41 — 2q + G-1) - 


In the limit the sum becomes the integral obtained in the proof of the lemma. 


The Lagrangian L is T — U. 

To write the Euler-Lagrange equation, let us find the momenta and forces. 

Obviously p(x) = Ou/Ot (assuming the string is homogeneous and has 
density 1 throughout). 

The potential energy in the pre-limit case equals E (q1 — q )?, where 
qı = u(z,) as in Fig. 6.7. 


Q 
Fig. 6.7. The discrete model 
Tı of the string 


To find the force we take the derivative with respect to q,, obtaining 


=F (q41 — Q) + F(q — G-1) = -F (Q41 — 24: + Q-1). 


The expression in parentheses on the right is precisely the second difference 
(Fig. 6.8), which becomes the second derivative 0?u/0x? upon passage to the 
limit. 


Fig. 6.8. The second difference 


Thus, after a second passage to the limit the Euler-Lagrange equation will 


have the form 
8u 07u 


ðt? Ox?.” 
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that is, the form of the wave equation. The coefficient k is directly proportional 
to its tension and inversely proportional to its density. We already have a 
method of solving it — finding its natural oscillations. It is necessary to find 
the principal axes of the quadratic form of the potential energy in the metric 
given by the kinetic energy. In the finite-dimensional case we wrote: 


q(t) = Ate . 


Thus for the string it is necessary to seek a solution of the form u = e”#€(z). 
We substitute this expression into the equation: 


zw eer] = kete" (x) ; 


thus arriving at the eigenfunction problem £” (x£) = —A€, À = w? /k. 
This application of the general theory of oscillations is called the method 
of separation of variables. The solution is: 


E= acos VAr + bsin Vz; €(0)=0>a=0; 


ED =0, sin YVA =0, vVM=nr: vi= 2, WNN. 

We see that there is a countably infinite set of principal axes and eigen- 
functions; as the number n increases the natural frequencies increase also. 
When an eigenfunction is multiplied by the real part of a complex exponen- 
tial function of time we obtain a standing wave (Fig. 6.9): the shape of the 
wave does not change, but oscillations harmonic in t occur. 


Fig. 6.9. Standing waves 


In the finite-dimensional case any oscillation could be decomposed into 
superpositions of natural oscillations. Such is the case also in the equation of 
the vibrating string, but it requires a separate proof. (This proof is provided by 
the theory of Fourier series, which was created for this purpose. For that reason 
the theory of Fourier series and its generalizations are also called harmonic 
analysis.) 

A justification for the possibility of applying the theory of oscillations to 
systems with an infinite number of degrees of freedom (like the string) was 
obtained by mathematicians only rather recently (in the late nineteenth and 
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early twentieth centuries), and it led to the creation of functional analysis and 
then to quantum mechanics. However, the most general algebraic structure of 
the theory of oscillations with both a finite and an infinite number of degrees of 
freedom is a much more fundamental theory (used by physicists with complete 
success long before its rigorous mathematical justification). 

The situation here resembles the situation in regard to justifying the theory 
of real numbers (obtained, strictly speaking, only in the nineteenth century). 
The ancient Greeks of Pythagoras’ time were compelled to conceal the the- 
orem on the incommensurability of the diagonal and side of a square, which 
had destroyed their belief in the power of numbers (which at the time were 
understood only as rational numbers). 

This did not prevent Newton from creating analysis without bothering 
about the details of justifying the arithmetic of real numbers (which, however, 
he knew very well). The author of the present course has attempted to teach 
the audience how to conjecture and predict bold generalizations (such as in 
the passage from oscillations with a finite number of degrees of freedom to 
the oscillations of continuous media) rather than the difficult skill of giving a 
rigorous justification of these results (including the unavoidable investigation 
of the minimum necessary smoothness of the objects in question). 

At present quantum field theory finds itself in a position analogous to 
that of the theory of oscillations of continuous media and the theory of real 
numbers. It provides very impressive mathematical results, but not their jus- 
tification. 


Problem*! (Sturm, Hurwitz). Consider a real 27-periodic function 


f(x) = ` (ax cos kx + by sin kz) , 


k>N 


whose Fourier series begins with the Nt” harmonics. Prove that the number of 
zeros of f on the circle {x mod 27} is not less than the number of zeros of the 
first harmonic with nonzero coefficient in the series. (That is, it is at least 2N.) 


Example. When N = 1, the number of zeros is at least two. This is Morse’s 
inequality: a function on the circle has at least two critical points. 
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From now on, the superscript * in similar situations denotes more tough problems 
whose solutions are not supplied therein. 


Lecture T 


The Theory of Oscillations. The Variational 
Principle (Continued) 


Thus the kinetic and potential energies of the string (Fig. 7.1) have the form 
bs k! 

T = = fù? dz, U == f (uz) dz. 
2 53 


Fig. 7.1. The vibrating string 


The Lagrangian L equals T — U. The variational principle ô f Ldt = 0 


OL OL 
implies the Euler-Lagrange equation Ta T which in our case has the 
q 
form of the vibrating string equation: 
Utt = küsz : 


These results can be generalized to the multidimensional case. Let the vi- 
brating body be of arbitrary dimension. For example, when n = 2 it is a 
membrane. The membrane, like the string, can be pictured as the limit of a 
system of vibrating balls connected with springs, as in Fig. 7.2. 

For a membrane or a multidimensional body the kinetic energy can be 
expressed in analogy with that of the string: T = 4 [f(du/dt)? dz. 

By the integral here we mean the corresponding multidimensional integral. 

We remark that one can consider vibrations not only in one direction, 
but also in more than one. In that case u becomes a vector-valued function. 
(Geometrically the equilibrium position corresponds to the zeroth section of 
the bundle, but the bundle itself can be arbitrary.) The kinetic energy will be 
written by the same formula, except that the integral will be the square of 
the norm of the velocity vector. 
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Fig. 7.2. Approximations 
of the string and 
membrane ANANAS ONIN OT OO" membrane 


How is the potential energy expressed? Here is the answer: 


U = 5 | (Yu)? dr, 
where Vu = ĝu = (2+, TES 2) is the gradient. This integral is called the 
Dirichlet integral. 

Dirichlet studied mostly number theory, but he is the one who noticed 
that harmonic functions minimize this integral (under suitable boundary con- 
ditions). This physical idea - the so-called Dirichlet principle — has turned out 
to be a very powerful tool both for proving the existence of solutions of the 
corresponding problems, but also for studying them and even for practical ap- 
proximate computation. It is interesting that attempts to prove the Riemann 
conjecture on the zeros of the zeta function are connected with this integral 
and its generalizations. The conjecture is that all the non-trivial complex zeros 
of the zeta function lie on a single real line. The idea (which goes back at least 
to Hilbert) is to find an oscillation problem through whose eigenfunctions the 
zeros of the zeta function can be expressed. The fact that the eigenvalues are 
real would then imply that all the zeros lie on a single real line in the complex 
plane. 


Problem. Find the minimum of the Dirichlet integral in the space of smooth 
functions on the sphere. 


SOLUTION. For a constant the integral equals zero, while for a nonconstant 
function it is positive. Consequently, the minimum equals zero and is attained 
only on constant functions. 


In the multidimensional case the Euler-Lagrange equations have the form 
2 2 2 
ve = kAu, where A = £5 freed a is the Laplacian in Cartesian coor- 


dinates. Actually the Laplacian depends not on the coordinates but on the 
Euclidean structure of the space, as the next result shows. 


Theorem 1. Au = div(grad u). 
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Before proving the theorem we recall the definitions of the concepts that 
occur here and verify that they are invariant. 

Let u be a function. Its differential du is a 1-form acting on a tangent 
vector €: du| (£) = (az,&), since any 1-form in Euclidean space is the inner 
product with a fixed vector at the given point. This vector is called the gradient 
of the function u at the point z: a, = grad ti Thus the gradient is determined 
by the Euclidean structure of the space (or more generally, by the Riemannian 
metric of a manifold). 

In any coordinates we have du = Ces dz, +--+ oe dës BY &. 
In orthonormal Cartesian coordinates this is the notation for the inner prod- 


uct of the vectors € and (g4,..., 2%), so that in orthonormal Cartesian 
coordinates the gradient of u has the components (2, ee ou), 


Exercise. Find the components of gradu in polar coordinates on the Eu- 
clidean plane and in spherical coordinates in three-dimensional Euclidean 
space. 


The divergence of a vector field is defined on a manifold on which a volume 
element is defined, in particular, on a Riemannian manifold. 

Consider a vector field v = D w(x) a2 and the corresponding equation 
£ = v(x) of the phase flow in n-dimensional space. 

Let us recall the definition of the flux of a field across a surface. We denote 
the volume element by r = 7”. The volume element 7” is an n-form. Consider 
the (n — 1)-form 7,7 corresponding to the vector v. This form is obtained by 
substituting the vector v in the first argument of 7” and leaving the other 
n — l arguments free. 

The flux of the field v across an (n — 1)-dimensional oriented hypersurface 
is the integral of this form over the surface. 


Example. Suppose T = dz1 ^- --Adz£n in orthonormal Cartesian coordinates. 
Then 


lyT = V1 d2 ^+- - A dEn — v dx, Adr3A--- Adan +- -EUn dr, A---Adry_] . 
In particular, when n = 3 for the velocity field of the system t = P, y = Q, 
z= R, we obtain 1,7 = P dy Adz + Qdz Adz + Rdz ^ dy. 


The flux has a clear meaning in fluid dynamics. It represents the amount 
of the phase fluid passing through the surface in unit time, as in Fig. 7.3. 


U7 


Fig. 7.3. The flux of a vector field 
through a surface 
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Consider the flux V, of the field v through a small sphere of radius € with 
center at x. The divergence of the field v at the point z is the limit lim V./T(E), 


where 7(€) is the volume of the ball. The meaning of the divergence is the 
“density of the source” at the given point. 

Thus the divergence is defined if a volume element is given. In particular 
the divergence is defined in Riemannian and Euclidean spaces. The divergence 
of a field v on a manifold with volume element 7 is connected with the exterior 
differentiation of (n — 1)-forms by the relation 


(divv)r = d(iyr) , 
which is completely independent of the coordinate system. 


Theorem 2. In orthonormal Cartesian coordinates the divergence of the 
field v is the trace of the Jacobian matriz of v: 

Ov} OUn 

divv = —+---+—. 

Ox, OLn 
PROOF. Consider the new position of the point x under motion along the 
phase flow of the field v over a small time interval of length £: x + x+ev(x)+ 
ole). Thus the volume element undergoes a small linear transformation, as in 


Fig. 7.4. 


Fig. 7.4. A small linear transformation of the volume element 


Lemma. In first approximation only the changes in vectors in their own di- 
rection contribute to the change in the volume of the parallelepiped they span. 


Exercise. Prove the lemma. 


Since the parallelepiped spanned by n vectors can be taken instead of a 
sphere in the definition of the divergence, when computing the flux one must 
find the increment of the volume of this parallelepiped under a small trans- 
lation along v. According to the lemma the principal term of this increment 
is the trace of the Jacobian multiplied by the volume of the original paral- 
lelepiped. Dividing by that volume and passing to the limit, we find that the 
divergence equals the trace. Theorem 2 is now proved. O 


In Cartesian coordinates we find by an elementary computation that 
Au = div(grad u). The Laplacian is defined by this formula in any Riemannian 
manifold. 
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Problem. Compute the Laplacian of a function defined on the circle in the 
Euclidean plane. 
2 


ANSWER. Au = oe where zx = cosy and y = sin. 
Ou u., ee 
The formula Au = z5 +-+ >a İS valid only in orthonormal Cartesian 
oxi Ore 


coordinates. The operator has different expressions in other coordinates, even 
when the space is Euclidean. 


Exercise. Write the Laplacian in polar, cylindrical, and spherical coordi- 
nates. 


Remark. As will be shown below, the geometric meaning of the Dirichlet 
integral is the principal term of the increment in area of the membrane. It is 
remarkable that it is a quadratic expression. 

This is completely obvious in the one-dimensional case: when the angle 
is small, the difference between the hypotenuse and a leg is of second-order 
smallness (verify this), as shown in Fig. 7.5. This simple fact has great conse- 
quences and important applications. 


b b=a0e+O(e*), c=atO(e?) 
a 


Fig. 7.5. The hypotenuse is practically equal to the leg 


Examples. 1. If you return home along a sinusoidal curve you increase the 
distance you must walk only slightly (by approximately 20%), since only the 
portions of the curve with steep slope make a significant contribution to the 
increase in distance, and these are only a small portion of the curve, as shown 
in Fig. 7.6. 


Fig. 7.6. A sinusoidal 
curve is not much 
longer than a 

straight line 


2. By inclining the motor slightly from the longitudinal axis of an airplane 
one can prevent the tail assembly from burning up in the wake, with a loss 
of thrust of second-order smallness, as shown in Fig. 7.7. For example, even 
at a large inclination of 6°, we obtain € ~ 0.1, which yields a loss in thrust of 
only 0.5%. 
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Fig. 7.7. Rotating 
the motor makes 
practically no change 
in the thrust 


3. By analyzing the results of Tycho Brahe’s astronomical observations! 
Kepler at first believed that Mars moved around the Sun in a circle, but that 
the Sun was not at the center. Indeed, if we consider an ellipse with semiaxes 
a and b and small eccentricity e, as in Fig. 7.8, we have b = avl — e? = 
a(1—e?/2+---). When e is small (about 1/10 for Mars), it is difficult to observe 
the difference between the ellipse and a circle, even though the difference 
between the center and the focus is noticeable. 


Fig. 7.8. An ellipse with small eccentricity 
is almost indistinguishable from a circle 


This can also be confirmed by the following simple experiment. Add one 
drop to a circular cup of tea near the center. After reflecting from the walls 
of the teacup, the waves will converge to a point symmetric with respect to 
the center, as in Fig. 7.9. 


Fig. 7.9. Convergence 
of waves in a disk 
and in an ellipse 


' These observations were made with his unaided eye. Even at the end of the 
seventeenth century, it was still remained to be proved that observations made 
through a telescope gave precision at least as good. 
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EXPLANATION. Ifa drop is added precisely at the focus in an elliptical teacup, 
the waves will converge at the other focus after reflection. A circle differs by 
only a small amount from an ellipse with small eccentricity. For that reason the 
waves converge at the second focus of this approximate ellipse (not precisely, 
but slightly spread out). With careful observation one can detect a second 
convergence of the waves at the original point. 

We can now explain why the difference is of second order of smallness. 
A line minimizes the length functional, so that the increment cannot be of 
first-order smallness. 


Let us now return to the geometric meaning of the Dirichlet integral. It 
is asserted that this integral is the principal part of the increment of area, 
that is, if So is the area of the membrane in equilibrium and Se its area at 
the displacement £u, we need to prove that the new area S equals So + 
5 f (uz)? dz + o(e”), as illustrated in Fig. 7.10. Let us look at the change in 
area of a small region. We choose a convenient coordinate system: one axis 
goes along gradu and the other is orthogonal to it, as shown in Fig. 7.11. 
Then the new area is the product of the new lengths, one of which has not 
changed, while the other is computed as in the one-dimensional case. As a 

result we have Se = So(1 + (€7/2)(Vu)? +---). 


Se 


B a a 


Fig. 7.10. Change in the area of a membrane 


Se 


Fig. 7.11. Computation 
of the increment 
in the area 


Ch EU 


It remains to integrate this increment over the entire membrane. 

We can carry out all computations in the pre-limiting system of balls and 
springs shown in Fig. 7.12. 

We construct a homotopy au from the equilibrium position to the dis- 
placement £u; in doing so we need to find the term of order £2 in the ex- 
pression for the potential energy, that is, the work done in making such a 
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(TY +h) 
p (J 
u(x — h, y) u(x + h, y) Pg Vs l 
"ula, y— h) 


Fig. 7.12. Computation of the potential energy of a membrane 


displacement. The side view in both orthogonal directions is the same as in 
the one-dimensional case (Fig. 7.12). Hence the balancing forces are already 
computed and consist of the second differences with a negative sign (with a 
suitable choice of units of measurement): 


—u(x — h, y) + 2u(x,y)— ulz +h, y), —uļz,y-— h)+2u(z,y)— ulz, y+ h). 


By passing to the limit we obtain terms of the form —uzz and 
Integrating with respect to a from 0 to €, we obtain 


—Uyy- 


2 
- f adudaw = “udu | ada = -2 uAu. 


Finally, by integrating over the region we obtain (up to a numerical multiple) 


u u 
v=- ff (Sat a) eddy. 


We now integrate by parts (first with respect to x) 


OU a geen -+ a 
T Dion T 


T2 


Tı 


We now take into account the fact that u = 0 on the boundary (the membrane 
is clamped). Thus the integrated term vanishes. The integration with respect 
to y is similar. As a result, we find 


U = Jf vwracay. 


Thus the potential energy is proportional to the increment of area. As we have 
seen, the Euler-Lagrange equation for this case is the wave equation. 
In the course of the proof we have proved the formula 


[ (vw? ae = - f Auude Gig) 
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In addition, we have computed the force, which turned out to be —Au, so 
that the wave equation can be interpreted as a Newtonian equation. But we 
can legitimately compute oe once we know T and U. The latter are quadratic 
forms on an infinite-dimensional function space (q¢ = uz, q = u). 
Let us find the variation 6 f 1 (du)? dy. Let u = uo(x) + €€(x), Elan = 0. 
We need to find the coefficient of £ in the expansion of this integral and 
represent it as an infinite-dimensional inner product: 


1 (dug EN”, fl (du \* Jug O€ 
1 Ouo 
-e feauar+ [5 (32) dz +- 


(in the integration by parts we used the fact that € | ap = 9). The first of these 
integrals happens to be the inner product. Hence the force is indeed —Au. 

All this reasoning is also valid for a Riemannian manifold. 

It is a remarkable fact that a large class of problems is described by the 
wave equation. The wave equation can be derived starting from the assump- 
tion that the given variational problem satisfies certain conditions of an ax- 
iomatic nature. (To be sure, the variational derivation of the problems of 
mathematical physics remains without any convincing explanation.) The ax- 
ioms are the following: 


1. The kinetic energy T is a quadratic form in the velocities, which entails 
the spatial homogeneity of the system. However, when the density is variable, 
T is still expressed by a similar form. 


We confine ourselves to the case of a homogeneous medium, in which T 


is the integral of ee We shall exhibit the axioms that U must satisfy in 


order for the problem to be described by a wave equation. 


2. The quadratic form U is additive and local, which physically means 
there is no distant action. The quadratic form can be written as U = $(Au, u), 
where A is a linear operator and the parentheses denote the inner product 
(that is, the integral of the product of the functions). An operator A is local 
if the value of Au at a point can be computed in terms of a finite number of 
derivatives of the function at that point. 


In our case only neighboring points interact, which leads to the appearance 
of only first-order derivatives uy in the Dirichlet integral (and only second 
derivatives in the corresponding symmetric operator A = —A). The inter- 
action of a large number of points leads to the appearance of higher-order 
derivatives. For example, the operator A?u appears in the equation for the 
problem of the bending of a thin plate. In the case of a membrane the increase 
in energy occurs only because of the stretching, not bending. For a plate en- 
ergy is expended in bending as well. This is also an example of a local system, 
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but with higher-order derivatives. If only first-order derivatives occur in the 
expression for the potential energy, it will have the general form 


Ou Ou Ou 9 
) i Sa, Da; + 2u Ore : 


In which cases does the corresponding symmetric operator become the 
Laplacian? When further restrictions are imposed on the operator acting on 
the functions in the Euclidean space R”: 


3. Homogeneity (translation-invariance): 
A(u(x + a)) = (Au) (x +a). 
4. Isotropy (rotation-invariance): 
A(u(gz)) = (Au)(gr) ; 


The quadratic form in the gradient vector defines a level ellipsoid in the 
tangent space at each point. The space is isotropic if this ellipsoid is a sphere, 
that is, the properties of the membrane are independent of rotation. Non- 
isotropic media are also encountered. 

It follows from the homogeneity and isotropy that b = 0, since the linear 
part of the operator is the inner product of 6 and u2, If b were nonzero, the 
properties of the membrane would change under some rotation. 

The term cu? may occur in a homogeneous isotropic medium. For example, 
it occurs in the problem of the oscillations of a balloon (thin shell) about an 
equilibrium position. The absence of this term is caused by the following 
restriction. 


5. The metric of the space (x, u) is invariant under translations along the 
u-axis. 


All these restrictions lead to a situation in which only the Laplacian re- 
mains in the equation. 

Now let us consider the problem of the equilibrium positions — the critical 
points of the potential energy. If ile o = 9, then 0 is a stable equilibrium 
position, a minimum of the Dirichlet integral. But other boundary conditions 
can be considered. For example, consider a string with endpoints clamped at 
arbitrary points, as in Fig. 7.13. 


Fig. 7.13. The Dirichlet problem 
for the vibrating string 


In this case the minimum of the potential energy is attained on a line. 
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We now pose the analogous problem for the vibrating membrane: Au = 0, 
u| 90 = ¥: This is the Dirichlet problem for Laplace’s equation, shown in 
Fig. 7.14. The solutions of the equation Au = 0 are called harmonic functions 
in the domain §2. We shall seek a solution of the Dirichlet problem in the class 
of functions that are continuous on the closure of 92 and twice differentiable 
in its interior. 


Fig. 7.14. The Dirichlet problem 
for the vibrating membrane 


an 


There is a general device for finding a minimum of the potential U. To 
do this one must move along the gradient, that is, along the vector field 
q = —VU. When the quadratic form U is positive-definite, we necessarily 
reach a minimum in this way. In our case the equation for gradient descent 
has the form (¢ =) a = k Au (= — grad U). In particular, for the string we 
obtain the equation ĝu = kuzz, which is called the diffusion equation or the 
heat equation. 

The heat equation is obtained in the problem of heat propagation. Consider 
elements being heated that are located at the nodes of a lattice. At each 
successive moment of time the value of the temperature at a node is some 
weighted mean of its values at five points, as in Fig. 7.15. For simplicity we 
draw the one-dimensional case (also Fig. 7.15). 

Taking the average, we obtain a new function, to which the same proce- 
dure is applied, and so forth. This is precisely an implementation of gradient 
descent. Thus the propagation of heat occurs just as a minimum of potential 
energy can be sought using the gradient method. 


Fig. 7.15. Equalization 
of the temperatures 
at the nodes 


If the equation is interpreted as a diffusion equation, u corresponds to the 
density of the distribution of particles. The Laplacian arises as the result of the 
invariance conditions for the system described above. Incidentally, in this case 
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the term cu in the equation (corresponding to cu?/2 in the potential energy) 
is often physically justified: it corresponds to the creation or annihilation of 
particles, depending on the sign of c. 

A method of finding stationary solutions (equilibrium positions) has thus 
been described. 

In special cases it is possible to find explicit solutions. We first study 
harmonic functions. Being solutions of the linear equation Au = 0, they form 
a vector space. Let us attempt to satisfy the boundary conditions. When 
n = 1 harmonic functions are linear functions and the space of them is two- 
dimensional. A basis is formed by the functions 1 and xz. One can pass a 
(unique) line through two given points and satisfy the boundary conditions. 

Now let n = 2. 


Theorem 3. Let us identify the real plane with the complex line. Let f(z) be 
a holomorphic function. Then Re f and Im f are harmonic functions. 


Exercise. Prove the theorem. Use the Cauchy-Riemann equations. 


Problems. 1. Let y(xz,y) be the oriented angle subtended by a fixed line 
segment in the plane when seen from the point (x,y). Prove that y(z,y) is a 
harmonic (multi-valued) function, except at the endpoints of the line segment. 


2. Construct a bounded function that is harmonic inside a disk, continuous 
everywhere except at two given points on the boundary, equal to 0 on one of 
the two arcs into which the two points divide the boundary, and equal to 1 
on the other arc. 


3. Carry out the construction of the preceding problem for the case when 
the circle is divided into n arcs and n boundary values are given. 


4. Find all the harmonic functions on the sphere, that is, solve the equation 
div grad u = 0, where u is a smooth function on the sphere in R”. 


HINT. If a function is not a constant, its Dirichlet integral is nonzero, and the 
function cannot be a critical point of the Dirichlet integral. Therefore Au Æ 0. 


5. A function is given at time 0 on the vertices of a cube. At each successive 
time (i+ 1) the value at each point is replaced by the arithmetic mean of the 
values at the adjacent vertices at the preceding time (2). Find the limit of the 
resulting sequence of functions as i — +00. 


HINT. The limiting function is “harmonic.” 


6. Let u be a homogeneous function of degree 0 on R”, that is, u(Ar) = 
u(x) for any x and any A > 0. Prove that 


r? Au = Au, 


where A is the Laplacian div grad on the unit sphere, that is, (Au)pa1 = 
div grad (u| _,). 
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7. Let u be a homogeneous function of degree k in R”, that is u(Ar) = 
AFu(x) for any x and à > 0. Prove that 


Au = r° Au — (k? + (n — 2)k)u, 


where Au is a homogeneous function of degree k coinciding with the function 
div grad(u| _,) when r = 1. 
In particular, for homogeneous functions of degree k on the plane (n = 2) 


Au =r? Au—k?u, 
and for homogeneous functions of degree k in three-dimensional space 
Au = r? Au — (k? + k)u . 


8. Find all homogeneous harmonic functions in R” \ {0} depending only 
onr. 


HINT. The case n = 2 is special. 


9*. Prove that the number of domains into which the zero level curve of 
the kt! eigenfunction of the Laplacian in a complex domain (with Dirichlet 
boundary conditions) divides this domain is at most k. This result holds on 
compact Riemannian manifolds of any dimension, for example, on spheres. 

For example, the first eigenfunction (having eigenvalue of smallest absolute 
value) has no changes of sign at all (it vanishes only on the boundary). 


HinT. The ratio of the Dirichlet integral to the integral of the square of an 
eigenfunction is the same on each of the N components of the complement of 
the zero level curve of the eigenfunction. Hence there exists an N-dimensional 
space of functions on which this ratio of quadratic forms is the same as for 
the eigenfunction. 

It follows from the theorems of Lecture 6 on the axes of an ellipsoid that 
the index of axis k is not less than the dimension N. 


10. Consider the eigenfunctions of the Laplacian on the n-dimensional 
torus: 
Oru 
ox? 
We denote by N(E) the number of eigenfunctions for which \ < E. Study the 
behavior of N(E) as E — oo. 


=-u, (z, mod a) ET”. 


HINT. The Laplacian commutes with translations. Therefore the eigenfunc- 
tions are exponentials 


erla) = e E (k, z) = X kj 2,/a, » kp ee. 


Consequently the question reduces to computing the number of points with 
integer coordinates in a large ellipsoid. 
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ANSWER. In the cotangent bundle to the torus T*T™” consider the domain 
§2(E) defined by the inequality }° £? < E (£; are the coordinates of the form 
>); da;). Then as E — œo 


N(E) ~ (2)~" Vol R(E) = const. E"/? 


Remark 1. The analogous Weyl formula holds on any manifold and for any 
“elliptic” operator, for example, for an operator of second or higher order 
with variable coefficients P(z,iZ) = Xu. Here the domain {2 is defined by 


the condition P(x, £) < E, where P is the sum of the terms of the polynomial 
P that are of highest order with respect to € (the “principal symbol”). 


Remark 2. The later terms of the asymptotic expansion of N (E) are difficult 
to study, even for the usual Laplacian when a; = 27. The problem is that 
the surface of the sphere $` k? = E may contain many points with integer 
coordinates (in other words, the eigenvalue E may be of large multiplicity). 
This possibility occurs, for example, in the case of the Laplacian on the sphere 
S?, for which the eigenvalue E has multiplicity of order VE (see Lecture 11). 


Lecture 8 


Properties of Harmonic Functions 


In this lecture we shall study the properties of harmonic functions, mostly in 
the plane, although some theorems will be proved in the n-dimensional case. 
We shall begin with the three-dimensional case. 

Consider the function u = 1/r, where r = \/x? + y2 + z2 and (x,y,z) are 
Cartesian coordinates. 


Exercise. Verify that the function u is harmonic. 


This simple example is actually very important. The reason is that the two 
most important forces studied in physics - gravitational and electrostatic — 
are structured as follows: they act along the line joining the particles and their 
magnitude is inversely proportional to the square of the distance between the 
particles. 

But in that case the potential whose gradient is the force of attraction 
to one particle is proportional to 1/r, which is a harmonic function. A linear 
combination of harmonic functions is harmonic, so that the potential of the 
force of attraction toward a finite number of particles is also harmonic. 

From a finite number of particles one can pass to a continuous mass (or 
charge) distribution with density p(x) in a region D. Then the potential of 
the field of attraction toward the body in question at a point x is computed 


by the formula 
{ p(y) dy 
Jz7—yl 
D 


as shown in Fig. 8.1, the function being a harmonic as a function of x in a 
region free of mass (resp. charge). Thus the fundamental laws of nature are 
connected with harmonic potentials. 

Let us attempt to understand why the function u = 1/r is harmonic, that 
is, why div grad u = 0, without doing any computation. 
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Fig. 8.1. A gravitational field 
D 


First of all, at the point xz, div gradu depends only on r, the distance 
from x to 0. Thus div grad u(x) = f(||x||). Now consider the region between 
concentric spheres of radius r and R, with r < R. Let us compute the flux of 
grad u across the boundary of this region. The gradient of u is directed toward 
O and hence is perpendicular to these spheres. The area of the outer sphere 
is (R/r)? times as large as that of the inner sphere. But the magnitude of 
the gradient on the outer sphere is less than on the inner sphere by exactly 
this same factor (since the force is inversely proportional to the square of 
the distance). Therefore the fluxes of the gradient across the two spheres are 
exactly equal, so that the total flux of the gradient across the boundary of the 
spherical shell is exactly 0. It follows that div grad(1/r) vanishes. Of course, 
the equality A(1/r) = 0 can also be verified by simply differentiating. 

We now turn to the case of the plane. Consider a holomorphic function f(z) 
on the plane or in a domain of the plane. Then from f(z) one can easily form 
the harmonic functions Re f(z), Im f(z), In|f(z)|, and arg f(z) (the last two 
being the real and imaginary parts of the holomorphic function In f(z)), and 
the like. The converse is also true. Every harmonic function can be obtained 
in this way. 


Theorem 1. Every harmonic function in a simply connected domain of the 
plane is the real part of a holomorphic function, which is defined in that do- 
main up to an additive purely imaginary constant. 


PROOF. We are seeking a holomorphic function f = u + iv, the harmonic 
function u being given and v to be found. By the Cauchy—Riemann equations 
Uz = Vy and uy = —vz, so that du = a, where a = (—u,)dax + (uz) dy. 
In a simply connected domain there exists a function v with this differential 


x 
provided da = 0. In the present case we can define v to be f a, and the integral 


To 
will be independent of the path within the domain, since by Green’s theorem 
the integral over any closed path is zero. Thus it remains only to verify that 
da = 0. But da = (uzr + Uyy) dz A dy = 0 since u is harmonic. Thus the 
function v exists and is determined up to a real additive constant. O 


We see that the theories of harmonic and analytic functions in the plane 
essentially coincide, the former being the real parts of the latter. This fact 
determines the role played by analytic functions in mathematical physics. By 
use of these functions one can obtain the solutions of many problems, for 
example, airflow around a wing (a Zhukovskii airfoil). 
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Theorem 2 (the mean-value theorem). The mean value of a harmonic 
function over a circle equals its value at the center (see Fig. 8.2): 


27 


z fude = u(0) . 
20 
0 


Fig. 8.2. The mean-value theorem dy 


We shall prove the analogous theorem in the multidimensional case below, 
and we shall see that this property can be taken as the definition of a harmonic 
function. 


PROOF. As we know, u = Re f, where f is holomorphic. Consider the Cauchy 
integral formula 
1 ese, 
0) = — | — dt. 
F(0) 2ri / t ue 
Y 


This equality holds for any curve y enclosing the point 0. 

In the real-valued theorem we shall make explicit use of the fact that the 
curve y is a circle. We parameterize a point of the circle as t = Re'?, so that 
dt _ iRe® dy 

t Rev 
integral, we find 


= idy. When we substitute these values into the Cauchy 


2r 


HO = fDi, f=utiw. 


0 


We obtain the equality we need separately for the real and imaginary parts. 
The theorem is now proved. O 


In the multidimensional case the theories of harmonic and analytic func- 
tions diverge. 


8.1. Consequences of the Mean- Value Theorem 


1. The Maximum Principle. A harmonic function has no extrema in the 
interior of its domain. More precisely, if u(xo, yo) = max u(x, y) and (xo, Yo) € 
D, then u = const. 
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PROOF. Suppose the point (xo, yo) at which the maximum occurs is an in- 
terior point of the domain. Consider a small circle with center at that point, 
as in Fig. 8.3. If the value at some point (x), yi) of the circle is less than at 
(zo, Yo), then by continuity the function on a whole arc is less than at (zo, yo), 
violating the mean-value theorem. Hence the function must be constant in a 
neighborhood of the point (xo, yo). In other words, the set of points of the 
domain at which u(zo, yo) = maxu(z,y) is open. But this set is obviously 
closed. Thus the set of interior maximum points is either empty or coincides 
with a component of the domain. As a result the function must be constant 
on the entire component of the domain containing (zo, yo). O 


u(zo, Yo) > ulz f yı) 


Fig. 8.3. The maximum principle 


(r1, 41) 


This proof relies only on the mean-value theorem. Thus the maximum 
principle will have been proved in the multi-dimensional case as soon as we 
prove the multidimensional mean-value theorem. 

From the point of view of analytic function theory the maximum principle 
is a manifestation of the local open mapping principle. Indeed, if the real part 
of an analytic function could attain a maximum at an interior point, the open 
mapping principle would be violated, as shown in Fig. 8.4. 


l 

f (z0, Yo 
( ) Fig. 8.4. The maximum 
principle and the open 


l 
l 
l mapping principle 
l 


max Re f 


2. Uniqueness. The continuous solution of the Dirichlet problem in a 
bounded connected domain is unique. 


PROOF. The difference of two solutions is a harmonic function in the interior 
of the region and identically equal to zero on the boundary. Since the domain 
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is bounded, the maximum and minimum of that difference must be attained 
(either in the interior of the domain or on its boundary). If the maximum or 
minimum occurs in the interior, the difference is constant, and hence equal to 
zero. If it occurs on the boundary, then both the maximum and minimum equal 
zero, and so the difference is again identically zero. Thus by the maximum 
principle the difference is identically zero in the domain. O 


We remark that the uniqueness of the solution of the Dirichlet problem 
holds for continuous functions satisfying the mean-value theorem, that is, 
there is at most one continuous function satisfying the mean-value theorem 
and having prescribed values on the boundary of a bounded domain. 


Remark. In unbounded domains the solution is not unique. For example, 
the function x is harmonic in the half-plane x > 0, continuous on the closure 
of the half-plane, and identically equal to zero on the boundary, but not in 
the interior. 


Theorem 3. A continuous function satisfying the mean-value theorem is 
harmonic. 


PROOF. Consider the Dirichlet problem (in the ordinary sense, that is, for 
a harmonic function) in a domain with boundary values given by the values 
of a function satisfying the mean-value theorem. (We shall prove below that 
such a solution exists.) Consider the difference of the harmonic function and 
the function satisfying the mean-value theorem. It vanishes on the boundary 
and satisfies the mean-value theorem on the interior. Hence it must vanish 
identically. O 


Remark. The fact that the kernel of the Laplacian consists of very smooth 
functions — the real parts of analytic functions — is inherent in the nature of 
this operator. For the wave operator (d’Alembertian) the picture is quite dif- 
ferent. A singularity in the data for the Cauchy problem for the wave equation 
propagates along characteristics, so that the solutions of the wave equation 
cannot be smoother than the initial data. 


We now turn to the solution of the Dirichlet problem. We first construct 
the solution for a disk. 


OUTLINE OF THE FIRST METHOD. We assemble a rather large supply of har- 
monic functions in the disk. On the boundary of the disk we consider the 
function cosny. Is there an analytic function with this real part? There is, 
namely, 


fn(z) = (=) = A j (cosny + isin ng) 
R R á mye 
The real part of the function if,(z) on the boundary |z| = R is —sinny. 


Thus we can solve the Dirichlet problem if the boundary condition has the 
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form cosny or sinny. But we know that any boundary condition can be 
expanded in a series of these functions. Hence we can solve the problem with 
any boundary condition. 

We now carry out the first method in the Dirichlet problem for the disk. 
We expand the boundary condition in a Fourier series: 


a 
f(y) = = + N (an cosny + bn sinng). 


Passing to polar coordinates, we write the harmonic extention of the func- 
tions cos ny and sin nọ over the disk: 


r n r n : 
—} cosna, —]) sinna. 
eS a 
As a result we obtain the solution in the form of a series: 
r 


n n 
u(r,a) = > + ‘> (an (5) cos na + bn (=) sin na) (x) 
Our aim is to represent the solution in integral form. The Fourier coeffi- 


cients are 
2r 


on = = | f(B)cosnf ds 


o~ 
= 
|l 
ape 


f f(8) sin nB a8 , 


ne f f(8) dB. 


We substitute these relations into (*), obtaining 


u(r.) = = fro [E (F)" cosn(a — 8) + 5] 48. 
0 


where 
z= (r,a), (Z) cosn(a — 3) = Re E , t= Rel. 


We now sum the geometric series: 


bE ey 4 sap a 
TE 2 1-2 2 t-z 2° 


obtaining 
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+z  (t+2z)(t—2z) 


t 
We now transform the kernel: ae Gane). 


The real part of the numerator is tt — zz = R? — r?, since the remaining 
expression tē — tz is purely imaginary. 

The denominator |t — z|? equals R? + r? — 2Rr cos(a — 8). 

Thus we obtain the Poisson kernel 


R? — r? 
R? +r? — 2Rr cos(a — 8) ` 


To summarize: 


27 


1 Per 
une) == a a a Ale: 
0 


Problem. Prove that this integral can be obtained by a conformal mapping 
of the solution we shall obtain below for the half-plane (second method). 


By the Riemann mapping theorem any simply connected bounded domain 
can be conformally mapped onto a disk, and hence a solution of the Dirich- 
let problem can be obtained in any such domain using the Poisson kernel, 
provided the conformal mapping is known. 


OUTLINE OF THE SECOND METHOD. 


Problem. Let u be a solution of the Dirichlet problem in a domain D1, and 
let f : Dı — Də be a conformal mapping of the domains. Construct a solution 
of the Dirichlet problem in the domain Də. Try to convince yourself that this 
solution works without computing the Laplacian. 


We know that a disk can be conformally mapped onto a half-plane. Let us 
try to solve the Dirichlet problem for a half-plane: Au = 0, u(z,0) = g(x). 

We begin by noting that for the simple boundary condition shown in 
Fig. 8.5 the solution is the function +y,(z, y), where 


Palz, y) = Imln(z — a) = arctan 


Indeed, this function is harmonic in the half-plane and satisfies the boundary 
condition of Fig. 8.5. We now consider the “corner” function 


Fab (X, y) a ~ Cooley) Fa Pal, y)) 


Fig. 8.5. A boundary condition 
for the Dirichlet problem | 
in a half-plane 
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Fig. 8.6. The “corner” function 


(see Fig. 8.6). It is harmonic in the half-plane and satisfies the boundary 
condition g(x) = 1 fora < x < b and g(x) = 0 for x < a and q > b. 
We approximate an arbitrary boundary condition g(x) (say, continuous 


and equal to zero outside a finite interval) by a piecewise-constant function, 


Fupg(a)(x — a 
as in Fig. 8.7, and the corresponding harmonic function is ` Paglu) = a) 


b-a 
g(z) 
Fig. 8.7. Approximation 
by a piecewise-constant 
boundary condition 
a b r 


We now pass to the limit as the partition is refined, obtaining a harmonic 
function with boundary condition g(x) written in the following integral form 


fim LAGU da . 


b-a b—a 


The function lim Falt, y) 


—>a — 


in this integral is called the kernel. Let us study 


its level lines. The picture of the level lines for Fy,(x, y) is shown in Fig. 8.8. 
A conformal mapping leaving the point 0 fixed and mapping a to oo maps this 
picture to the simple picture of the level lines for the corner. After dividing 
by the length of the interval and passing to the limit, we obtain the picture 
of the level lines shown in Fig. 8.9. 


Q 
Fav(x, y) = z 


Fig. 8.8. Level lines 


—> for the function Fab 
a/ 7 
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Fig. 8.9. Level lines of the kernel 


It is clear that some computations are required here: the derivative of the 
function + arctan -+ with respect to the parameter a is 


d - 
— | — arctan 
T 


1 1 y 1 y 

da eal T y \ (z-a)? n(r-a} +y?’ 
1+ (=) 
r—a 
Finally, 
e 
_t y 
=i | goias) da. 
— 0 


The picture of the level lines for the kernel shown in Fig. 8.10 is well-known 
from physics — it is the picture of the equipotential lines of a so-called dipole 
(the lines of force of a dipole, which are orthogonal to the equipotential lines, 
form the same kind of picture). 


Fig. 8.10. A dipole 


8.2. The Mean-Value Theorem in the Multidimensional 
Case 


The Maximum Principle in the Multidimensional Case. We shall give 
a geometric explanation of the fact that a harmonic function cannot have 
a strict maximum at an interior point of its domain. A harmonic function 
minimizes the Dirichlet integral f (Vu)? dz. 

D 
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If there were a strict interior maximum, the graph would be bell-shaped 
in a neighborhood of that point, as shown in Fig. 8.11. Cutting off the top of 
the bell with a horizontal plane, we would obtain a function with a smaller 
Dirichlet integral. Thus a function with an interior extremum cannot minimize 
this integral. To be sure, this reasoning does not encompass functions with 
nonstrict maxima, as in Fig. 8.12. Nevertheless, we have proved that a function 
equal to zero on the boundary of a bounded domain and minimizing the 
Dirichlet integral vanishes identically on the domain. 


fae Fig. 8.11. A function that does not 


minimize the Dirichlet integral 


Fig. 8.12. A nonstrict maximum 


The full maximum principle follows from the multidimensional mean-value 
theorem. 


PROOF OF THE MEAN-VALUE THEOREM. On the sphere of radius R with cen- 
ter at 0 in R” consider the element of solid angle 92,,_1, that is, the projection 
of the element of surface area on the concentric sphere of radius 1, as shown 
in Fig. 8.13. We set f(R) = f u(r) Qn-1. Then 

|z|=R 


ð 
P= | Emai f (Vundan, 
lz|=R |c|=R 


Fig. 8.13. On the mean-value theorem 
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and this is the flux of the gradient of u across the sphere divided by R™~!. 
But if u is harmonic, that flux is 0. Thus f’(R) = 0, and so f(r) is constant. 
But for small values of R it is obvious that this integral is approximately u(0). 
Hence f(#) = u(0), which was to be proved. O 


Corollary. u(0) = a f  u(x)dz, where Vr is the volume of a ball of ra- 
dius R. R |z|<R 


Indeed, to integrate over the ball we can integrate first over a sphere of 
fixed radius and then with respect to the radius. 


Corollary. The maximum principle holds for a harmonic function in n-di- 
mensional space. 


Corollary. The Dirichlet problem has a unique solution in the n-dimensional 
case. 


Remarks. 1. When deriving the Poisson formula we did not verify conver- 
gence out to the boundary. Of course our reasoning shows that if there exists 
a formula for the solution of the Dirichlet problem, it must have the form 


2r 


1 R 7? 
u(r,a) = || wae dg . (+) 
0 


Nevertheless, it is still necessary to verify carefully that for every con- 
tinuous function f the function u defined by formula (**) indeed solves the 
Dirichlet problem. Try to carry out the necessary estimates independently. 
(In case of difficulty the reasoning can be found in the textbook of H. Car- 
tan Elementary Theory of Analytic Functions of One and Several Variables, 
Chap. IV, §4, par. 4.) 


2. Of course, the second derivation of the Poisson formula (for the half- 
plane) has the same gap as the first - we explain how to write the Poisson 
kernel, but we do not verify that the resulting solution works for every con- 
tinuous function on the boundary. 


3. In our second derivation of the Poisson formula we first solve the Dirich- 
let problem for very simple, but still discontinuous initial data. The solu- 
tions we found (connected with the “corner” function) are bounded. Only for 
bounded solutions can one guarantee the uniqueness of the solution of the 
Dirichlet problem with discontinuous boundary data, but we have not stated 
or proved this theorem. All this emphasizes the heuristic character of our 
derivation of the Poisson formula. 


4. Here is an example of the nonuniqueness of the solution of the Dirichlet 
problem for discontinuous initial data. Consider the Dirichlet problem in the 
half-plane y = 0 with initial condition u(z,0) = 0. There are solutions that are 
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discontinuous at oo. Consider any entire function f(r) that is real-valued on 
the real axis (for example, f(z) = z or f(z) = exp(z) and the like). Then the 
function u = Im f(z) gives a solution of our problem. We note that all these 
solutions, except u = 0, are unbounded. Prove this using analytic function 
theory. 

Prove also a more general fact: two bounded functions that have only a 
finite number of points of discontinuity in a closed disk, are harmonic on the 
interior of the disk, and coincide on the boundary coincide at all their points 
of continuity. 


TERY 
Ny 
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Lecture 9 


The Fundamental Solution for the Laplacian. 
Potentials 


The interaction of mathematics and physics is proceeding along strange paths 
at present. Dirac, one of the greatest physicists of the twentieth century, stated 
the strategy of theoretical physics as follows: “In setting about to construct 
a physical theory, one must reject all previous physical models, as well as the 
‘physical intuition’ based on them, which is nothing more than a collection of 
preconceived points of view.” We must, according to Dirac, simply choose a 
beautiful mathematical theory and systematically find physical interpretations 
for its consequences, not fearing any contradiction with previous theories. 


I learnt to distrust all physical concepts as the basis for a theory. 
Instead, one should put one’s trust in a mathematical scheme, even if 
the scheme does not appear at first sight to be connected with physics. 
One should concentrate on getting an interesting mathematics. 


P. M. Dirac 
(Cited in the book N. Wiener, by P. Massani, Birkhäuser, 1990, p. 6.) 


Amazing as it is, the whole theoretical physics of the twentieth century (in 
contrast to the physics of the nineteenth century) confirms Dirac’s correctness. 
It has been remarked that every new physical theory refutes all preceding 
theories, while the mathematical models remain. 

But today we shall take up a mathematical theory that, to the contrary, 
long remained unrecognized by mathematicians, although physicists had long 
made free use of it - the theory of so-called generalized functions. The most 
important example of a generalized function is the Dirac 6-function. We first 
study it on the line, then in R”; it could be defined on any manifold. This is 
the mathematical analog of such physical concepts as a point charge and a 
point mass. 

The physical “definition,” against which orthodox mathematicians protest, 
is as follows: the 6-function equals 0 everywhere except at zero, where it equals 
infinity. Its integral over the entire line is 1. Physicists have an excellent ability 
to work with such definitions, which have no explicit mathematical meaning. 
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We shall adopt the following definition, which is simple in applicatior 
makes it possible to explain the physical meaning easily, and is less tediot 
than the definition in the “correct mathematical theory of generalized func 
tions.” 

Our 6-function of the variable x is denoted d(x) and concentrated at ( 
When 6(x) occurs in a formula, it is to be replaced by a 6-shaped functio 
6-(x), and the limit is to be taken as £ tends to zero. 

A 6-shaped function is an ordinary smooth function with the followin 
properties: 


1. Its integral over the entire line equals 1. 

2. It is nonnegative. 

3. It is concentrated in the interval (—e¢,¢), that is, it vanishes outside thi 
interval. 


(These last two conditions could be weakened.) 
The typical form of the graph of a 6-shaped function is shown in Fig. 9.1 


Fig. 9.1. A 6-shaped function 


In passing to the limit one can fix the 6-shaped function. (The numbei 
€ occurs as a parameter in the formula.) One may also leave it unspecified 
provided the family 6,(x) satisfies conditions 1)-3). 


9.1. Examples and Properties 
1. f d(x) dz = lim f ðe (x) dz = 1. 
2. ô(x) = 0 when z # 0, since lim 6-(x) = 0 for x £ 0. 


3. Let us compute f f(x)ó(x)dz, where f(x) is a continuous function. 

The main contribution to the integral f f(x)ôe(x) dz is given by the expres- 
sion f f(0)d-(x) dz in Fig. 9.2, which tends to f(0) as € — 0; the difference 
between it and the whole integral also tends to zero as € — 0. (This is an 


exercise.) As a result, f f(x)d(x) dz = f(0). 
4. Similarly, f f(y)6(x@ — y)dy = f(z). 


The reasoning is the same, taking into account that ô- (x—y) is the translate 
of the function s(x) from a neighborhood of 0 to a neighborhood of y, as in 
Fig. 9.2. 


9.2. A Digression. The Principle of Superposition 79 


f(x) Aa 


x ly x 


Fig. 9.2. Computing the convolution with the 6-function 


The formula of item 4 also has the following remarkable interpretation: an 
arbitrary function f is a continuous linear combination of translated 6-func- 
tions 6(- — y) concentrated at the points y with coefficients f(y), so that the 
6-functions 6(- — y) form a “continuum size basis” in the space of functions 
on the line. It is also useful to interpret 6(-— y) as the density of a unit point 
mass or charge concentrated at y. 


5. Homogeneity. 


Definition. A function f is homogeneous of degree d if f(Ax) = At f(x) for 
any positive A. 


Let us find out whether the 6-function is homogeneous, and if so, of what 
degree. As an example, let us find (2x). Substitute 2x in place of x in the 
function ôe(x). The value of the integral is then cut in half. Passing to the 
limit, we find (2x) = 36(2). The reasoning is analogous for any multiple. 
Thus, the function (x) is homogeneous of degree —1, just like the function 
1/x. This similarity is not accidental. In a certain sense these functions are 
related, to be sure, only in the one-dimensional case. 

In the n-dimensional case the function d(x) is homogeneous of degree —n. 
One can also verify this as follows: 6(x1,...,%n) = 6(%1)---d(%n), since the 
corresponding product of 6-shaped functions gives a 6-shaped function in n- 
dimensional space. Since the factors are homogeneous of degree —1, the prod- 
uct is obviously homogeneous of degree —n. 


9.2. A Digression. The Principle of Superposition 


Although this principle is the basis of all of so-called linear physics, it is in 
essence a simple fact of linear algebra. 

Consider a linear transformation L : V — W mapping one vector space 
into another, and consider a nonhomogeneous linear equation Lu = f. The 
equation has a solution if and only if f € Im L, where Im L, the image of L, 
is a subspace of W. The solutions form an affine subspace in V parallel to 
Ker L, where Ker L, the kernel of L, is a subspace of V. To be specific, let up 
be a particular solution of the nonhomogeneous equation. Its general solution 
can then be represented as u = Up + Un, where up, is the general solution of 
the homogeneous equation Lu = 0, as shown in Fig. 9.3. 
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U = Up + Uh 
/ 


Fig. 9.3. Solutions of a nonhomogeneous linear equation 


L 4 Im L 
V W 


We now state the principle of superposition itself. Let u} and uz be re- 
spectively solutions of the equations Lu = fı, Lu = fo. Then œu, + ague is 
a solution of the equation Lu = a fı + Q@2f2. Here a, and ag are arbitrary 
constants. 


Physically, if we throw one stone into water we obtain a certain picture 
of the waves; if we throw a second stone, we obtain a second picture of the 
waves. If we throw both stones together, the resulting picture of the waves 
is the same as if we added the disturbances from throwing the first stone to 
those from throwing the second. 


The proofs of all the assertions listed above are obvious by the linearity of 
the operator. 


We now combine the principle of superposition with the possibility of 
expanding any function into a superposition of 6-functions. We obtain the 
following. If we want to solve a linear differential equation with an arbitrary 
function f on the right-hand side, it suffices to learn how to solve this equation 
with the 6-function on the right, then apply the principle of superposition, 
except that instead of sums one must write integrals. 


Let us apply these considerations to the Poisson equation Au = f. We 
must first solve the equation Au = 6(z). Of course, the solution of this equa- 
tion is not unique. Indeed, the kernel of A consists of all harmonic functions, 
so that an arbitrary harmonic function can be added to the solution. 


Of all the solutions of the equation Au = 6(x) one can exhibit one reimark- 
able one called the fundamental solution. Let us find the grounds on which it 
can be exhibited. 


The Laplacian is invariant with respect to the group of Euclidean motions 
(rotations and translations). Indeed, the Laplacian is the divergence of the 
gradient and hence is defined by the Euclidean structure alone. On the other 
hand, the 6-function is also spherically symmetric. Therefore a solution of 
this equation is again obtained by rotating any solution of Au = 6 around 
the origin. The arithmetic mean of two solutions is also a solution. In general 
any linear combination of several solutions whose coefficients add to 1 is again 
a solution. For that reason, one can average a solution over the group of all 
rotations about the origin, which is a compact group. 
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For example, when n = 2, the group SO(2) of rotations of the plane is 
isomorphic to the circle S!, and the averaging measure is so dy. 

When n = 3, the group of rotations SO(3) is isomorphic to the real three- 
dimensional projective space RP?. It has a two-sheeted covering by the three- 
dimensional sphere (the group of unit quaternions), which in turn is isomor- 
phic to the special unitary group SU(2), also known as the spin group of 
order 3, as in the following diagram: 


S? & SU(2) =Spin3 


>| 


SO(3) = RP? 


In particular, a metric can be carried over from S? to SO(3) under the two- 
sheeted covering. This metric is rotation-invariant, so that there is a measure 
with which one can carry out averaging. 

After averaging we obtain a solution that is rotation invariant. In other 
words, the value of the solution u(x) depends only on the distance from zx to 
the origin. For such functions the operator A becomes an ordinary second- 
order differential operator. For that reason the equation for u has the form 
u” + A(r)u'+ B(r)u = 0 for r > 0. (Here the right-hand side is zero, since the 
6-function vanishes for r > 0.) 

The coefficients A and B can be computed explicitly, if one does not mind 
passing to polar coordinates in the Laplacian, but one can also get by without 
any computations. First of all, we know the obvious solution u = 1. Hence 
B = 0, so that the operator has the form u” + A(r)u’. 

We now apply homogeneity considerations. The Laplacian maps homoge- 
neous functions into homogeneous functions of degree two less. Since 6(z) is 
homogeneous of degree —n, a homogeneous solution u of the equation Au = 6 
must be homogeneous of degree 2 — n. Consequently it must be proportional 
to 1/r"~?. In particular, when n = 3, it is proportional to 1/r. (We remark 
that when n = 2 this formula does not work — there is no homogeneous non- 
constant solution in this case.) 

Let us now compute the Laplacian: A(1/r”~?) = div grad(1/r"~). 

The gradient of our homogeneous function is a vector field directed toward 
the center along radii; its length equals the absolute value of the derivative 
with respect to r: 


(To be verified: when n = 3, the gradient of the potential 1/r is a Newto- 
nian force, inversely proportional to r?.) 

To compute the divergence we remark that it is rotation invariant. We 
compute the flux of the gradient field across two concentric spheres, as in 
Fig. 9.4. 
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2-n 


Fig. 9.4. Computation of div grad r 


The flux across a sphere equals the magnitude of the field multiplied by 
the surface area of the sphere. The area of a sphere of radius R is w,_;R"?, 
where wn—1 is the area of the unit sphere (4r for n = 3, 27 for n = 2). Asa 
result the flux is (2 — n)wn-1. In particular, for n = 3 it is —47z. 

Our field thus has the remarkable property of “incompressibility”: its flux 
is the same across any sphere centered at 0. A spherically symmetric field is 
unambiguously determined by this condition (up to a constant factor). For 
example, in three-dimensional space only a Newtonian field has this incom- 
pressibility property. On the plane there also exists an O(2)-invariant incom- 
pressible field: the length of a vector of this field is inversely proportional to 
the distance from the center. This is a gradient field, but its potential is not 
a homogeneous function; it is In L, We thus arrive at the remarkable formula 
A a jyt 
r? r 


Remark. This formula could have been guessed, and along with it the fun- 
damental solution of the equation Au = 6 in the plane, by considering the 
dimension of the plane n = 2 as the limit of the dimensions n = 2 + e: 


1 1 1 
= = = exp (cm2) =1l+eln—+4+O(e?). 
r r r 


pn—2 


Let us return to the computation of the divergence. The total flux of the 
field across the boundary of the layer between concentric spheres is zero. But 
this flux is the integral of the divergence over the layer. Thus 


T2 


Ve div grad(r?-") dr = 0 


rı 


for any rı and re, so that the divergence itself must be zero. Finally, 


1 , 1 
A (==) = div grad (=) =0 


for r > 0. But then one can also find the coefficient A(r). Thus we have 
avoided a lengthy change of variables. 


Claim. A(1/r”7?} 4 0 throughout space. 
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Indeed, consider the flux of the gradient field across a sphere with center 
0, which equals (2 — n)wn-1. Moreover, the sphere can be replaced by the 
boundary of an arbitrary region. If this region contains 0, the integral of 
the divergence over the region equals (2 — n)wn-1, while the integral is 0 if 
the region does not contain 0. But this means, by the very definition of the 
d-function, that the divergence of the field equals 


(2 — n)wn—14(2) . 


For example, when n = 3, we have A(1/r) = —476(x), and so forth. 

Of course this result can also be verified using é-shaped functions. But we 
see how much simpler it is to work directly with the ideal limiting object, the 
6-function itself. 

When n = 2, the function In(1/r) will serve as the fundamental solution 
(potential). Indeed, the derivative of this function with respect to r is —1/r, 
and from that point on all our reasoning goes through. 

We remark that on the level of fields the result is the same in all dimensions: 
the field is inversely proportional to r”~!. The case n = 1 can also be included 
in the general picture: the fundamental solution of the operator d?/dz? is the 
function |tz| with a suitable t, and the magnitude of the field is proportional 
to 1/r!—!, that is, it is constant, as in Fig. 9.5. 


Fig. 9.5. The fundamental solution for the Laplacian when n = 1 


We see that t = 1/2, so that the fundamental solution is |z|/2. It is in- 
teresting that the case n = 1 fits into the general picture even on the level 
of coefficients. To be specific, the general formula for a fundamental solution, 
which is i 


uo = = 
i (2 = Jun- r2 


becomes r/wo = r/2 when n = 1. Indeed, the surface area (zero-dimensional 
volume) wọ of the zero-dimensional sphere S°, which consists of two points, 
is 2. 

Thus, the general formula for the fundamental solution holds in all dimen- 
sions except 2, where the logarithm appears: ug = — +- In =. 


Remark. The appearance of the logarithm here has its own profound causes: 
the logarithm is a so-called associated homogeneous function. The logarithm 
is associated with power functions in exactly the same way as the adjoined 
vectors are associated with eigenvectors in the theory of Jordan cells. You have 
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already encountered this phenomenon in the theory of resonance when solving 
homogeneous linear equations with constant coefficients: when the eigenvalues 
become multiple, quasi-polynomials are adjoined to the exponentials. 


The physical meaning of the fundamental solution is the potential of a unit 
point charge at the point 0. The gradient of this potential is the field created 
by this charge. 

It is physically interesting to consider fields with n = 3. But fields in 
smaller dimensions can also be obtained from three-dimensional pictures hav- 
ing suitable symmetry. 

Consider, for example, a uniformly charged line in RÌ. Let us find the 
force created by this charge at a given point as the resultant of point charges 
smeared over the line. The result will be the same as in consideration of a 
plane passing through our point perpendicular to the charged line, and the 
whole charge must be placed at 0, as in Fig. 9.6. 


Fig. 9.6. A line generates the same field of charges in each plane orthogonal to it 
that a point charge generates 


WARNING. It is important to add forces, not potentials, since the correspond- 
ing integral for the potentials diverges. We obtain another remarkable and 
amazing formula: 


1 
2ln — + const . 


+00 
f dz 
— Oo 
(It is forbidden ask the value of the constant!) 


Similarly, if the charge is uniformly distributed on a plane, we obtain the 
same picture of forces on each line perpendicular to the plane that we would 
obtain from a point charge when n = 1. 

A similar reduction can be carried out starting with any dimension. 

We now consider a continuous function f(z) defined in a bounded domain 
§2. The solution u of the Poisson equation Au = f can be found by the 
principle of superposition, representing f as a superposition of translates of 
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the 6-function. Let up be the fundamental solution, that is Aug(x) = 6(z). 
In other words, ug(x) is the potential created at the point x by a unit charge 
located at the origin. If a charge of magnitude f(y) is located at y, it creates 
the potential f(y)uo(x — y) at x. 

The superposition of all these fields is the field with potential given by the 
formula 


ulz) = | fly)uo(e») du 
Q 
This function is called the Poisson integral. By the principle of superposition 


Au(e) = J F(v)6(e — y)dy = f(x) . 
N 


Therefore u(x) is a solution of Poisson’s equation. 


Theorem. Let the continuous function f(x) have support in a bounded do- 
main 922. Then there exists a solution of the Poisson equation Au = f given 
by the formula 


u(x) = | flyjuo(z — y) dy, 
/ 


where uo is the fundamental solution for the Laplacian. 


Once an explicit formula for the solution is written, it can be verified 
directly that the equation is satisfied, without resorting to ð-functions. 


Exercise. Prove the theorem directly. 


The solution is not unique, since any harmonic function can be added to it. 
A unique solution can be exhibited by imposing a condition at infinity. 

When n > 2 the fundamental solution vanishes at infinity, and hence the 
entire Poisson integral vanishes. This condition identifies the solution uniquely. 
Indeed, the difference of two such solutions would be a harmonic function 
vanishing at infinity. By the maximum principle, it would then have to be 
identically zero. 

When n = 2, a solution that is unique up to an additive constant can be 
singled out by the condition |u(x)| < C In |z]. 

Now let us consider charges concentrated only on the boundary 02 of the 
domain §2. We denote the surface charge density at a point y of the boundary 
by f(y). We set 


u(r) = J tujul — y) dy . 
an 
This function is harmonic everywhere outside the boundary of the do- 
main §2 and vanishes at infinity. (When n = 2, |u| grows no faster than the 
logarithm.) 
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What happens at points of the boundary itself? Consider an element dS 
of the boundary. We construct a cylinder over dS along the normal with small 
height £, as in Fig. 9.7. 


+ 


n 
Aaa e 
Cay z Fig. 9.7. The field created by a single layer 
n= 


Consider the flux of the gradient of u across the boundary of the region G 
just created. One would naturally expect the flux across the lateral surface to 
be small as £ is, since this lateral surface has small area. (We leave this point 
unproved for the present.) We first compute the total flux: 


On 


pedo = f(n,Wu)do= f Audr= f fas. 
ðG OG G 


(A2)NG 


(The second of these equalities uses Stokes’ theorem; the third takes ac- 
count of the charges in G.) On the other hand, the flux across the base is 
Cae + 2) dS up to a correction that is small compared with dS. Neglect- 
ing the flux through the lateral surface, which is small when € is small, we 
conclude that the equality T + oe = f holds at each point of 022. 

The sum on the left is the jump in the normal derivative if we regard the 
normal as oriented in one direction. Thus, the function u(x) is harmonic in 
the interior and exterior of the surface, and the jump of its normal derivative 
equals the density of the charge distributed on the surface. This function is 
called a single-layer potential — the charges are distributed in a single layer 
over the surface. 


Example 1. Let us compute the potential of a sphere with a constant charge 
density. Up to a constant factor we obtain u(x) = f rere as in Fig. 9.8. 
lWi=R 

This integral over the sphere can be computed explicitly, although the com- 
putation is complicated. (Newton was able to handle this.) 

But the answer can also be found without integration, using the symmetry 
of the sphere. We know in advance that u = u(r) is a harmonic function inside 
the sphere and outside it. Hence we can seek a solution in the form a + b/r. 
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y u 
Fig. 9.8. The potential Ar R 
of a uniformly charged ° ! 
sphere X i 
R r 


(The constants can be different inside and outside the sphere.) Inside the 
sphere b = 0, since if not, the function would not be harmonic at 0. Thus u 
is constant inside the sphere, and therefore equal to its value at the center, 
namely —4r R. The force field inside the sphere is therefore zero. 

As it happens, Newton discovered this by balancing forces, as shown in 
Fig. 9.9. His reasoning is also valid for an ellipsoid. More precisely, for an 
infinitely thin homogeneous layer between similar ellipsoids with a common 
center, where symmetry considerations (the dependence of the potential only 
on the radius) do not carry through. 


Fig. 9.9. There is no attraction 
inside a sphere or an ellipsoid 


Problem. Prove that the attractive forces of opposite elements of the layer 
between concentric similar ellipsoids are negatives of each other. 


Let us return to our sphere and find the constants a and b for the exterior 
region. On the outside a = 0 (the potential is zero at infinity). The coefficient 
b is the total charge and equals 4r R. One can also find a from consideration 
of the continuity of the function u at the points of the sphere itself under 
approach from within and without, since the integral converges uniformly. 
The graph of the potential is shown in Fig. 9.8. 


Example 2. Let A be a symmetric operator in Euclidean space and I the 
identity operator. The operator A — AI defines a one-parameter family of 
quadratic forms ((A — AJ)z,x), called a Euclidean bundle. The operator 
(A— AI)! (the resolvent) is also symmetric and defines the nonlinear bundle 
of quadratic forms ((A—AJ)~! 2, x) (dual to the Euclidean bundle in the sense 
of projective duality). Consider the family of quadrics ((A — AZ)~!z,x) = 1. 
For example, when n = 2, the bundle has the form of Fig. 9.10. All the curves 
here have the same foci. 
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Fig. 9.10. The curves of the bundle 
for n= 2 


The equations of the quadrics of the bundle can be written in an eigenbasis 
of the operator A as 


2 2 
Tí £ 
a gh da aera, tee 
ai- À Qn — À 
where a1,... , an are the eigenvalues of the operator A. Quadrics corresponding 


different values of À are said to be confocal to each other. 


Problem. Prove that precisely n mutually confocal quadrics of the bundle 
pass through each point of n-dimensional space (corresponding to some values 


Aj,.--,An). These quadrics are pairwise orthogonal at each point of their 
intersection. 
The values Aj,...,An are called the elliptic coordinates of the point. 


When n = 3, an ellipsoid, a hyperboloid of one sheet and a hyperboloid of 
two sheets pass through each point of space, as shown in Fig. 9.11. 


Fig. 9.11. The surfaces of the bundle 
for n = 3 


Problem. Consider an ellipsoid and an infinitely near concentric similar el- 
lipsoid. Consider the potential created by a uniform distribution of charge in 


the layer between these two ellipsoids. Such a potential is given by the den- 


dz, A\---Adzx > d . 
sity eae Atal that is, wee where ds is the Euclidean area and f the 
quadratic form that gives the ellipsoid. It is called a homeoid density. Then 


1. The electrostatic field inside the ellipsoid is zero. 
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2. On the outside the lines of force (phase lines of the vector field given by the 
gradient of the potential) are coordinate lines of the elliptic coordinates 
described above. 

3. All the level surfaces of the potential in the exterior are ellipsoids confocal 
with one another. 


We recall some simple propositions of electrostatics. 


The potential of a steady charge distribution of a conducting surface is 
constant on the surface. Indeed, if the potential were not constant the poten- 
tial gradient (that is, the electrostatic field) would have a component in the 
direction of the surface, which would cause the charges to move. 

The homeoid distribution just defined is the charge distribution on the 
surface of a conducting (“metallic”) charged ellipsoid. 

It is interesting that the charge density on a conducting surface is greater 
where the curvature of the surface is greater. Inside a conductor there cannot 
be any charges of a steady distribution. Indeed, the gradient of the poten- 
tial of a steady charge distribution equals zero inside a conductor (otherwise 
the charges would begin to move). Hence the potential inside a conductor is 
constant and consequently the charge density there is zero. 

In exactly the same way the potential of a smooth conducting surface 
bounding a bounded domain is constant in that domain if there is no charge. 

Indeed, the potential is harmonic inside the surface and constant on the 
surface. Consequently, by the maximum principle it is also constant inside the 
domain. 

Simply put, an electrostatic field cannot penetrate to the interior of a 
metallic cavity (“screening”). Analogous facts hold in R” for any n. 


Problem 1. A unit charge is placed at a point at distance r from the center 
of a noncharged conducting circle of radius R in R?. Find the electrostatic 
field created, its equipotential lines, and its lines of force. 


HINT. Start with the case r = 0. 


Problem 2. Carry out the computations for the density of an equilibrium 
distribution of charges for a “charged” conducting ellipse or square in the 
two-dimensional problem. What is the singularity of the charge distribution 
at the vertices? 


9.3. Appendix. An Estimate of the Single-Layer 
Potential 


We shall prove that the electrostatic field created by a smooth single-layer 
charge distribution on a smooth bounded surface is bounded all the way to 
the boundary. 


90 9. The Fundamental Solution for the Laplacian. Potentials 


Hence it follows in particular that the flux across the lateral surface of a 
cylinder of small height € (in the constructions on p. 86) is small compared 
with £, a fact we made use of in proving the formula for the jump of the 
normal derivative of the single-layer potential. 

Consider the family of normals to the surface of the layer. A sufficiently 
small neighborhood of this surface is fibered into segments of these normals. 
It suffices to prove that the force is uniformly bounded inside such an r-neigh- 
borhood (since outside the neighborhood the force is majorized by pS/r?, 
where p is the maximal charge density and S is the area of the surface). Con- 
sider a segment of the normal passing through a point P on the surface. In a 
neighborhood of the point P we introduce Cartesian coordinates with origin at 
P: z along the normal to the surface and x in the tangent plane. If the radius 
r of the neighborhood of P is sufficiently small, the equation of the surface of 
the layer in our coordinates can be written as z = h(x), |h(x)| < C|x|?. 

To prove that the force is bounded it suffices to obtain a uniform estimate 
of just the force created by charges in the sphere of radius r just constructed. 
(The force due to other charges, say in the r/2-neighborhood of P, does not 
exceed 4pS/r?.) 

A real difficulty occurs in the estimate of the force created near the point 
P by charges in the neighborhood. Here the boundedness is obtained only 
because opposite parts of the neighborhood pull in opposite directions. 

We begin by considering the special case when the layer in our sphere is a 
horizontal plane (z = 0) and the charge density is constant. 

In this case the horizontal components of the forces created on the z-axis 
by charges at opposite points x balance each other precisely. The resulting 
force is vertical and is given at the point Z as an integral (when the density 
is one): 


r/Z 
F f Z 2nz dx 5 f EdE on (1 Z 
(x? + Z2)3/2 (€2 + 1)3/2 Vr? + Z2 : 
0 0 
where x = |x|. (For simplicity we assume Z is positive; the force does not 


change when Z changes sign.) 

In any case the force is bounded by 27 (corresponding, as we already know, 
to the field of a uniformly charged plane). 

We now replace the charge distribution along the surface z = h(x) with 
density p(x) in the region z < r by a distribution of constant density p(0) in 
the disk x < r on the plane Z = 0. (Here Z is the coordinate of the point on 
the z-axis where we wish to estimate the magnitude of the field.) 

The vertical component of the force at Z is now given by the integral of 


(Z + h(x))p(x)J(x)x dz dy 
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where p is the density and J is the Jacobian, which is equal to the ratio of 
the element of area of the layer over x to the element of area x dz dy in the 
plane of x (with polar coordinates z, y). 

The quantities p and J are bounded in this region. In addition, the ex- 
pression in the denominator is 


— 


T? + (Z +h)? =x? + Z? + 2Zh(x) + h?(x) > 5 (e E) 
for sufficiently small r, since |Zh| < 7 + hk? and |h| < Cr?. 

For that reason the entire term proportional to Z in the integral can be 
estimated from above by a constant, just as in the case of a homogeneous 
plane. 

The second term (proportional to h) in the integrand is bounded. (Its 
numerator is less than Cız? and its denominator is greater than Coax.) Thus 
the entire vertical component is bounded. The horizontal component, that is, 

xp(x)J(x)xz dz dy 
(z? + (Z + A(x))?)3/?” 

First of all, we may assume that pJ = 1. Indeed, the functions p and 
J are smooth, so that the replacement of p(x)J(x) by p(0)J(0) changes the 
integrand by a bounded term (< C3z? in the numerator and > Coz° in the 
denominator). 

It remains to compare the integral for pJ = 1 with the same integral when 
h = 0 (which is equal to zero, as we know). 


the integral of is more difficult to estimate. 


Lemma. 


1 1 


(x? + Z2)3/2 (x? 4 (Z + h)?)3/2 


C4 


— r2 ® 


PROOF. Let Z = Ax. Then the left-hand side is 


1 1 1 
z3 | (1 +A2)3/2 (1+ (A+ p)?)8/? 


? 


where h = uz, so that |u| < Cz. 


The function ———-~-—~ has a derivative that is bounded by some con- 
(12) 37? 
stant Cs on the entire line. By the formula for finite increments 


< 
(1+ X?)8/2 (1 + (A + p)?)3/2 rile 


| 1 1 


Since |u| < Cz, the quantity being estimated in the lemma does not exceed 


which was to be proved. O 
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Finally, for the integral that gives the horizontal component of the force 
when Jp = 1, we obtain the estimate 


/ f GET KTP 


(The same integral equals zero when h = 0, and the difference of the integrands 
for h = 0 and for any h is estimated by the quantity C4/x? in the lemma.) 


P 
2 
zr C 
< 2r f =“ dr = 2rrC4 . 
T 
0 


Remark. By similar elementary estimates one can also prove that the single- 
layer potential is continuous in the entire space and that the limits of its 
derivatives along exterior and interior normals exist and are continuous, nor 
merely the formula for the jump in the normal derivative to which we have 
confined ourselves here. 


Remark. In the study of potentials there is a very useful convention due to 
Faraday that “4r electrostatic lines of force emanate from each unit charge” 
(whatever that may mean). It is assumed that the lines begin in positive 
charges (or at infinity) and end at negative charges (or at infinity). 

In this case the field intensity turns out to be equal to the “density of the 
lines of force.” For example, for single-layer charges half of the p4r dS lines 
of force emanating from the charges of the region dS go in one direction and 
the other half in the opposite direction, from which the formula for the jump 
follows. The factor 47 occurs here because physics considers the potential of 
a unit charge of the form u = 1/r, while we have defined the single-layer 
potential as the integral of translates of the fundamental solution of Au = ô. 


Lecture 10 
The Double-Layer Potential 


We have thus found the fundamental solution for the Laplacian in R”: 
1 


“i — Dwar? where w,,_, is the surface of the unit 


for n > 2, up = 
sphere in R”; 


1 
for n = 2, ug = —— ln =; 
F 


for n = 1, uo = —. 

Incidentally the signs are determined by the condition Aug = 6 as follows: 
when the singularities near zero are smoothed out, the smoothed function 
must have a positive second derivative. The coefficient is determined from the 
condition that the flux of gradu across the unit sphere must be equal to 1. 

Using fundamental solutions we constructed the following potentials: the 
solid potential, which gives a solution of the Poisson equation, and also the 
single-layer potential. We now consider a double-layer potential. 

We consider a charge distribution consisting of two layers on a hypersur- 
face, one layer positive and the other negative, with a certain density p(x) 
(p is a function on the hypersurface), as in Fig. 10.1. 


Fig. 10.1. A double layer T 
on a hypersurface 


Here we have in mind the following limiting passage. Let l be the distance 
between charges along the normal drawn to the hypersurface at the point q; 
the charges +Q = p/l and —Q = —p/l have magnitude of order 1/l; they 
create a certain field. We compute the limit of the field as | — 0. This limit 
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is called the dipole field, and its potential is the dipole potential. The dipole 
moment p = Ql is conserved in the limiting passage. The direction joining the 
infinitely near charges that make up the dipole is called the dipole aris. The 
potential of a double layer of density p distributed on a hypersurface is the 
integral of the dipole potentials of the charges located on the hypersurface 
whose axes are the normals to the hypersurface, the dipole moment in an 
infinitely small region ds on the hypersurface being pds. 

We now carry out the computation of the dipole potential in terms of the 
fundamental solution ug: 


(l)e eC) 


= pA + (1) Soleo: 
Ng 


Here n, is the outward normal to the hypersurface at q. 
We have found the element of the double-layer potential formed by a dipole 
located at q; altogether the double-layer potential at x is 


Ouo(x — q) 


d 
Ong q 


u(x) = i p(q) 


10.1. Properties of the Double-Layer Potential 


1. Geometric meaning of the double-layer potential. Instead of the 
fundamental solution normalized by the condition Aug = 6, it is custom- 
ary to use the function ug = 1/r”~? (when n > 2) and uo = In(1/r) when 
n = 2. The double-layer potential is a harmonic function in both the exte- 
rior and interior regions bounded by the surface, since this was true for the 
limiting passage 


Au == ():, 


R”\Ssr-t 


There may even be several exterior regions if the surface is not connected. 


Example. Let n = 2 and p = 1 on a connected closed curve. Then u(x) = 
const, the constant being different outside, inside, and on the curve. 


We shall now prove this remarkable fact. We first compute the dipole 
potential. The center of the dipole can be taken to be at the origin. The 
derivative of the fundamental solution in the direction of any vector v is 


2 Iv] cosy 
R 


1 1 
Lla- = —- Lyr = 
r r 
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Indeed, Lyr can easily be computed for any dimension: 


see + O), vir 
L, E E area) ok ET 


2yr? + +2 r 


as shown in Fig. 10.2. 


Fig. 10.2. The derivative of the radius-vector 
with respect to the vector v 


The derivative of the fundamental solution for any n > 2 is computed 

similarly: 
1 (2 — n)(r,v) 
pn-2 rn i 

We remark that the potential of an isolated charge decreases at infinity 
like 1/r?~?, while that of a dipole decreases faster, like 1/r”~!. (The positive 
and negative charges are said to “screen” each other.) 

Let us draw the level lines of the dipole potential (for n = 2). If we direct 


the x-axis along the vector v, we obtain the function = = + — 7: This is 
r tf +Y 

a very important harmonic function in the plane with the origin removed. It 

has a singularity at 0. Its level lines are conic sections given by the equations 

xz? +y? = cx. They form a family of circles passing through the origin and 

tangent to one another at the origin. These circles are perpendicular to the 


axis of the dipole at its center, as shown in Fig. 10.3. 


Fig. 10.3. Level lines 
of a dipole potential 


Problem. Draw the level lines of the potential of a pair of oppositely signed 
charges of equal magnitude on the plane. 


Hint. The difference of the logarithms of two numbers is the logarithm of 
their ratio. The geometric locus of the points whose distances to two given 
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points are in constant ratio (Fig. 10.4) is a circle (why?). As the charges 
coalesce (l — 0), the system of circles in Fig. 10.4 becomes the system of 
circles in Fig. 10.3 with a common tangent at the center of the dipole. 


Fig. 10.4. The level lines of the potential 
of a pair of opposite charges 


In our coordinate system a level line of the potential of a dipole, which is 
a circle, is tangent to the y-axis. If the curve at whose points we are studying 
the dipole potential is perpendicular to the axis of the dipole and has finite 
curvature at its center, it has second-order tangency with its osculating circle, 
which belongs to our family of level lines of the dipole potential. It follows that 
the limiting value of the dipole potential as a point approaches the center of the 
dipole along this curve passing through the center is completely determined 
and equal to the value of the dipole potential on the osculating circle, that is, 
it equals half the curvature of the curve. 


Corollary. The integrand in the integral that defines a double-layer potential 
restricted to the torus {x € S, q € S} can be extended by continuity even to 
the diagonal x = q (where the integrand is undefined). 


Consider an element of the curve and the angle it subtends at a given 
point x. The result is a function that is harmonic with respect to x at all 
points except the points of this element of the curve, as shown in Fig. 10.5. 


| x 
l 
ly Fig. 10.5. The angle subtended by an element 
of the curve is a harmonic function 
ds 
i . dssiny 7 ce 
Indeed, this angle is ——— = ds———, where x = rsiny and 
r tety. 
Y = r cos g. 


From this the geometric meaning of the (dipole) potential becomes clear: 
it is the element of angle subtended by an element of the curve when seen from 
a given point. Moreover, this holds in any dimension: the contribution of the 
neighborhood ds of the hypersurface to the integral computed at a given point 
x is the element of solid angle subtended by the element ds at the point x. 
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Problem. Prove that for any dimension the solid angle subtended by an 
element of the surface is a harmonic function. 


HINT. Use the formula for the derivative of the fundamental solution along a 
vector field. 


For n = 2 and n = 3 the dipole potential, defined as the derivative of 


1 1 
In — or — respectively, coincides exactly with the element of solid angle. In 
r r 


l l i i 1 
greater dimensions n the coefficient n — 2 appears if one uses ug = ere 


Actually the equality holds only up to a sign depending on the dioe of 
orientation of the hypersurface over which the integration of the angle element 
is taken. The hypersurface is oriented as the boundary of the oriented “inside” 
region. 

If p = 1, in the case n = 2 we obtain, after integration of the angle 
element, the “angle subtended by the boundary curve from the point x.” 
More precisely, with a given choice of orientation of the curve the potential is 
—2r times the number of revolutions of the vector from x to a point on the 
curve as the entire curve is traversed. When n = 3, we obtain similarly the 
oriented complete solid angle subtended from x by the surface, as shown in 
Fig. 10.6. 


Fig. 10.6. The complete solid angle 


This reasoning has replaced some rather lengthy computations. When 
n = 2, we obtain for a double-layer potential of density 1 distributed along a 
connected closed curve the value —27 in the interior, 0 in the exterior, and 
—m on the curve at each point where a tangent exists. 

When n = 3, we obtain respectively —47, 0, and —27. The values in any 
dimension can be obtained analogously. 


Remark. In the n-dimensional case the level hypersurfaces of the dipole po- 


= const. 


tential are given by the equation n 

The equipotential surfaces of a dipole have a tangent hyperplane at its 
center perpendicular to the axis of the dipole, pass through the center of the 
dipole, and are tangent to one another at this point, but they have a flattening 
at this point, as shown in Fig. 10.7. When n is even, the equipotential surfaces 
are analytic, even at the center of the dipole. When n is odd, the equipotential 
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Fig. 10.7. Equipotential surfaces 
of a dipole when n > 2 


surfaces have the same kind of singularity at the center of the dipole as the 
graph of the function |x|" of n — 1 variables. 

When n > 2, the curvature of every equipotential surface of a dipole equals 
zero at the center of the dipole. 

A surface in general position tangent to the equipotential surfaces at the 
center of a dipole will therefore have points of intersection with arbitrarily 
small equipotential hypersurfaces of the dipole when n > 2. For that reason 
our dipole potential is very large (in absolute value) on some of them. Hence, 
on a surface in general position passing through the center of a dipole and 
perpendicular to its axis the potential of a double-layer element will have a 
singularity that becomes infinite. 

Because of this, in general, the theory of the solvability of boundary-value 
problems for Laplace’s equation in the multidimensional case is not completely 
analogous to the Fredholm theory normally discussed in textbooks for the 
case n = 2. To be specific, when n > 2, the theory of integral equations with 
continuous kernels used in the Fredholm theory is no longer comprehensive. 
To be sure, the singularity is a rather weak one; it is integrable, so that one 
can apply the theory of operators with integrable kernel. 

Now suppose the density p is not constant. Let us investigate the properties 
of a double-layer potential. 


1. The harmonic property of u in the interior and exterior continues 
to hold. 


2. It is easy to find the difference in the values of u on the boundary of 
the domain and its limiting values as the boundary is approached from the 
interior and exterior. It is the same as it would be if the density were constant, 
namely (when n = 2): 


u_(Xo) = u(Xo) — Tp(Zo), U4 (Lo) = u(xo) + Tp(Zo) , 


where u_(Zo) is the limit from the interior and u,(z9) is the limit from the 
exterior, as shown in Fig. 10.8. 

To prove this we consider an auxiliary double-layer potential with constant 
density p = const = p(zo); the assertion is true for this density. We then 
consider the potential with density p—p(q), which equals 0 at q. We shall prove 
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u(xo) mp(xo) mp(zo) 


U4. (zo) —eoŇ 
u- (zo) u- (zo) u(zo) u4 (xo) 


Fig. 10.8. The jump of the double-layer potential 


that it has a zero jump. Instead of a formal proof, I shall explain why this is 
so. If the density were zero in an entire neighborhood of q, that neighborhood 
would not create a double layer. The rest of the surface would create a double- 
layer potential in that neighborhood, which is a harmonic function without 
singularities. Therefore the jumps at q would be zero. In fact the density is 
not zero in a neighborhood, only at the point q itself. For that reason, in order 
to prove that the jump vanishes, we need to carry out estimates analogous to 
those in the appendix to the preceding lecture. These elementary estimates 
show that the jump in the potential is zero. A sufficient condition for this is 
that the density be continuous at q. 


3. The normal derivative of the double-layer potential (the normal being 
oriented in the same direction for points in both the interior and the exterior) 
has no jump. 

This is certainly true in the case of a constant potential — the derivative is 
zero both from within and without. In the general case the derivative of the 
potential is the force (the “tension” of the field). Consider a cylindrical region 
N near the element ds of the surface, as shown in Fig. 10.9. Let us compute 
the flux of the force field gradu across the boundary of this region. As in the 
preceding lecture, we take account of only the normal component: 


i se 
On, On) 


M+ 
ðn 
grad u 
Fig. 10.9. Computation of the jump Q 
of the normal derivative of a na 
double-layer potential 


By Stokes’ formula, this flux is 


J (raduna) ds = J civ gradude = J Aude = const f Qax = 0, 
£2 2 


eRe] NQ 
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since this integral is the total charge in the region under consideration, and 
that is zero (the positive and negative charges cancel each other). Thus 


We shall use the properties of a double-layer potential below. At the moment 
we shall study the properties of the Laplacian. 


2. Oscillations of spherically symmetric bodies. The fundamental so- 
lutions have important applications to the problem of the vibrations of spheri- 
cally symmetric bodies. Consider, for example, a circle in the plane (a periodic 
string). The circle x? +y? = 1 is a Riemannian manifold. Hence the Laplacian 
div grad is defined on it. On the circle there is a standard polar coordinate g. 
In terms of this coordinate the operator A = div grad on the circle can be 
ð 

expressed in the usual way: A = Bie 

In higher dimensions the problem of a vibrating string is replaced by the 
problem of a vibrating sphere. There is no single coordinate here. Instead of 
the trigonometric functions (sin kọ, cos kọ) that describe the natural vibra- 
tions of a string, the vibrations of a sphere are described by the so-called 
spherical functions, which arise in all problems of mathematical physics hav- 
ing spherical symmetry. The place of the theory of Fourier series of functions 
on a circle is taken by expansions in spherical harmonics. 

From the practical point of view the problem of the gravitational potential 
of the Earth, which influences the motion of satellites, is of interest. Satellites, 
in turn, can be regarded as indicators for judging the distribution of mass in 
the Earth. As we shall soon see, the contributions of the higher-order har- 
monics decrease rapidly with distance from the Earth, so that at points not 
too near the Earth the potential is closely approximated by the sum of a 


C22 
spherically symmetric principal part — anda dipole perturbation =. 


Let us consider the problem of the eigenfunctions of the Laplacian on the 
n-dimensional sphere. For n = 1 these are the usual trigonometric functions. 
2 


d 
By solving the periodic boundary-value problem for the equation <a = Àu, 


one can find all the harmonics. These are just combinations of sin ky and 
cos ky with integers k. (The eigenvalue is \ = —k?.) 

What is the situation in higher dimensions, for example, on the two- 
dimensional sphere z? +y? +2? = 1? Here we become acquainted with spherical 
functions, which from this point of view are a generalization of the trigono- 
metric functions. 

Consider a function u(xz,y,z) that is constant along each ray emanat- 
ing from the origin. Such a function is said to be homogeneous of degree 0. 
(We recall that a function is homogeneous of degree k if u(Az) = A*u(z) 
for every A > 0.) 
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At the origin the function may be undefined. For example, one homoge- 
neous function of degree 0 is 


+9 see ee 
da 2) = 222 + 3y2 + 422 ` 
Problem. Find Au|,, and compare it with div grad(u| c>). 


This last operator is called the spherical Laplacian. I shall denote it by A. 
As it happens, the following identity holds in a space of any dimension n 
for a homogeneous function u of any degree k: 


Õu =r? Au — Au, A=k?+k(n—-2). 


The operator on the left-hand side of this equality is the spherical Laplacian, 
extended to the entire space by homogeneity of degree k. On the unit sphere 
r? = 1 it is simply the spherical Laplacian: 

Aul 52, = div grad (ul) . 


r= 
Let us now consider some special cases: 

1. For k = 0 we have the identity 

Au=r? Au. (10.1) 


This identity is easily explained. The gradient of a homogeneous function u of 
degree 0 restricted to the unit sphere coincides with the gradient of the restric- 
tion of the function itself to that sphere, since the gradient of a homogeneous 
function of degree 0 is tangent to the sphere. It is also easy to deduce by 
considering fluxes that the divergences of the gradients also coincide. That is, 
for homogeneous functions of degree 0 the usual Laplacian and the spherical 
Laplacian coincide on the unit sphere: (Au) aa = (Alaa The identity 
(10.1) can be obtained from the extension of the spherical Laplacian to the 
entire space preserving homogeneity. Indeed, the operator A lowers the de- 
gree of homogeneity by 2. It follows that Au is a homogeneous function of 
degree —2. In order for the extension to have degree 0, the function Au of 
degree —2 must be multiplied by r?. 


2. Let n = 2 and k £0. Our identity assumes the form 
Au = r? Au—k?u. 


Let us apply it to the function u = Re z*, which is harmonic in the plane and 
homogeneous of degree k. (In particular, when k = 2, u = cos 2y on the unit 
circle. ) : 

Our identity now assumes the form Au = —k?2u. 
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For the time being we shall not take up the proof of this general identity. 
We note only that it is based on Euler’s formula for homogeneous functions 
of degree k (whose proof is left as an exercise): 


Corollary. Suppose a function is homogeneous and harmonic. Then, re- 
stricted to the unit sphere, it is an eigenfunction of the spherical Laplacian: 


Au = —Au. 


We shall see below that all the eigenfunctions of the spherical Laplacian 
can be obtained in this way. 


Definition. The spherical functions on the sphere S”~! are the restrictions 
to the sphere of homogeneous harmonic polynomials in R”. 


Problem. Find the dimension of the space of spherical functions that are the 
restrictions of the harmonic polynomials of a given degree k of homogeneity 
in R°”. 


For example, for n = 3 we obtain: 


k 0 1 2 
dimension 1 3 5 
(basis x,y,z) (basis ry, yz, zz, 
Peay y e 


We shall now pause to study an interesting application of spherical func- 
tions to a topological problem. 


Maxwell’s Theorem. When n = 3, all spherical functions of a given degree 
k can be obtained by successively differentiating the potential 1/r along suitable 


1 
constant vector fields Ly,- Ly, — and restricting the result to the unit sphere. 
r 


(The result of differentiating in this way is called a multifield potential.) Here 
the k fields v1,..., Up are uniquely determined (up to nonzero factors) for a 
given nonzero spherical function of degree k. 


The dimension of the space of spherical functions of degree k is 2k + 1. 
By considering this space up to a nonzero multiplicative constant and omit- 
ting 0, we obtain the projective space RP?*. 

On the other hand, by Maxwell’s theorem all the spherical functions of 
degree k can be obtained, up to a nonzero multiple, from 1/r by differentiating 
along k constant vector fields of length 1 in RÌ. The constant vectors of length 
1 are uniquely determined up to sign by the function itself. The differentiations 
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along constant fields commute, so that the result is independent of the order 
of differentiating. 
Thus we have constructed a one-to-one mapping 


RP? x... x RP?/S(k) > RP” , 


where S(k) is the symmetric group of permutations of k factors. 

The first of these spaces is called the symmetric kt? power of the projective 
plane RP?, and is denoted S*RP?. The Maxwell mapping S*RP? — RP?* 
just constructed is a homeomorphism. Maxwell’s theorem seems to be the 
most elementary proof that these spaces are homeomorphic. 

A related theorem from algebra is Viete’s theorem. The Viete mapping 


roots — elementary symmetric polynomials in the roots 


defines a homeomorphism of the symmetric kt” power of C and C*: SFC = CF. 

The corresponding projective theorem is: S*(CP!) ~ CP*. The complex 
projective line CP! is the Riemann sphere S?, so that Viéte’s projective the- 
orem provides a homeomorphism S*(S?) ~ CP. 


Lecture 1] 


Spherical Functions. Maxwell’s Theorem. 
The Removable Singularities Theorem 


We consider a homogeneous function F of degree k in R”: 
F(A) =* F(x) VA>O. 


Such a function may be undefined at 0. (A homogeneous function may be 
defined in a homogeneous domain, for example, in some solid angle with origin 
at 0.) 

We have defined a modified spherical Laplacian A for each k. It maps 
homogeneous functions of degree k to homogeneous functions of the same 
degree. We recall its definition: a function is restricted to the unit sphere, 
then the divergence of its gradient is taken, and the result is extended to all 
of space except the origin as a homogeneous function of degree k. 


Theorem. 
AF =r?AF-AF, A=k*?+k(n-2). 


PROOF. We compute AF = div grad F. Let f be the restriction of F to the 
unit sphere: f = Flaen f: S”! — R. Then (grad EY css = grad f+ 3E 2 
as shown in Fig. 11.1. (Of course, we are assuming here that the tangent space 
to the sphere at the point is imbedded in the tangent space to the ambient 
space.) 


Fig. 11.1. The spherical and total gradients 


By homogeneity F(rq) = r* f(q), so that ÌE = kr*—! f(q). It suffices for us 
to find the gradient on the unit sphere, since it is homogeneous of degree k— 1. 
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Our field decomposes into tangential and normal components, the tangential 
component having no flux across the top or the bottom of a test surface and 
the normal component having no flux across the lateral walls, as in Fig. 11.2. 


Fig. 11.2. Computation of the flux 
through a test surface 


Let us compute the flux of the gradient across the test surface: 


J (era Fyn) dS = [ed F,n) dS tolte KIr tki, 
OG T 


where the first term is the flux across the lateral surface, the second is the 
flux across the top, and the third is the flux across the bottom. As a result, 
we obtain: 


[tered Fin)dS+oa(n+k — 2)ekF + o(e). 


T 


Here eo is the volume of the test region. The functions F and f are equal 
on the unit sphere. Dividing the flux by the volume of the test region and 
contracting the test region to a point on the unit sphere, we find that at this 
point 

AF =Af+kF(n+k—2). 


Indeed, the flux of grad F across OG is the integral of div grad F with respect 
to G, while the flux of grad F across T is £ x (the flux of grad f across 0o)+0(e). 
This last flux is the integral of div grad f over ø. 

Thus, our identity is proved on the unit sphere. When it is extended to a 
point at distance r from the origin, the F on the right-hand side is multiplied 
by r*, while AF is multiplied by r*~?. That is why the factor r? arises in the 
first term of the identity being proved. 

The modified spherical Laplacian is homogeneous of degree k by definition. 
Therefore if the identity holds on the unit sphere, it holds everywhere (except 
at the origin). Thus the identity is proved. O 


Our theorem is a special case of the following simple, but useful remark on 
the Laplacian of a function on a submanifold of Euclidean space. In computing 
the value of the Laplacian at a point one can replace the submanifold by its 
tangent space at that point and the function by the corresponding function 
in the tangent space. 
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Consider an m-dimensional submanifold (“surface”) 
y= f(z), ceER™”, yeR' 


in the Euclidean space R™+! with the metric dx? +dy?. We call the x-subspace 
horizontal and the y-subspace vertical. We shall take the Cartesian coordinates 
of x as the coordinates of the point X = (x,y), y = f(x), on the surface. 

A function U defined on the surface can be written in these coordinates 
as a function on the horizontal plane 


Lemma. Assume that the tangent plane to the surface at the point 0 is hori- 
zontal (that is, that (df/dx)(0) = 0). Then the Riemannian Laplacian of the 
function U at this point on the surface coincides with the Euclidean Laplacian 
of the function u on the tangent plane. 


PROOF. We write the metric on the surface as a Riemannian metric dS? 
on the horizontal plane. It follows from the Pythagorean theorem that the 
difference dy? between this metric and the Euclidean metric in the horizontal 
plane is a second-order infinitesimal in |z]: 


dS? — dx? = O(|z|?). 


(Here and below smallness of a quadratic form or differential operator means, 
of course, that the coefficients are small.) Therefore the operators GRAD and 
DIV of the Riemannian gradient and the Riemannian divergence correspond- 
ing to the metric dS? on the horizontal plane differ from the Euclidean grad 
and div by infinitesimals of second order: 


GRAD - grad = O(|z|?) , DIV -div = O(|z?]) . 
Applying these operators successively to the function u, we find 
DIV GRAD u = div grad u + O(|z|?) + div O(|z]”) . 
This last term is O(|z]|), so that 
(DIV GRAD u)(0) = (div grad u)(0) . o 


Remark. The translation of a function to the tangent plane does not neces- 
sarily have to be done using orthogonal projection: one can use any family of 
smooth curves. All that matters is that the curve passing through the point 
is orthogonal to the surface at that point. 


Example. Consider the unit sphere y = v1 — 2? in n-dimensional Euclidean 
space and the point z = 0 on it (m = n — 1, l = 1). Let us consider the 
function U on the sphere and compute its spherical Laplacian at x = 0. 
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By the lemma it equals the Euclidean Laplacian A(x)u (the sum of the 
second partial derivatives of the corresponding function u(x) with respect 
to Ti). 

We shall denote the homogeneous extension of U of degree k to Euclidean 
space by U. On the tangent plane y = 1 to the sphere at the point x = 0 this 


function becomes 
5 k - . 
a(x) = (VIF z?)"u(3), |č- z| = O(ļzļ). 


Therefore the Euclidean Laplacian of the extension at the point x = 0, 
y= lis 


2 d2 
AU = ae y" u(0) + (Azii)(0) = k(k — 1)u(0) + Ai. 
But @(x) = u(x) + £ 2?u(0) + O(|z|?). Therefore at the point x = 0 


A,(t)+ A,ut+k(n—-1)u. 
Finally, at the point in question, 
AU = Azu + [k(n — 1) + k(k — 1)Ju(0) . 


Since the point chosen differs in no way from the others, we have proved 
our formula AF = r?AF — AF, A = k(n +k — 2) for r = 1. By homogeneity 
it holds everywhere. 


Corollary. If the function F (defined in R” \ {0}) is harmonic and homoge- 
neous of degree k, it is an eigenfunction for the (modified) spherical Laplacian: 
AF = -AF (that is, its restriction to the sphere is an eigenfunction of the 
spherical Laplacian). Conversely, an eigenfunction of the spherical Laplacian, 
when extended from the unit sphere as a homogeneous function of degree k, is 
harmonic everywhere except at the origin. 


We remark that for a given eigenfunction of the Laplacian on the unit 
sphere one can find two harmonic extensions, since the quadratic equation 
A = k? + k(n — 2) determines two values of k. 


Example. Let n = 2, that is, we are considering functions on the plane 
and their restriction to the circle. We have A = k*. The function zx is a 
harmonic homogeneous function of degree 1. But there also exists a harmonic 
homogeneous function of degree —1, whose restriction to the unit circle is the 
same as that of the function z. 

In polar coordinates the formulas of these conjugate functions: F = 
reosp = z, F = (cosy)/r = x/r?, form a dipole potential (the kernel of 
a double-layer potential). Similarly, for the conjugate exponents of homogene- 
ity k and —k: 
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Let n = 3, that is, we are considering functions of three variables and their 
restrictions to the two-dimensional sphere. We have A = k? + k. Therefore the 
exponent conjugate to k is Å = —1 — k. : 

For k = 0 the conjugate functions are F = 1, F = 1/r. 

For k = 1 we have F = z, F = z/r®. As in the two-dimensional case, this 
is a dipole potential (the kernel of a double-layer potential). 

For k = 2 we can find a harmonic quadratic form, for example, F = 
pee a ee) Bl aS 5 bet 1), 


The functions z and f are spherical, that is, they describe the natural 
vibrations of a sphere. Moreover, they are so-called zonal functions: they are 
invariant under rotation about the z-axis and reverse sign when crossing the 
parallels that divide the sphere into separately vibrating zones, and the di- 
viding parallels remain fixed, as shown in Fig. 11.3. 


z z 
Fig. 11.3. Zonal functions ey es 
describe the natural = 
vibrations of a sphere Eee, = 
eee 


Problem. Does there exist, for every integer k, a zonal spherical function 
that is the restriction to the sphere of a polynomial of degree k in z? 


In the general case A = k? + k(n — 2), and the homogeneity exponent 
conjugate to k is k = 2 — n — k. The values of the conjugate functions at 
the points rq (q being a point of the unit sphere and r the distance from the 
origin) are F(rq) = r* f(q) and F(rq) = f(q)/r®t"~?. 

These values can be connected using an inversion. When this is done, the 
image of qr is q/r, as shown in Fig. 11.4. 


Fig. 11.4. Inverse points q 


We have: F(q/r) = f(q)/r*, so that F(rq) = F(q/r)/r"~2. (We remark 
that when n = 2 the inversion continues to be harmonic, and the coefficient 
drops out of the last formula.) 
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For any n the exponent of homogeneity k has completely disappeared 
from the formula connecting F(q/r) and F (rq). Hence this formula holds 
for linear combinations of homogeneous functions with different exponents 
of homogeneity. But any harmonic function can be approximated with such 
functions. Thus under this transformation any harmonic function F maps to 
a harmonic function F(x) = F(x/|z|?)/|x|"~?. (This can also be verified by 
direct computation. ) 


Proposition. There exist homogeneous harmonic functions of degree k only 
for integers k. (Eigenfunctions of the spherical Laplacian exist only for A 
corresponding to integer k.) 


PROOF. We begin by studying homogeneous harmonic polynomials. The 
space of homogeneous polynomials of degree zero is one-dimensional, those 
of degree 1 have dimension n. We shall find the dimension of the space of 
polynomials of degree 2 and higher in the three-dimensional case for the sake 
of simplicity (see Fig. 11.5): 


degree 0 12 3 
dimension 1 3 6 10 
y? 
ry’ Fig. 11.5. Determination of the dimension 
of the space of homogeneous polynomials 
ey of degree three 
r’ 2? 
; k+2)\(k+1 . 
In general dim S*R? = E +2 = CELLE The Laplacian maps this 


space to the space S*~?R? of homogeneous polynomials of degree 2 less than 
the degree of the polynomials of the original space. The harmonic homoge- 
neous polynomials of degree k form a vector space, namely the kernel of the 
mapping 

Avi aS eR. 


Claim. The Laplacian maps S*R? onto S*~?R3. 
Indeed, any polynomial is the image of some polynomial under the Lapla- 
cian. (It suffices to verify this for monomials. For monomials in one variable x 


this is obvious: 0?u/Ox? = x? when u = x?+?/(p + 2)(p + 1). For the mono- 
mial «?y? in x and m additional variables y; we begin with the monomial 
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u = rPt?y7/(p + 2)\(p + 1). We have Az, yu = z?y? +v, where v = Ayu is a 
polynomial whose degree in y is less than |q|. Therefore z?y? will be in the 
image of Az, if all polynomials of degree less than |q| in y belong to that 
image. This makes it possible to prove what is required by induction on |q]. 
(If Az yw = v, then Az plu — w) = r?y?.) 
Thus the dimension of the space of spherical functions of degree k equals 
the difference of the dimensions of S*R? and S*~?R?, that is, 
EFIT I RSD ppt a. o 
2 2 
Corollary. For each nonnegative integer k there exists a (2k+1)-dimensional 
vector space of spherical functions on the sphere S*. These are eigenfunctions 
of the spherical Laplacian with eigenvalue —A, where A = k? + k, and they 
are the restrictions to the sphere of harmonic homogeneous polynomials of 
degree k in R3. 


PROOF. The fact that the restrictions are eigenfunctions follows from the 
identity just proved. The dimension of the space of harmonic homogeneous 
polynomials of degree k in R has just been computed. The dimension of the 
space of their restrictions to the sphere is the same, since a homogeneous 
polynomial that vanishes on the sphere is identically zero. C] 


We have now found the dimension of the space of harmonic homogeneous 
polynomials of degree k for n = 3. In the general case the dimension grows 
like k™~?. 


Claim. Spherical functions with different eigenvalues are pairwise orthogonal. 


Indeed the natural vibrations with different natural frequencies are always 
orthogonal, like the principal axes of an ellipsoid with different lengths in 
Euclidean space. 

Hence it follows, for example, that f (3z? — 1)dz = 0. Of course, this 

2 


integral can be computed explicitly using symmetrization: the three integrals 
for x, y, and z are all equal, and their sum equals 0. 


Theorem. Every spherical function (eigenfunction of the spherical Lapla- 
cian) is the restriction to the sphere of a homogeneous harmonic polynomial 
in the ambient space. 


PROOF. The eigenvalues of the spherical Laplacian are nonpositive, since the 
potential energy — the Dirichlet integral — is nonnegative. 

Eigenfunctions that are harmonic on the sphere correspond to the eigen- 
value 0. A function that is harmonic on the sphere is constant (by the maxi- 
mum principle), since the sphere is a closed connected manifold (compact and 
without boundary). Indeed, let us cut a small hole out of the sphere. Then 
the maximum will be reached on the boundary of the hole. Contracting the 
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hole to a point, we verify that at every point on the manifold the value of a 
function harmonic on a closed connected manifold equals the maximum of the 
function on the whole manifold. 

Now let us consider any negative eigenvalue —A, where A > 0. For any 
A > 0 there exists an exponent of homogeneity k > 0 such that A = 
k(k +n — 2), as shown in Fig. 11.6. 


Fig. 11.6. The eigenvalue as a function 
of the degree 


2-n k 


An eigenfunction of the spherical Laplacian with eigenvalue —A can be 
extended to the entire space with 0 removed as a harmonic homogeneous 
function of degree k (by our fundamental identity). 

Since k is positive, the extended function is bounded in a neighborhood 
of 0. By the removable singularities theorem to be proved below, a bounded 
harmonic function in a neighborhood of 0 remains harmonic when it is ex- 
tended to the origin by continuity. For the extended function we have u(0) = 0, 
since u(rq) = r®u(q), k > 0. 

The extended function is harmonic at 0 and consequently infinitely dif- 
ferentiable. But a homogeneous function of degree k that is infinitely differ- 
entiable at 0 is necessarily a homogeneous polynomial of degree k (and k 
is necessarily a nonnegative integer). This can be seen by Taylor’s formula. 
Thus, any spherical function can be extended to a homogeneous harmonic 
polynomial in the entire space. O 


Maxwell’s Theorem. When n = 3, all spherical functions of degree k can 
be obtained by successive differentiation of the function 1/r along k directions. 


ProoF. First of all, it is obvious that the derivative of a harmonic function 
along a constant vector field is harmonic. For example, the partial derivatives 
of a harmonic function are harmonic. l 

Second, on the sphere the derivatives coincide with certain homogeneous 
polynomials. This is true for the original function 1/r: that is, 1 is a homoge- 
neous polynomial of degree 0. 

Suppose that after a differentiations the result is a function of the form 
F/r°, where F is a homogeneous polynomial of degree d. When we differentiate 
along a constant vector field the homogeneity is preserved, while the degree 


11. Spherical Functions. Maxwell’s Theorem. Removable Singularities 113 
decreases by 1. Let Vy be differentiation along the constant field v. Then 


F VF -rt—Fa-r*-?(v,r) (VF)? — aF(v,r) 
Oat eg ee 

The numerator here is a homogeneous polynomial of degree d+ 1. Therefore 
the resulting harmonic function, after being restricted to the unit sphere, co- 
incides with a homogeneous polynomial. This polynomial is harmonic. Indeed, 
after differentiating we obtain a harmonic homogeneous function of negative 
decree. Its exponent of homogeneity is conjugate to the exponent of homo- 
geneity of the numerator, which is a homogeneous polynomial of positive 
degree. When a harmonic function is restricted to the unit sphere, the result 
is a spherical function. If we extend it with the conjugate positive exponent 
of homogeneity, the result is a polynomial and harmonic in the entire space 
(by the preceding theorem). Thus, the restriction of the derivative in question 
to the sphere coincides with the restriction of a homogeneous polynomial that 
is harmonic in the entire space. 

It can be proved (see Appendix A) that the space of functions produced 
in this way is a vector space (which is not at all obvious). 

To prove that all spherical functions can be obtained in this way it suffices 
to show that no eigenspace of the space of spherical functions can be invariant 
with respect to all rotations. Then the subspace obtained by the Maxwell 
construction (which is invariant under rotations) must coincide with the entire 
space of spherical functions. 

To prove that the representation of the group of rotations of the sphere S? 
by linear transformations of the (2k + 1)-dimensional space of spherical func- 
tions is irreducible, that is, that there are no nontrivial invariant subspaces, 
it suffices, for example, to verify that the entire space consists of linear com- 
binations of a single function and its rotations. In the role of such a spherical 
function that generates the whole space we can take, for example, a zonal 
spherical function. Let us study these functions in somewhat more detail. 

There are two ways of extending spherical functions to homogeneous har- 
monic functions. Let us consider the case of the Maxwell construction, in 


which all differentiations are carried out with respect to the same direction, 


, a k 
for example, along the z-axis. Restricting the function (2)" + to the sphere, 


we obtain the spherical functions. When k = 0, we have 1; when k = 1 we 
obtain z/r?; each time we shall obtain expressions of the form 


(ea r? — oF) ‘ae . 
Oz 


If F = F(z,r), this form of F is preserved under these differentiations, so 
that a function of z alone results on the unit sphere. Hence, in the space of 
spherical functions of any degree k there is a function depending on z alone. 
It is the restriction to the sphere of a certain polynomial in z, which is called 
a Legendre polynomial. Legendre polynomials of different degrees are pairwise 
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orthogonal on the interval [—1,1] of the z-axis, since the integral over the 
sphere of their product, which equals zero, transforms easily into an integral 
with respect to z, by Archimedes’ theorem that a sphere without poles and 
a cylinder are symplectomorphic. (The element of area of the sphere equals 
the element of area of the circumscribed cylinder under projection onto the 
cylinder along horizontal radii, as shown in Fig. 11.7.) O 


ea. 
fos 
/ = -1 
l l ' lz 
alaa a E E 


Fig. 11.7. Archimedes’ symplectomorphism of a sphere onto a cylinder and the 
graphs of the first three Legendre polynomials 


a 
XR 


Corollary. The Legendre polynomial of degree k has k distinct roots in the 
interval (—1,1). 


Indeed, the number of zeros is at most k, since the polynomial has degree k. 
Assume that the number of roots is less than k. Then the number m of roots 
on (—1, 1) at which the polynomial changes sign is also less than k. We form a 
linear combination of Legendre polynomials of degree at most m having a zero 
of order one at these m zeros. There exists such a linear combination, since the 
linear combinations of Legendre polynomials of degree at most m constitute 
all polynomials of degree at most m. (The space of polynomials has degree 
m + 1, and the first m + 1 Legendre polynomials are linearly independent.) 
We multiply this linear combination by our Legendre polynomial of degree k. 
The product does not change sign on the interval (—1,1), which contradicts 
the orthogonality of the Legendre polynomial of degree k to all Legendre 
polynomials of lower degree. 

Thus the number of distinct roots of the Legendre polynomial of degree k 
on the interval (—1, 1) is k. Consequently these roots are simple roots. 

Every Legendre polynomial describes a zonal natural vibration of the 
sphere, as shown in Fig. 11.8. | 

Besides the Legendre polynomial f of degree k, same eigenvalue also has 
the associated functions f;(z)cosjy, f;(z)singy (j < k), which also describe 
natural vibrations, as in Fig. 11.9; each piece of the grid vibrates indepen- 
dently. 

The associated functions can be obtained from the Legendre polyno- 
mial by rotating the sphere about the z-axis or the y-axis through a small 
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LEA 


Fig. 11.8. The vibrations of a sphere 
corresponding to a Legendre polynomial 


ne 


Fig. 11.9. The vibrations of a sphere 
corresponding to the associated functions 


angle œ. The Legendre polynomial thereby becomes a nearby spherical func- 
tion. An associated function measures the difference between the original 
function and the translated function. To be specific, the derivative of the 
translated function with respect to @ is an associated spherical function: 
(2(0/0z) — z(0/Ox)) f(z) = zf'(z) = f'(z)V1 — 22 cosy. Repeated differ- 
entiation produces cos2y and so forth. These computations prove, by the 
way, that all spherical functions are generated by a single zonal function and 
consequently prove that the representation is irreducible (has no invariant 
subspaces). 

In particular, it follows from this that the Maxwell construction produces 
all spherical functions. 

It can be proved (see Appendix A below) that the space of kt” derivatives 
of the function 1/r with respect to k constant (and consequently commuting) 
vector fields is a vector space. (The derivative is a linear function of the field.) 
The dimension of this vector space of derivatives is at most 2k + 1, since a 
nonzero kt" derivative is determined by the directions of k vectors (k points 
on the 2-sphere) and one common coefficient. 

On the other hand, we constructed above a linear transformation from 
this vector space of dimension 2k + 1 onto the whole (2k + 1)-dimensional 
space of harmonic homogeneous polynomials of degree k in RÌ. Thus this 
transformation is an isomorphism. Consequently the Maxwell construction 
not only gives all spherical functions with a given k, it gives them each exactly 
once, up to a multiple. This proves the following topological proposition: 


S*(RP?) = RP? | 
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A theory of spherical functions exists in any dimension, and there are 
zonal functions in every dimension. When n = 3, these are the Legendre 
polynomials. Let us see what results when n = 2. 

The vibrations of a circle are described by the eigenfunctions 


cosky, sin kg. 


The existence of zonal vibrations means that a linear combination of these two 
functions can be formed that depends only on z. But that is obvious: such a 
function is the function z = cosy, and also the function cos kọ for any k, as 
in Fig. 11.10. 


Fig. 11.10. Zonal 
vibrations of a circle 


Corollary. cos ky = Tk(cos y), where Tẹ is a polynomial. 


These polynomials are called the Chebyshev polynomials. They are or- 
thogonal on the interval [—1,1], with a weight, however (compute which 
weight) and have zeros with properties analogous to those of the Legendre 
polynomials. 

Incidentally, the role of the associated functions here is played by the 
functions sin ky resulting from zonal functions by differentiation with respect 
to the angle of revolution of the zonal function. 

The majority of the so-called special functions of mathematical physics 
arises in the problem of the vibrations of bodies with various symmetries in 
a space of some dimension. 

For example, the problem of the vibrations of a disk lead to the so-called 
Bessel functions. Bessel was a German astronomer who encountered these 
functions in Fourier-series expansions of the gravitational potential of the 
Sun restricted to a Keplerian ellipse. 

The fact that such important physical problems as the problem of vibra- 
tions of a disk and the attraction of planets lead to the same mathematical 
theory is a marvelous manifestation of the universality of mathematics and 
the unity of all reality. 

To conclude the proof Maxwell’s striking theorem, it remains only to prove 
the (also striking) removable singularities theorem used in it. 


e 
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Theorem. Let n = 2. If a function in R” is harmonic in some punctured 
neighborhood of a point and bounded, it can be extended by continuity to that 
point and the extended function is harmonic in the entire neighborhood of the 
point, including the point itself. 


This theorem has a simple physical meaning. A harmonic function de- 
scribes the stationary state of a membrane stretched over some curve. For 
example, the membrane could be stretched onto two hoops. The theorem as- 
serts that the membrane cannot be propped up at one point by a needle. The 
needle would pass right through, as in Fig. 11.11. 


Fig. 11.11. A membrane can be propped up by a hoop, but not by a needle 


In the case of a thin plate, which is described by the equation A?u = 0, 
the corresponding assertion is not true, that is, a plate can be supported at 
one point. 


PROOF. Consider a circle with center at the point 0 in question. (By trans- 
lating and dilating coordinates we can assume that the circle is defined by the 
condition r = 1.) We construct a harmonic function in the disk that coincides 
with the given function on the boundary. We already know that such a func- 
tion exists. (An explicit function is given by the Poisson integral, exhibited in 
Lecture 7.) We must prove that the difference between the original function 
and the constructed function, which we denote u, is zero in the punctured disk. 

Consider the function wo = C ln(1/r). It is harmonic in the punctured disk 
and equal to zero on the circle r = 1. We choose C so that this function is 
greater than u when r = €. 

The function u is bounded in the punctured disk; suppose |u| < M 
there. So that the function uo will be equal to M when r = e we choose 
C = M/I\n(1/e). Then u = up = 0 on the unit circle and u < up on the circle 
of radius £. Both functions are harmonic in the annulus between the circles 
and u < uo on both boundaries. By the maximum principle the function ug — u 
is nonnegative everywhere in the annulus, as shown in Fig. 11.12. Thus 


pere E, 
=  In(1/e) 
The function —u also satisfies such an inequality. 

Let £ tend to 0. Then we obtain u = 0 at all points of the disk, as in 
Fig. 11.12. O 
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Fig. 11.12. The barrier method of removing a singularity 


An analogous proof carries over to the case when the singularity at zero 
is weaker than logarithmic. The method just described for removing the sin- 
gularity is called the barrier method. The essence of the method is to find 
a positive harmonic function away from a singularity tending to infinity as 
the singularity is approached. By multiplying such a “barrier” function by a 
suitable constant, we can make it as small as desired at any point away from 
the singularity. If the function being studied tends to infinity more slowly 
than the barrier function as the singularity is approached, the singularity is 
removable: on the boundary of an annulus the function in question is less (not 
greater) than the barrier function, so that it is less (not greater) than the 
barrier function everywhere and hence equal to zero. 

The removable singularities theorem holds in any dimension; the singu- 
larity must be weaker than that of the fundamental solution. The proof is 
exactly the same as in the case n = 2, except that C'ln(1/r) must be replaced 
by C(r™-2 — 1). 


Corollary. A bounded harmonic function in the plane is constant. 


Remark. This result, due to Liouville, also holds in R” for any n, but the 
simple proof given below is valid only for n = 2. 


PROOF. By use of stereographic projection the Riemann sphere can be cov- 
ered by two flat charts related to each other by the transformation w = 1/z. 
Consider the function u(z) in the chart w, where it is harmonic and bounded in 
a neighborhood of zero. By the removable singularities theorem it is bounded 
and harmonic on the entire Riemann sphere and hence constant. O 


Remark. It suffices to have boundedness from only one side. For example, 
every positive harmonic function on the whole plane is constant. 


Problem. If a harmonic function on the plane grows no faster than a poly- 
nomial, then it is a polynomial. 


HINT. The derivatives of such a function also grow no faster than polynomials. 
Hence, by the Cauchy—Riemann equations it follows that the holomorphic 
function f whose real part is the original harmonic function grows no faster 
than some polynomial. Therefore the function f(z)/z% = w™ f(1/w), which 
is holomorphic in some punctured neighborhood of z = œ, (that is, w = 0) 
is bounded. It follows from the removable singularities theorem that it is 
holomorphic at w = 0, so that f is a polynomial of degree at most N. 
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Problem. Does a function that is harmonic and bounded in the complement 
of a line segment in R’ have a removable singularity? 


HINT. Take the potential of a charge of density 1 on the line segment as the 
barrier function. 


Problem. Prove Liouville’s theorem in R”: a bounded harmonic function is 
constant. 


HINT. Use Poisson’s formula or an expansion in spherical functions on a large 
sphere. 


Lecture 12 


Boundary- Value Problems for Laplace’s 
Equation. Theory of Linear Equations 
and Systems 


12.1. Four Boundary-Value Problems for Laplace’s 
Equation 


Consider the compact connected smooth hypersurface S"~! in R”, which di- 
vides R” into two regions: the interior (bounded) region G and the exterior 
(unbounded) region G’. Suppose a continuous function f : S"~! — R is given 
on the boundary. The Dirichlet problem for Laplace’s equation is to find a 
function u in the closure of the region G (G’) for which the following condi- 
tions hold. 


1. The function u is harmonic: Au = 0 in the region G (the interior problem) 
or G’ (the exterior problem). 

2. The function is continuous in the closure of the region: u € C(G) (resp. 
u € C(G')). 

3. The function satisfies the boundary condition u| gun = f- 

4. In the case of the exterior problem certain additional conditions are con- 
sidered at infinity; these conditions affect the existence and uniqueness of 
the solution. They are differently posed in different textbooks. The most 
common condition (which looks strange at first sight) is 


jul<C asrt=>œ, n=2; 


u-0 asr—-o, n>2. 


The Neumann problem is this problem with Condition 3 replaced by 
ax = f; in that case in Condition 2 it is necessary to require C -smoothness 
or some other sufficient condition for the normal derivative to exist. 

In the case of the Neumann problem conditions like Condition 4 are usually 
imposed as well. Physically the Dirichlet boundary condition (Condition 3) 
corresponds to a membrane clamped along its rim and the Neumann condition 
2u — 0 to a free membrane. 


ðu 
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By combining the interior and exterior versions of the Dirichlet and Neu- 
mann problems we obtain four types of boundary-value problems. Let us take 
each of them up in turn. 


1. The Interior Dirichlet Problem 


There exists a unique solution to this problem. 

For example, when n = 2, the Poisson integral gives the solution of the 
Dirichlet problem for the disk. By the Riemann mapping theorem any region 
can be conformally mapped onto the disk, so that it is possible to obtain 
a solution of the Dirichlet problem there also. Unfortunately, the proof of 
the Riemann mapping theorem is based on the existence of a solution of 
the Dirichlet problem; and besides, this technique does not apply in higher 
dimensions. Nevertheless, the result holds in any dimension: There exists a 
unique solution of the Dirichlet problem. 

This result remains valid in the case of a planar or multidimensional region 
bounded by several connected surfaces. The uniqueness follows immediately 
from the maximum principle (see the preceding lectures). The idea of the proof 
of existence is to minimize the Dirichlet integral under prescribed boundary 
conditions. This minimum can be achieved, as we already know, only at a 
solution of the Dirichlet problem. It is in fact attained, as is physically obvious. 
(It is the stationary state of a membrane stretched over the boundary.) But 
the proof is nontrivial and subtleties are encountered. For example, for a 
membrane propped up at one point the lower bound of the Dirichlet integral 
is not attained (see the removable singularities theorem in Lecture 10). We 
shall not prove that the minimum exists. 


2. The Exterior Dirichlet Problem 


This problem can be reduced to the interior problem. To do so one must 
perform an inversion with center at some point of the region. If this point is 
taken as the origin, the inversion is given by the formula x ++ x/|z|*?. When 
n = 2, harmonic functions transform to harmonic functions. For n > 2 the 
function u(x/|xz|?)/x2"~? will be harmonic. We study the new function in the 
bounded region into which the exterior of G maps under the inversion, as 
shown in Fig. 12.1, and pose the interior Dirichlet problem for this function. 

One can see that a singularity may arise at zero. But if the original function 
is bounded outside of G, the singularity is removable. Thus there exists a 
unique bounded solution of the Dirichlet problem. When n > 2, this solution 
tends to zero at infinity, since the bounded solution is divided by the quantity 
|z|"~? when it is returned to the exterior region, and that quantity grows at 
infinity. 
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Fig. 12.1. Transformation of a region under inversion 1 


An unbounded solution is not unique. For example, when n = 2, if the 
boundary is the circle r = 1, one can add any multiple of the fundamental 
solution (which grows at infinity and equals zero on the boundary) to the 
bounded solution. One can also add other harmonic functions, for example 
rcosy, (cosy)/r and the like. If we omit Condition 4, we obtain an infinite- 
dimensional space of solutions. 

When n > 2, Condition 4 is chosen so that the transformed function 
u(x/|x|?)/x”7ž, when divided by the fundamental solution 1/|z|"~? tends to 
zero as |x| — 0. The removable singularities theorem is then applicable to 
the transformed function, and the exterior Dirichlet problem reduces to the 
interior problem. That is how the condition u — 0 as x — 0 is obtained. 


3. The Interior Neumann Problem 


Suppose first that n = 2. A harmonic function is the real part of a holomorphic 
function. The imaginary part of this holomorphic function is the harmonic 
conjugate of the real part. The Cauchy—Riemann equations imply that the 
derivative of the unknown function along the normal to the boundary equals 
the derivative of the conjugate function in the orthogonal direction, as shown 
in Fig. 12.2. 


Fig. 12.2. Reduction of the Neumann problem 
to the Dirichlet problem 


We thus find that the derivative of the conjugate function v in the direction 
tangential to the boundary is prescribed. If v itself could be recovered from 
this derivative along the boundary, the Neumann problem for u would have 
been reduced to the Dirichlet problem for v. But it is not always possible to 
recover v from its derivative along the boundary, since topological obstacles 
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may intervene. For example, if the boundary is topologically a circle with 
coordinate s, a necessary and sufficient condition for v to be recoverable from 
f = ðv/ðs is that $ fds = 0. In the case of a nonsimply connected region, 


an annulus, for set it remains true that only one condition is required 
for solvability: the integral of the boundary function over the boundary of the 
annulus must be equal to zero. But there are not two conditions asserting 
that the integral over each circle must vanish, as one might have expected. 
(It is instructive to think about why this is so!) Thus, the following result is 
obtained. 


Theorem. The interior Neumann problem is solvable if and only if 
¢ fds = 0. The solution is unique up to an additive constant. Thus 
dG 


the solution space is one-dimensional. 


Remark. Consider a linear transformation A: X — Y. In the space X there 
is a subspace Ker A, the kernel of A, and in Y there is a subspace Im A, the 
image of A. 

The number k = dim Ker A is the number of linearly independent solutions 
of the homogeneous equation Au = 0, and the number s = codim Im A is the 
number of independent conditions one must impose on the right-hand side of 
the nonhomogeneous equation Au = f to guarantee solvability. 

We also use the notation c = codim Ker A and i = dimImA. We have 
c+k=dimX and i + s =dimY. In addition the equality c = i holds. (The 
quotient space over the kernel is isomorphic to the image.) Thus k — s = 
dim X — dim Y. 

We see that the number k—s, called the index of A is actually independent 
of the transformation and is determined by the spaces alone. 


Sometimes this observation can be generalized to infinite-dimensional 
spaces. This can be done provided A is sufficiently well approximated by 
finite-dimensional transformations. In our applications X and Y will be func- 
tion spaces, and we shall require sufficiently good approximation by finite 
sums of Fourier series. 

Let us examine our boundary-value problems from this point of view. 

For the interior Dirichlet problem s = 0, k =0,s—k=0. 

For the interior Neumann problem s = 1,k = 1, s—k=0. 

In both cases, it turns out to be a question of mapping the same spaces. 

The necessity of the condition $ fds = 0 for solvability of the interior 


ðG 
Neumann problem can also be obtained starting from considerations of min- 
imizing the Dirichlet integral. For n = 2 it is obvious, for example, that the 
normal derivative cannot be positive on the entire boundary without violating 
the maximum principle, as shown in Fig. 12.3. If the derivative is positive, a 
singularity that is at least logarithmic arises, as shown in Fig. 12.3. 
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Fig. 12.3. The normal derivative of a harmonic function cannot be positive 


Problem. Integrate by parts in the Dirichlet integral f (Vu)? dz in the mul- 
G 


tidimensional case and verify that the condition $ fds = 0 is necessary for 
OG 
solvability of the interior Neumann problem. 


4. The Exterior Neumann Problem 


In the two-dimensional case the exterior problem can be reduced to the interior 
problem. (This cannot be done in the multidimensional case.) We remark that 
under the condition 0u/On = 1 on the boundary circle r = 1, the interior 
problem cannot be solved, but the exterior problem has the solution In(1/r). 
Any function on the boundary can be represented as the sum of a constant 
and a function with zero integral over the boundary. By inversion the exterior 
problem can be reduced to the interior problem: the derivative along the 
exterior normal becomes the derivative along the interior normal. This is true 
for an arbitrary region: inversion is a conformal mapping, and the normal to 
the boundary line becomes the normal to its image under inversion. 

The final result is as follows. Without Condition 4 the solution is not 
unique; under Condition 4 the solution is unique and exists if and only if 


$ fds =0. 
8G 


The same is true in higher dimensions. Here s = 1, k =1,k—s=0. 


12.2. Existence and Uniqueness of Solutions 


The proofs of all these results on existence and uniqueness of the solutions of 
boundary-value problems are constructed by the following method. We seek 
a solution in the form of a potential with an unknown density. An integral 
equation is obtained for this density and shown to be solvable. 

Consider, for example, the interior Dirichlet problem. We seek a solution 
in the form of a double-layer potential with an unknown density p on the 
boundary. We denote the value of the potential on the boundary by Ap, where 
A is a linear (integral) operator on the space of functions on the boundary. 
By the potential jump theorem for the double layer, the limiting value of 
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the potential inside the region is Ap + Ap, where X is a constant depending 
only on the dimension. (It equals 7 in the two-dimensional case and 27 in the 
three-dimensional case.) The Dirichlet boundary condition u| aq = f therefore 
assumes the form of the equation (A + AF)p = f for the density p. 

If the space of densities were finite-dimensional, we would have the so- 
lution immediately: p = (A + AE)~!. The existence of the inverse operator 
is guaranteed then by the absence of a nonzero solution in the homogeneous 
equation. 

The uniqueness of the solution of the equation (A+AE)p = 0 follows from 
the maximum principle, since a function that is continuous in the closure of the 
region G and harmonic in its interior and vanishes on the boundary must be 
identically zero. Therefore if the space of functions p were finite-dimensional, 
the proof of the existence of a solution of the Dirichlet problem would be 
proved. 

Actually this space is infinite-dimensional, but the operator A can be 
approximated well enough by finite-dimensional operators to preserve the 
conclusion. 

The reason the operator A is nearly finite-dimensional is as follows. The 
operator maps a density p to the value of a potential on the boundary. 
Consider a point z on the boundary, as shown in Fig. 12.4. 


Fig. 12.4. A high harmonic makes only a small 
p contribution to the value of the potential 
T at a boundary point 


We expand the density in a Fourier series. A high harmonic of this se- 
ries oscillates rapidly. Physically this indicates the presence of nearby charges 
with opposite signs. They act oppositely to each other, so that such a har- 
monic makes only a small contribution to the result at the point x. Therefore 
the operator A can be well approximated by finite-dimensional operators. As 
a result, the index of the operator A + AE equals zero, as for an operator 
on a finite-dimensional space, although this operator acts on the infinite- 
dimensional space of densities. l 

In the multidimensional case for a spherical boundary one can use expan- 
sions in spherical functions instead of Fourier expansions. For other boundary 
manifolds local expansions analogous to Fourier series can be used near each 
point. 

With that we conclude our brief discussion of the questions of existence 
and uniqueness of solutions of the basic boundary-value problems for Laplace’s 
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equation. A detailed proof of the results stated here requires the machinery 
of the theory of Fredholm integral equations, which we do not have the time 
to develop in these lectures. A brief exposition of this theory can be found, 
for example, in the textbook of G. E. Shilov An Introduction to the Theory of 
Linear Spaces. 


12.3. Linear Partial Differential Equations and Their 
Symbols 


Let us return to the general theory of differential equations. I begin by recalling 
the general concept of linearization. Let us start with an ordinary differential 
equation given by a vector field v in the phase space. Let x9 be an equilibrium 
point of the differential equation t = v(x), that is, v(zo) = 0. Then in a 
neighborhood of zg one can pose the problem of small oscillations about the 
equilibrium position, which are described by the linearized system: y = Ay, 
where A = gu | as shown in Fig. 12.5. 


T=) 7 


Ay 


Fig. 12.5. Linearization 
in a neighborhood of 
an equilibrium point To yY x 


Replacing the equation t = v(x) by the equation y = Ay is called lineariz- 
ing the system at xo. The linearized vector field “lives” in the tangent space 
to the original phase space at the equilibrium point. It is independent of the 
coordinate system used in computing the Jacobian matrix ge ee 

In the case of partial differential equations the problem of small oscilla- 
tions of a continuous medium (a string, a membrane, and the like) originates 
in this way. It is the linearization of the corresponding equation in the dynam- 
ics of continuous media, whose phase space is the infinite-dimensional space 
of positions and velocities of points of the medium. In linearized problems 
involving small oscillations of a continuous medium the role of the operator A 
is played by linear partial differential operators, which may be of arbitrarily 
high order. 

The unknown function may be vector-valued rather than scalar-valued. 
(Picture, for example, a string in three-dimensional space. Its deviation from 
the equilibrium position at each point is a vector having two components, 
which are the coordinates in the direction perpendicular to the string.) 
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To describe such situations in general one must consider the vector bun- 
dle whose base is the manifold of independent variables. The values of the 
unknown functions (“fields,” as the physicists call them) define a section of 
this bundle. In the local theory one can write everything in coordinates and 
assume that x € R” and u € R”. In physics the components of a vector u are 
called fields; here r is the number of fields and n = D + 1, where D is the 
physical dimension. The extra 1 is reserved for time. 

A linear differential operator in the coordinate system (21,...,2n) has the 
form of a polynomial in the operators 0/Ox, with coefficients depending on z, 
that is, the form 


P(,2) = Y ag(2)d*, 


|a|=0 
where 


la| 
a = (@i, Aan) EZI, lal=art+---+a,, = o Ta 
The number m = max |a| is called the order of the operator. Such an operator 
acts on scalar-valued functions (r = 1). If there are several fields (r > 1), then 
“matrix-valued” differential operators occur with values in the vector bundle 
over the same base as the original fields. The dimension of a fiber is the same 
as the number of original fields, namely r. That is, the number of equations 
equals the number of unknown functions. In local coordinates the operator 
u |= v can be written as a system of equations: 


vı = P (8, £)ui +--+ + Pip(O,z)ur , 


Ur = P, (8, 2)uy +.: ++ P02) a; ; 


If we denote the whole matrix of operators P (8, x£) by P(O, x), we obtain the 
matrix notation v = P(ð, x)u. Here P is a matrix of polynomials in 0/0z, 
with coefficients that depend on z. 


Example. The Cauchy-Riemann equations have the form 


TE O, —Oy u\ 
Oy, ór v 
Here r = 2, the matrix is skew-symmetric, the degree of the polynomials is 
m = 1, and the coefficients are constant. 


Consider, more generally, the important special case when the differen- 
tial operators P are translation invariant with respect to z, that is, all the 
coefficients a(x) are constant. 
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The origin of this case is as follows. If it is known that the fields vary 
steadily, in first approximation on the local level one can apply the operation of 
freezing the coefficients, that is, passing to equations with constant coefficients 
equal to the values of the given coefficients at the point in question. For slowly 
varying fields freezing gives a rather good approximation. 

Thus, let us consider operators with constant coefficients: A = P(0). As it 
happens, such equations can be solved, and a purely algebraic theory results. 
Suppose for simplicity that r = 1. (We are considering a single equation with 
one unknown function.) 

Consider the translations T; : R” — R”, Tx = x + s. The operators 
A = P(ð) and T act on functions and commute with each other: AT, = T, A. 
Hence they have common eigenvectors. What eigenfunctions do the transla- 
tions have? Let us assume that our functions are complex-valued. 

Functions on the circle can be represented by Fourier series. As the segment 
rolled into a circle grows longer, the frequencies occurring in the representation 
grow closer together. To be specific, for a segment of length L the frequencies 
are k = 2rl/L, as shown in Fig. 12.6. 


Fig. 12.6. Closely packed 


frequencies occur in the ees ~nre HHH HH ome 
Fourier expansion 0 L 


One should expect that as the line segment is stretched to the entire line 
the series would become an integral. Then the imaginary exponents form 
a “continuum size basis”: e'**, and k ranges over the dual space (R”)*. 
(In physics k is called a wave vector.) The functions ek = eF? are called 
(harmonic) planar waves; they are eigenfunctions of the translation operator 
and differentiation: 


Tek = Ake ; Ak = elke f Oek = (ik) “ex ; 


We see that Ae, = P(ik)e,, that is, each harmonic planar wave is an eigen- 
function of any linear differential operator with constant coefficients, the eigen- 
value being a polynomial in the wave vector. 

The computations are analogous in the case of r fields: after multiplying 
the exponential by the amplitude w € R”, we obtain P(0)we, = (P(ik)) weg. 
Here P(ik) is a matrix whose entries are polynomials. 


Exercise. Find this matrix for the Cauchy-Riemann operator. 


f ik; —ikə 
ANSWER. ce ie | 


For any wave vector k and any differential operator P(O) with constant 
coefficients there are, in general, r eigenvectors w, of the operator P(ik). 
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Each of them defines a harmonic planar wave with wave vector k of the form 
w ex, which is an eigenvector of the matrix differential operator A = P(0). 
Physically, planar waves w;e, with a given wave vector k have common fronts, 
but different “polarizations” wj € R”. 

When r = 1, we have a single differential operator P(O), and the poly- 
nomial P(ik) in the wave vector k is called the symbol of this operator. The 
highest-degree term of the symbol is called the principal symbol. 


For example, the symbol of the Laplacian in R” is —k? — --. — k? and 
coincides with its principal symbol. The symbol is usually denoted by the 
letter ø. For example, for the Laplacian o = —k?. 


For a matrix operator the matrix of symbols is called the matrix symbol 
of the operator. 

The symbol of the system is the determinant of the matrix symbol. For 
example, the symbol of the Cauchy—Riemann is 


iki —ikg\ 2 12 
det (i 2) = ki — k5, 


that is, it coincides with the symbol of the Laplacian on the plane. 

The principal symbol of a system is the principal homogeneous part of 
its symbol. In wave physics the corresponding concept is usually called the 
dispersion relation. 


Remark. The principal symbol is invariantly determined, that is, it is inde- 
pendent of the coordinate system used to construct it. It is a homogeneous 
polynomial in the components of the cotangent vector k of the space of inde- 
pendent variables. Invariance of the principal symbol means that its values on 
each cotangent vector is independent of the coordinate system. The complete 
symbol is not invariant. For example, terms of first order can be added after a 
change of variables in the Laplacian. One can see that the principal symbol is 
invariant from the following considerations. (For simplicity we limit ourselves 
to the scalar case r = 1.) We apply the operator to the planar wave ex: 


P(8,x)ek = Om(k) + Om-1(k) +e, 


where the coefficients of the polynomials may depend on z. Let k — oo. If 
Om # 0, this is the principal term of the asymptotics of the value of the opera- 
tor on a high-frequency planar wave. Making a smooth change of variable x in 
ex (x), it is easy to verify that this asymptotics depends only on the vector k, 
and not on the coordinate system. 


The algebraic properties of the principal symbol exert a decisive influence 
on the properties of the corresponding differential equation, and in this sense 
the theory of differential equations can be regarded as a branch of algebraic 
geometry. 
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Example. Suppose the function f on the manifold can be expanded in a 
series of exponentials (for example, on the circle or the torus): f = $ freg. 
We seek a solution of the equation Au = f in the form of a series u = ẹ` ukek. 
We obtain —k?ug = fk, uk = —fe/k?. It is obvious that the series Y` ugeg 
converges; moreover the smoothness of u has increased by two orders compared 
with that of f. 


The same reasoning goes through for any linear operator of order m with 
constant coefficients if |am(k)| > C|k|™. 

Such an operator is said to be elliptic. We see that the infinite-dimen- 
sionality of the problem actually disappears due to ellipticity. The solution 
method is purely algebraic. 

An operator with variable coefficients is elliptic if it becomes elliptic under 
freezing of the coefficients at each point. 


Example. The Laplacian div grad on any Riemannian manifold (for example, 
on the surface of a sphere), is elliptic. (Prove this!) 


For an elliptic operator P(0,x) with variable coefficients the preceding 
formulas do not give an exact solution of the equation Pu = f. However, 
by freezing the coefficients, one can obtain a very good approximation to the 
solution of the elliptic equation Pu = f on a compact manifold. 

Choosing f to be “wave packets” like f = e~ $8" e(z), one can use the 
approximate solution u = f/om(k) when k is large. The rapid decrease of u 
as k — oo is guaranteed because the problem is “almost finite-dimensional.” 
Because of this property the theory of elliptic equations and boundary-value 
problems is nearly as close to finite-dimensional linear algebra as the corre- 
sponding theory for the Laplacian, with which we became acquainted earlier, 
even in the case of variable coefficients. 

If the operator is not elliptic (for example, the wave operator), this method 
does not go directly through, but it can be modified. 


Problem 1. Find the manifold of zeros of the principal symbol of the wave 
equation 


u o?u n Ou 
Ot? ðr?  ðy2` 
ANSWER. This is a quadratic cone in R3 (the “light cone”). 


Problem 2. Suppose the quadratic form $` apgkpkq is positive definite. For 
which values of the velocity c and the wave vector k does the wave equation 


O7u o?u 
Z? > ar OLpOLg 


have a solution in the form of a planar wave e’*~** moving in the direction 
of the vector k with velocity c? 
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Problem 3. Find the manifold of zeros of the principal symbol of the system 
of Maxwell equations 
a) in a vacuum: 


ee = curl H , oF = — curl E, 


ot ot 


b) in a homogeneous nonisotropic medium. 


HINT. See the book of Courant and Hilbert, Equations of Mathematical 
Physics. 


Problem 4. Find the manifold of zeros of the principal symbol of the system 
of Dirac equations for four complex-valued functions u; of the four variables z;: 


4 
X an(0/Or% — aku — pbu =0, 


k=1 
where 
0001 0 0 0 =i 
[001 0 00 i 0 
a“a=jioirio 0]? 2=1o -i 0 0]? 
100 0 i 0 0 0 
0 Or 1 0 1 0 0 0 
jo 0 0 =-1 {o ie 0 090 
eama a MOF o MSR FOS + eat. Qi ls 
0 -1 0 0 0 0 0 -1 
10 0 0 
01 0 0 
b=] o 9 -1 9 
io 0 -1 


The general outline of the theory of propagation of waves defined by sys- 
tems of linear partial differential equations is as follows. 

The set of zeros of the principal symbol is a field of cones in the cotan- 
gent spaces of space-time. Passing to the projectivization (or spherization), 
we obtain the “Fresnel hypersurface” of retardations in manifolds of contact 
elements of space-time. l 

The geometric optics (that is, the theory of rays and fronts) defined by 
this hypersurface (see Lecture 2) describes approximately (under some condi- 
tions) short-wave asymptotics (also called quasiclassical in physics), that is, 
the behavior of waves whose lengths are small in comparison with the geomet- 
ric dimensions of the system, for example, compared with distances at which 
the coefficients of the system vary noticeably. 
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A fundamental condition is that of hyperbolicity, which consists of the 
following. Consider an algebraic hypersurface of degree d in m-dimensional 
real projective space. The hypersurface is said to be hyperbolic in relation to 
a point if each real line passing through this point intersects the hypersurface 
in d real points. If these points are pairwise distinct, the surface is strictly 
hyperbolic. 


Example. An ellipse is strictly hyperbolic relative to its interior points, non- 
strictly hyperbolic relative to its boundary points, and not hyperbolic relative 
to its exterior points. 


The condition of hyperbolicity of a system means that the hypersurface 
of retardations in the projective space of contact elements at each point of 
space-time must be hyperbolic relative to the time point. 

Here the time point is the projectivization of the vector dt of the cotan- 
gent space to space-time. (The corresponding contact element is the tangent 
hyperplane of an isochrone.) 

A strictly hyperbolic surface of even degree 2n consists of n ovaloids dif- 
feomorphic to a sphere and situated one inside another from the one closest to 
the time point to the one most remote from it. Physically, these components 
correspond to the different “modes” or types of waves that can propagate in 
the given medium. For example, in an elastic medium there are longitudinal 
and transverse waves. Longitudinal and transverse waves propagating in the 
same direction have, in general, different velocities. 

The ovaloid closest to the time point corresponds to the slowest wave, 
the next to a faster wave, and so on, up to the one that arrives first, the 
fastest wave, which corresponds to the outermost ovaloid. This is clear from 
the order in which the points where the corresponding cones intersect the 
vertex O with the timelike world line g = c located alongside O. The point of 
intersection nearest to the isochrone corresponds to the smallest arrival time 
t of the perturbation from O at c, and consequently to the fastest wave. 

It should be kept in mind that two waves of a given type correspond to each 
direction; one propagates forward, the other backward. (The corresponding 
wave vectors are oppositely directed.) 


Example. For the wave equation 2n = 2 there is one mode, and to each 
direction in physical space there correspond exactly two waves, propagating 
in opposite directions. 


The basis of the short-wave asymptotics and, in particular the transition 
from the “physical” optics of the wave equation to the geometric optics of the 
eikonal equation is beyond the scope of these lectures. I remark only that a 
systematic implementation of this program leads to the appearance of inter- 
esting topological invariants in the asymptotic formulas, the so-called Maslov 
indices, which describe “the loss of one-fourth of the wave when a ray passes 
through a caustic” and which manifest themselves in quantum mechanics as 
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the correction of 1/2 in the Bohr-Sommerfeld quantum condition. The general 
statement of the quantization conditions thus leads to a topological object, 
the Maslov characteristic class on Lagrangian submanifolds of the symplectic 


phase space. 


Problem. Extend Weyl’s formula (the asymptotics of the number of eigen- 
functions with eigenvalues less than E) to the case of hyperbolic systems. 


Appendix A 


The Topological Content of Maxwell’s 
Theorem on the Multifield Representation 
of Spherical Functions 


In this appendix we present the proof of Maxwell’s theorem on the multi- 
field representation of spherical functions. In particular it is proved that the 
functions that admit such a representation form a vector space. This is the 
proposition that was left unproved above (Lecture 11, p. 113). 

Simultaneously we shall prove a number of interesting topological and 
algebraic corollaries of Maxwell’s theorem, which seem to have been proved 
first by Sylvester in a little-known note that also contains both the basic 
ideological principle of the treatise of Bourbaki and a warning against the 
dangers of abusing the formalization of mathematics. 

We recall the statement of the theorem (pp. 102 and 112). 


Theorem 1. The restriction to the sphere of an n‘*-order derivative of the 
function 1/r along n constant (translation-invariant) vector fields in RÌ coin- 
cides with a spherical function of degree n. Any nonzero spherical function of 
degree n can be obtained in this way using some set of n nonzero vector fields. 
For a given function these fields are uniquely determined (up to a nonzero 
constant factor and a permutation of the fields). 


The spherical functions of degree n form a vector space of dimension 2n+1. 

The set of functions represented by the multifield construction described 
in the theorem appears a priori to be essentially nonlinear. It follows from the 
theorem that the range of the corresponding multilinear mapping is a linear 
manifold. The uniqueness asserted in the theorem can be restated in purely 
topological terms. 


Theorem 2. The configuration space of n (possibly coincident) indistinguish- 
able points on the real projective plane (that is, the n? symmetric power 
Sym” (RP?)) is diffeomorphic to the real projective space of dimension 2n: 


Sym” (RP”) x RP?” . 
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The symmetric powers of nonorientable surfaces were computed (indepen- 
dently of Maxwell) by J. Dupont and G. Lusztig [4]. 


Remark. Theorem 2 is related to Viéte’s projective theorem 
Sym” (CP!) = CP” 
and is in a sense the quaternion version of that theorem. 


By regarding the Riemann sphere CP! as a two-sheeted covering of real 
projective space, we obtain as a corollary an algebraic mapping r : CP” — 
RP?” of multiplicity 2”, which is a higher-dimensional generalization of the 
classical theorem asserting that 


CP?/conj x S4. 


A.1. The Basic Spaces and Groups 


Consider the n-dimensional arithmetical quaternion space H” = @ Hy with 
its usual i-complex structure (i(ae + bi + cj + dk) = ai — be + ck — dj). Left 
multiplication by j operates on Hy, while preserving complex lines. It maps 
each line to its Hermitian-orthogonal line and acts on CP} = (Hj \ 0)/C* as 
an antiholomorphic involution op having no fixed points. 

Consider the Coxeter group B(n), which operates on the product (CP!)" 


by permuting factors and applying the mappings o, on certain factors. 


Theorem 3. The space of orbits of the action of B(n) on (CP!)” is 2n-di- 
mensional real projective space: 


(CP!)"/B(n) = RP? . 


According to Viete’s theorem, the space of orbits of the permutation group 
S(n) is complex projective space: 


(CP')"/S(n) = Sym"(CP!) = CP”. 


Thus we have obtained a natural mapping r : CP" — RP?” taking the orbit 
€ of the subgroup S(n) into the orbit r(é) of B(n) containing £. 

The group B(n) also contains the interesting subgroup Zo x S(n), con- 
sisting of permutations and permutations together with the antiholomorphic 
involutions op on each factor. The product of involutions Op OporatES on 
(CP!)"/S(n) as an involution o € Zo. 

The group embeddings S(n) — Zo x S(n) = B(n) generate mappings of 
the spaces of orbits 


(CP!)"/S(n) S (CP')"/(S(n) x Za) & (CP1)"/B(n) 


of multiplicities 2 and 2”~! respectively. 
The involution ø : (CP!)"/S(n) — (CP!)"/S(n) is the pre-image aWl(.), 
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Theorem 4. For even n the involution o operates on (CP')"/S(n) = CP” 
as the complex conjugation conj. 


Thus, for even n we obtain real algebraic mappings 
CP” & CP"/conj & RP?” 
of multiplicities 2 and 2”7! respectively. 


Remark. When n = 2, the second space is smooth (see, for example, [1] and 
[2]; this was probably known even before the publication of [2]): 


CP?/conj = S4. 
In this case the multiplicity of 8 is 2. 


The involution representing the two pre-images operates on S4 as the 
antipodal involution. 

I am grateful to S. Donaldson for this remark, which shows that the 
(strange) theorem of Maxwell is in a certain sense a multidimensional analog 
of the (no less strange) theorem that CP?/conj ~ S4. 

We shall show below that Maxwell’s theorem provides an explicit formula 
for a diffeomorphism CP?/conj — S4. 


Remark. In the majority of cases when we assert that two manifolds “coin- 
cide” we shall exhibit explicitly only a real algebraic homeomorphism between 
these manifolds. The pedantic verification that these homeomorphisms can be 
smoothed is in some cases left to the reader (see, however, Sect. A.4). 


A.2. Some Theorems of Real Algebraic Geometry 


Consider a real homogeneous polynomial f of degree n in the three variables 
(x,y,z). Maxwell’s theorem has the following strange algebraic corollary. 


Theorem 5. Any real homogeneous polynomial f of degree n has a unique 
representation as the sum of two such polynomials, one of which is the product 
of n linear real factors, and the other is divisible by x? + y? + z?. 


In particular, any real algebraic curve of degree n in the metric projective 
plane determines n real “principal axes” that are invariantly connected with 
the curve. 


PROOF. The real equation x?+y?+z? = 0 defines a real curve S in CP? having 
no real points. This curve, which is called an imaginary circle, is rational and 
is a topological sphere. Complex conjugation maps S into itself and operates 
on it as an antiholomorphic involution having no fixed points. (It is the usual 
antipodal involution of S?.) 


138 A. Topological Content of Maxwell’s Theorem 


The equation f = 0 defines a real algebraic curve K of degree n in CP?. 
The complex conjugation conj : CP? — CP? maps K into itself and permutes 
the 2n points of the intersection of K with S. 

Each pair of conjugate points of the intersection determines the line joining 
them. (The points are distinct, since conj has no fixed points on S?.) This line 
is real (since conj permutes two points belonging to it) and can be defined by 
an equation ax + by + cz = 0 with real coefficients. 

Consider the product g of n real linear functions obtained in this way. We 
shall show that f is proportional to g along S. 

The homogeneous polynomial f equals 0 at the 2n common points of the 
curves K and S. The curve S is rational. Let us choose a rational parameter 
t (for example, x = 2t, y = t? — 1, z = i(t? + 1)). Choosing coordinates in a 
suitable way, we can exclude the case when t = œ is one of the points of inter- 
section of S with K. The polynomials f (x(t), y(t), z(t)) and g(x(t), y(t), z(t)) 
of degree 2n have 2n common roots, and g # 0. Consequently, along S we 
have the equality f = cg, c = const. 

Thus, the homogeneous polynomial f — cg equals zero on S, that is, it has 
the form (x? + y? + z?)h(zx,y,z), where h is a real homogeneous polynomial 
of degree n — 2. The theorem is now proved. UO 


Remark. The uniqueness of the decomposition f = cg + (x? +y? + z7)h can 
be proved independently of the proof of existence. Assume that a second de- 
composition f = ¢'g' +(x? +y? +z?)h' exists. Then c’g’—cg = (x7 +h? +2?)h", 
where g and g’ are both products of n linear factors. Suppose g = 0 on the 
line l. On that line the polynomial c’g’ has n real roots and equals the right- 
hand side of the preceding formula, which has at most n — 2 real roots. Thus 
c’g’ = 0, and h” = 0 on l. Consequently, we can divide all three polynomials 
by the equation of l and prove uniqueness by induction on n. (Everything is 
obvious for n = 1, since cg = c’g’ by virtue of the fact that h” = 0.) 


Consider the projective space RP% of real algebraic curves of degree n 
(N = n(n + 3)/2). 

The real algebraic curves consisting of n real lines form a closed real alge- 
braic variety T of dimension 2n in this projective space (the image of (RP?)” 
under a multilinear mapping). The curves containing S form the projective 
subspace P = RP™, M = (n—2)(n+1)/2. In these terms Theorem 5 assumes 
the following form: 


Theorem 6. The varieties T and P in RPN are linked with linking coefficient 
1 in such a way that through any point of RP™ not belonging to the union of 
T and P there passes a unique real projective line joining T and P. This line 
intersects T (just as it intersects P) in a single point. 


It is now easy to prove Theorem 2. 
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PROOF OF ‘THEOREM 2. Consider the space H of dimension M +1 containing 
M-dimensional projective space P. Any such space intersects T. Indeed, let 
us choose a point O in H not belonging to P or T. The line joining O with T 
and P is contained in H and intersects T. 

The point of intersection is unique. Otherwise at a point of H there would 
be a point of O through which pass two lines joining T and P. 

Thus we have constructed a homeomorphism between T ~ Sym”(RP7?) 
and a manifold of spaces H of dimension M +1 containing P. This last mani- 
fold is RPN-M-1, which (at least on the topological level) proves Theorem 2, 
since N—-M—1=2n. O 


A.3. From Algebraic Geometry to Spherical Functions 


The derivatives of a harmonic function along constant vector fields are ob- 
viously harmonic functions. Thus all the repeated derivatives of the function 
1/r are harmonic except at 0. The following lemmas were proved in Lecture 11 


(p. 112). 


Lemma 1. The n“ derivative of the function 1/r along n constant vector 
fields has the form P/r?2"*1, where P is a homogeneous polynomial of degree n. 


Lemma 2. The homogeneous polynomial P of Lemma 1 is a harmonic func- 
tion. 


This lemma follows from the classical theorem on inversion, which we now 
recall. The harmonic function P/r?”+! is homogeneous and of degree —(n +1). 
On the unit sphere it coincides with P. It follows that the extension of this 
function from the unit sphere to the whole space as a homogeneous function 
of degree n is also harmonic. This extension is precisely P. 

To prove the theorem on inversion we consider the spherical Laplacian A, 
which is the Laplacian div grad on the unit sphere. We extend it as a homo- 
geneous function of degree k to R” \ {0} in such a way that any homogeneous 
function F of degree k is mapped by it into a homogeneous function AF of 
the same degree. 

According to Lecture 11 (p. 105) for any homogeneous function F of 
degree k we have 


AF =r°AF-—AF, where A= k? +k(m-—2). 
The following assertions are consequences of this formula. 


(i) A bounded harmonic homogeneous function of degree k on the unit 
sphere in R” is an eigenfunction of the spherical Laplacian correspond- 
ing to the eigenvalue — A. 

(ii) An eigenfunction of the spherical Laplacian corresponding to the eigen- 
value —A, when extended from the sphere to a homogeneous function 
of degree k on R™ \ {0}, becomes harmonic. 
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(iii) For any degree of homogeneity k there exists in R™ a dual degree 
k' = 2—m-—k such that a homogeneous harmonic function of degree k 
remains harmonic it we restrict it to the sphere and then extend it 
as a homogeneous function of the dual degree. For m = 3 the duality 
condition has the form k + k’ = —1. 


In particular, by Lemma 1, P/r?"*! has degree k = —(1 +n), m = 3, 
that is, the dual degree is k’ = n, from which Lemma 2 follows. 


Lemma 3. Any spherical function of degree n on S? can be represented in 
the form f(X,Y,Z)(1/r), where f is a homogeneous polynomial of degree n 
and X = ð/ðzx, Y = O/Oy, Z = 0/0z, and r? = xr? +y? + 2?. 


PROOF. The space of harmonic homogeneous polynomials of degree n is a 
vector space. It contains the subspace of harmonic polynomials whose restric- 
tions can be represented as described in Lemma 3 (by Lemmas 1 and 2). This 
space is obviously rotation invariant. But the representation of SO(3) in the 
space of spherical functions of degree n is irreducible. (Any function can be 
obtained from the Legendre polynomial of degree n by using rotations and 
forming linear combinations of rotated functions.) 

Thus, the subspace of Lemma 3 coincides with the whole space of spherical 
functions of degree n. O 


Lemma 4. Any spherical function of degree n can be represented in the form 
n 


fr(X,Y,Z)(1/r), where fr = [| (aX + BY +%Z) is the product of n real 
1=1 
linear factors. 
PROOF. According to Theorem 5 there exists a decomposition f(z, y,z) = 
fr(z,y,2z) + 9(x,y, z)(x? + y? + z?). Applying this decomposition to the rep- 
resentation of Lemma 3, we obtain 
F(X, Y,Z)(1/r) = fr(X,Y,Z)(1/r) +0, 


since X? + Y? + Z? = A and A(1/r) = 0. Thus, any spherical function has 
the multifield representation of Theorem 1. O 


Lemma 5. The multifield representation is unique. That is, the polynomial f; 
is uniquely determined by the spherical function. 


PROOF. Assume that fr and f+ are completely decomposable polynomials 
such that 
fr (X,Y, Z)(/r) = fr(X,Y,Z)(1/r) . 


According to Lemma | (applied n times), 
fr(X,Y, Z)(1/r) = (cfr(z,y,z)+r°g)r C+D, ce £0, 
f(X,Y,Z)(1/r) = (P fray, 2) +r?g) rD, of £0. 


Then fr(x,y, 2) — fr(2,y,z) = r?h(z,y, z). By Theorem 5, however, this can 
happen only if fr = f+. Theorem 1 is therefore proved. C 


A.4. Explicit Formulas 141 
A.4. Explicit Formulas 


Left quaternion multiplication by j maps the vector (z, w) = zj+we of C? = H! 
to (w, —Z). This vector is Hermitian-orthogonal to the original. 

We thus find an explicit formula t +» —1/t for an anti-holomorphic invo- 
lution of CP! having no fixed points. 

We can use a pair of Hermitian-orthogonal lines in C? to parametrize 
points in RP?. In that way we deduce Theorem 3 from Theorem 2. Further, 
we shall write out explicit formulas for the diffeomorphism 


Sym” (RP?) = RP?” 


of Theorem 2, using sets of n pairs of Hermitian-orthogonal lines in C? as 
coordinates on Sym” (RP?). 

We begin with the case n = 1. 

To each pair (z, w) and (w, —Z) of Hermitian-orthogonal vectors in C? we 
assign a quadratic form on the dual space. This form is a product of the linear 
functionals represented by these vectors: 


f = (zx + wy) (we — Zy) = for? + fixy + foy’ . 


Here (x,y) are coordinates on the dual plane of C?. Thus the coefficients of 
the quadratic form f have the form: 


foo=zw, fr=ww-2z, fo=—Zw. 
We remark that the coefficient fı is real-valued, while fo = — fo. We shall 
regard fo and fı as coordinates on the space R3 of forms f. 
Multiplying the original vector (z, w) by a complex number leads to multi- 


plying the resulting vector by the square of the absolute value of this number. 
Thus we obtain the mapping 


F:CP! — S? = (R? \ {0})/Rt , 


which maps the complex lines from C? to real rays in R3. 

The mapping (fo, f1) maps the point (z,w) on the complex line to a com- 
pletely definite point of the ray. For example, choosing z = t, w = 1, we obtain 
fo=t, fi =1-—tt. 

If the original point on the line is normalized by the condition |z|? + |w]? = 
1, its image lies on the ellipsoid |2 f9|? + |f1|? = 1 (which we can also call a 
sphere, regarding 2fo and fı as coordinates). 

Thus the Riemann sphere CP! = S3/S1 maps to the unit sphere in R? with 
coordinates 2fo and fı by the diffeomorphism (2fo, f1). Regarding t = z/w as 
a coordinate on CP", we obtain a mapping of the t-plane to the unit sphere 
in RÌ, which is precisely the stereographic projection. The formulas that we 
write out below are in this sense a multidimensional generalization of the 
stereographic projection. 
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By replacing the original line in CP? by its Hermitian-orthogonal line, we 
reverse the sign of the mapping F. Indeed, replacing z by w and w by —z 
changes the sign before both fo and fı. Thus, F transforms the involution 
j : CP! — CP! (which maps every line to its Hermitian complement) into an 
antipodal involution of the sphere S? in R. 

We now apply the analogous construction to the n‘” symmetric power of 
the space RP”. We begin with the nt” symmetric power of CP!. By definition 
a point of the complex manifold 


th 


Sym” (CP!) = CP” 


is a (non-ordered) set consisting of n lines in C?. 
Choosing n nonzero vectors (z%, wz) and multiplying the corresponding 
linear forms on the dual plane, we obtain a binary n-form 


H(z,y) = [| (eet + wey) = hor” +: + hny” . 
k=1 


The coefficients in this form are (homogeneous) coordinates on CP” = 
Sym"(CP!) (which define a smooth holomorphic structure on this space). 

If wk Æ 0, one can set wk = 1. As affine coordinates we obtain the basic 
symmetric functions of the variables {zx} 


ho = on(z), ES hn-1 = 01(2), (hna = 1). 


In what follows we shall use these a, as local coordinates on Sym” (RP?). 

We begin with n pairs of Hermitian-orthogonal lines in C?. We choose 
a representative of each pair in such a way that neither of the lines chosen 
coincides with any of those not chosen (at the given point, and consequently 
also in some neighborhood of it, where our coordinate system is valid). In 
order to assure that the nondegeneracy condition above holds in repeating 
the pair, it suffices to choose the same line each time. 

We denote by (zx, 1) the vectors that define the chosen lines (k = 1,..., n). 
As local coordinates on the real manifold Sym”(RP?) we shall use (the real 
and imaginary parts of) the n complex numbers 


O12 issn a eee Ones ce eens 


The orthogonal lines are defined by the vectors (1, —Z,). We construct the 
symmetrized 2n-form 


n n 
f= | [Cre +y) [[(2- zry) = for?” +-+ fany” 
k=1 k=1 


We shall see below that the coefficients fp are polynomials in o and a. 
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Theorem 7. The mapping F : C” — R?"*! taking the point (o1,...,0n) to 
(fo,---,fon) defines (locally) a diffeomorphism of the manifold Sym” (RP?) 
into the space RP?” of rays in R?"+!. In coordinates the mapping F is given 
by the following explicit formulas: 


fo = On , 
fi = On-1 — On ) 
fo = On-2 — On-101 + On_-202 , 


..esesssoessssosoosssossesssoseooo‘wolno‘ols‘l‘loi‘’ 


PROOF. It is not difficult to see that 


n 


(zk£ +y) = On” +0n-1£° y+ +y”, 
1 


Eea 
il 


(x2 — Zy) = 2" — Gyn" !y + Gon” Ay? +--+ (-1)"G yy” . 


ms 


> 
il 
_ 


Multiplying these two polynomials, we obtain (according to F. Aicardi) the 
formulas given above for the coefficients of the product. 

We also find that fon—~% = (—1)"~* fx. In particular, the middle coefficient 
fn is real-valued. 

It remains to be verified that the Jacobian is nonzero in this region. This 
can be proved without computation. The Jacobian we are interested in has 
order 2n + 1. One of the columns is the vector 


D2 (fos Foes see T a 
The other 2n columns are its derivatives 
(O0@/001, 0©/001,...,06/O0n, O®B/OGp) . 


We represent this nonholomorphic Jacobian as the value at the point r = o of 
the nonholomorphic Jacobian 7 (o,7) defined by the following construction. 
Consider the product 


n n 
[[ (er +y) [| (@ - wey) = For?” +--+ + Fony” . 
k=1 k=1 


The coefficients Fk are polynomials in o1(z), ..., On(z) and 7 =0,(w), ..., 
Tn = On(w). We denote by W = (Fo,..., Fən) a vector-valued function of the 
variables o and 7, and we consider the determinant of the matrix 


(Y,0W/00),0W/O7,,...,U/Ion,0V/OT_) . 
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This determinant T(o,7) is nonzero at the point o(z), r = o(zZ). Indeed, 
it follows from our non-coincidence condition that no zę coincides with any 
—1/z,. Hence, (o1,...,0n) and (71,...,7) form a local holomorphic coor- 
dinate system on CP?” = Sym?” (CP!) in some neighborhood of the point 
in question o, T = G. That is, the determinant T(o,7) is nonzero at that 
point. But the determinant we wish to prove nonzero coincides with the value 
T(o,c@), since fon-k = Ty Fe. Thus, we have proved that the mapping 
Sym”(RP?) — RP?” is a local diffeomorphism. 

From Sect. A.2 we know that this mapping is a homeomorphism. We have 
now proved Theorem 2 completely. O 


A.5. Maxwell’s Theorem and CP?/conj ~ S* 


The explicit formulas in Sect. A.4 also make it possible to construct a diffeo- 
morphism of the sphere $°” into the orbit space described below. 

We begin with the complex manifold (CP')” of ordered sets of n lines 
in C?. On this manifold we consider the following smooth (nonholomorphic) 
action of the Coxeter group D(n). An element of D(n) operates as a permu- 
tation of factors together with the replacement of an even number of lines by 
their Hermitian-orthogonal complements. 


Theorem 8. (CP!)"/D(n) = S?", and the corresponding diffeomorphism is 
locally given by the formulas of Theorem 7. 


PROOF. Permutations do not change the binary 2n-form f. Replacing a line 
by its complement changes the sign of the corresponding factor (z,x + y) x 
(x —Z,y). Consequently, an even number of such replacements does not change 
the sign of f (while an odd number does change it). O 


It is interesting that the relation (CP')"/D(n) ~ S?” thereby proved 
can be regarded as an informal extension of a theorem of Chevalley: The 
orbit space of the action of a real (2n —1)-dimensional group on C2” (which 
must be regarded as a generalization of the Coxeter group) is the smooth real 
space R2*!, 


Example. When n = 2, we obtain 
(CP')?/D(2) = S84, 


where the group D(2), which consists of 4 elements, operates on pairs of lines 
from C? as permutations and (possibly) replacement of both lines by their 
Hermitian complements. 

But (CP!)?/5(2) = Sym?(CP!) = CP?. Consequently, (CP!)?/D(2) = 
CP? /(j), where (j) is the replacement of both lines by their complements. 

Consider the complex manifold Sym” (CP!) = CP” of unordered sets of 
n lines in C?. The operation j of replacing each line by a line Hermitian- 
orthogonal to it, which acts on CP”, is an (antiholomorphic) involution. 
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Theorem 9. For even n the involution 7 : CP” — CP” coincides with com- 
plex conjugation (in some coordinates). 


PROOF. The coefficients of the form 
Tr n 
Haans [erk S a g" 
k=1 k=0 


are natural (homogeneous) coordinates on CP” = Sym”(CP!). Let us de- 
scribe the action of 7 on the coefficients hk. 
Let us transfer the action of 7 to the dual space: 


Det — Zey = zk (~Y) + we(Z) . 
Thus, the transferred form is given by the formula 
Ha —-:(z,y) = H, w(—9, T) . 


In terms of the coefficients we obtain the following expression for the trans- 
formed form: 


> ha(—gyr-# (a) = DO(-1)" Fay — 
k=0 k=0 
Thus, the action of 7 on the coefficients of the form H is given by the 
formula (jh) = (—1)fhn-p. Since n is even, we also obtain (jh)n-k = 
(—1)*hy. The required coordinates are hy + hn—, and i(hg — hn-p) for even 
k, and i(hk + hn-k) and hk — hn-k for odd k. (Naturally, we never consider 
hn/2 = hn/2-) 2 


When n = 2, our results reduce to an explicit formula for the classical 
diffeomorphism CP? /conj ~ S4. “Maxwell’s theorem” 


(CP')"/D(n) = Ss?” 


extends this diffeomorphism to higher dimensions. 


A.6. The History of Maxwell’s Theorem 


Maxwell’s own version of this theorem can be found in his major book Elec- 
tricity and Magnetism, Vol. 1, Chap. IX, Sects. 129-133 (pp. 222-233 in [5]). 
Sylvester criticized his reasoning, saying in [7]: 


I am rather amazed that such a distinguished author has given no 
attention to the fact that there always exists one and only one real 
system of poles corresponding to a given harmonic... 
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With all due respect to the great talent of Professor Maxwell, I must 
confess that the derivation of the purely analytic properties of spheri- 
cal harmonics from “Green’s theorem” and the “principle of potential 
energy,” ... as done by him, seems to me an incorrect method, of the 
same type and just as convincing... as if the rule for extracting the 
square root were deduced from Archimedes’ law of floating bodies. 


Sylvester proposed his own approach, which seems to be equivalent to 
Theorem 5: 


The method of poles for representing spherical harmonics invented 
or developed by Professor Maxwell actually reduces to choosing a 
suitable canonical form for a ternary quantic, the sum of the squares 
of whose variable (here, differential operators) equals zero. 


Nevertheless, Sylvester was not concerned with developing this question 
in detail (“being very pressed for time and in haste to catch a boat sailing for 
Baltimore in 24 hours”). The details of the proof appeared in [6] and, later, 
in [3]. 

Sylvester mentioned the connection of his theory of integrals of products 
of spherical harmonics with Ivory’s theorem on the attraction of an ellipsoid 
and proposed several generalizations of these ideas. 

It appears that neither the algebraic nor the philosophical ideas in this 
Note of Sylvester’s were understood or developed by the mathematical com- 
munity. The pages containing the note were not included in his Collected 
Papers in the library of the Paris Ecole Normale Supérieure. 

The note contains the following (anti-Bourbaki) paragraph: 


It is by no means unusual for a mathematical investigation... that 
a part is in some sense larger than the whole — the reason for this 
striking intellectual phenomenon is that, in regard to mathematics all 
quantities and connections must be regarded (as experience teaches) 
as being in a state of continuous change like the flow of a stream. 


This general philosophy led him to the conclusion that “... a general 
proposition should be easier to prove than any special case of it.” 

From the fact that Sylvester arrived at this last conclusion we can deduce 
that the truth of this important principle of him (borrowed by Bourbaki nearly 
a century later) does not entail the necessity of a lamentable petrification of 
the “flux” of mathematics. ' 
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Appendix B 


Problems 


B.1. Material from the Seminars 


Problem 1. Given a vector field v with a singular point 
of focus type (see the figure), find all solutions of the 


L 
equation L,u = 0. ay . \ 
Problem 2. The vector field v is given by A 


T=y, 
ST 
An initial function for the Cauchy problem is prescribed on the line z = 1. 


a) What is a necessary condition on the initial function al for the 
existence of a solution of the Cauchy problem? 


b) Is the solution unique? 
Problem 3. Given the equation 
Ut + Uur =0, (*) 
a) find its characteristics. 
b) Show that Newton’s equation 


A2y = 

Ore 
reduces to (*) under the hypothesis that y(t) is the position of a particle on 
the line and u(t, x) is the velocity of the particle that is at the point x on the 
line at time t. 
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Problem 4. a) Investigate the existence and uniqueness of the solution of the 


equation 
2 2 
Ox Oy 


with the initial condition Ps = 0 in the regions y > x? and y < 2”. 


b) Find the curve (see the figure) bounding 
the region in which the solution is unique. Prove 
that this curve is the set of centers of curvature 
of the parabola y = 2”. 


c) Do the same, but with the initial condition 


IERE. = 
[\ e+ =l 
d) Consider a small C%-perturbation of the 
ellipse. The perturbed curve y is shown in the 
figure, and l1, l2, l3, and l4 are the lengths of the 
y X arcs of the enveloping system of normals to y. 
Prove that lı + l3 = lg + l4. 


Problem 5. Does there exist a unique solution of the Cauchy problem 


ðu 20 
r(e +y2)— Hy = 


Ox oo «| 


= = 


in a neighborhood of the point (zo, 0) on the x-axis? 


Homework 


Problem 6. Let a = du — pı dz, — p2 dz2. Prove that no three-dimensional 
integral surface exists. 


Problem 7. Find the characteristics in the (2n + 1)-dimensional space of 
l-jets of a linear homogeneous equation, a linear nonhomogeneous equation, 
and a quasi-linear equation. 


Problem 8. Investigate the existence and uniqueness of the solution of the 
problem yur = Tuy, 7 = cosy in a neighborhood of the point (1, yo)- 


Problem 9. Find the maximal t for which the solution of the equation 


Ou Ou : 
uz—-=sing, ul|,_.=0 


att Ox 


can be extended over [0, t). 
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The Wave Equation, Waves, the Korteweg-de Vries Equation 


Problem 10. Find the characteristics of the wave equation 


Pu u 

— =c’—~, TER. 

ot? Ox? 
Problem 11. What form does the wave equation assume if the characteristics 
are used as coordinates? 


Problem 12. Find the general form of the solution of the wave equation. 


Problem 13. The initial condition uba = u(x) is prescribed for the wave 
equation. Find the solution with a given initial condition. Is it unique? 


Problem 14. Transform the Korteweg-de Vries (KdV) equation u; = 6uu, — 
Ure into an ordinary second-order differential equation using the substitution 
u = p(x — ct) (we are seeking a solution in the form of a moving wave). 


Homework 
p; O 
Problem 15. Consider the change of coordinates z’ = a yj. Express = in 
u u 8u = 
terms of — and ———— in terms of ————., and conversely. 
Oy" Ox*OxI Oy* Oy) 


Problem 16. Draw the phase portrait of the equation 

~ = 3y?+Cyp+K 

(the KdV equation after the change of variables u = y(x — ct)). 

Problem 17. On the phase plane (see Problem 16) find the solution in the 
form of a wave for which t > +œ, y - 0. 


The Wave Equation 


Problem 18. Prove d’Alembert’s formula 


z+at 
u(t,2) = 5(ole at) +y(e+at)) +5 f vay, 


which gives the solution of the equation of the vibrating string Un = a?uzz 
with initial conditions u(z,0) = y(x), uz(z,0) = y(x). 


The following problems (Nos. 19-24) involve a semi-infinite string (x > 0) 
with a free left-hand endpoint (u,(0,t) = 0) or a fixed left-hand endpoint 
(u(0, t) = 0). 


152 B. Problems 


Problem 19. Let the left-hand endpoint of the string be fixed, and suppose 
initial conditions y(x) and y(x) are given for x > 0. How should the functions 
p and w be extended to the set x < 0 in order that a bounded solution of 
the resulting Cauchy problem will coincide with the solution of the original 
Cauchy problem in the region {(z,t) : x > 0}? 


Problem 20. Answer the same question for the string with a free left-hand 
endpoint. 


Problems 21—24. Draw movies, that is, representations of the solutions 
of the Cauchy problem for different values of t > 0, for the given initial 
conditions. 


21. (z): , w(x) = 0, and the endpoint x = 0 is fixed. 


22. Same data, but the endpoint x = 0 is free. 


23. y(x) = 0, w(z): , and the endpoint x = 0 is fixed. 


24. Same data, but the endpoint x = 0 is free. 
Problem 25. Draw movies for the initial conditions as in problems 21-24, 
but for a string of finite length (0 < x < 1), for the following cases: 

a) both endpoints free, 

b) both endpoints fixed, 


c) one endpoint fixed and the other free. 


Homework 


Problem 26. Find the general solution of the Cauchy problem for a string 
of finite length with boundary conditions Als = f(t), ee = 0, if 

a) ulio = Melo 

b) ul = (x), ut| 9 = y(x). 
Problem 27. Suppose that under the conditions of the preceding problem 
f(t) is periodic with period T. 

a) Is the solution u(x,t) a periodic function with respect to t for arbitrary 
y and w? 

b) Verify that there is no periodic solution when T = (p/q)T, where p and q 
are integers and 7 = 2l/a. 
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Test 


Version 1 

1. Solve the Cauchy problem cuz + uy = 0, ul =o = sing. 

2. For which values of t can the solution of the Cauchy problem u; + uuz = 
—23, aN = 0 be extended to the entire half-open interval [0, t)? 


3. Draw movies for the following problems: 
i 
a) Ut = Urr, Y = 0, Y: "E , the endpoint xz = 0 is fixed, the 
1/2 


endpoint x = 1 is free. 


b) Utt = Uzr, ahan SW ps , urlo =0, ul = 0. 


4, Find the solution of the Cauchy Problem ut = uzz with initial condition 
ulia = sin? T; us|, —o = sing. 


1 1 1 
5. Is the series 1 + = cos 2x + = cos 3z +---+—cosnx +--- the Fourier 
n 
series of a continuously differentiable function? 
6. Let {yx} be an orthonormal system. The series $` (f, px )yrR converges 
to f in Lə. Express ||f||z, in terms of the Fourier coefficients (Parseval’s 
equality). 


Version 2 
, ðu ðu 3 i 
1. Are all solutions of the equation o = By on R^ functions of ry? 
z 
2. Find the solution of the equation yuu; +xuu, = ry whose graph passes 
through the curve z = y? +u? = 1. 


3. Draw movies for the following problems. 


b) Utt = Urz, P TA w = 0, the endpoint x = 0 is fixed, and 


the endpoint x = 1 is free. 


4. Solve the Cauchy problem uz; = Uzr with the initial conditions ül = 


COs T, thii = cos? z. 


5. Is the series 1 + cos2r +--+ cosnz +--- the Fourier series 


n 
of a continuously differentiable function? 


ninn 


6. Let y € C%(S'); more precisely, the function y is holomorphic in 
|Imz| < @ and p(z + 27) = (z). Prove that lo] < cev*, 
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Conservative Systems. Lissajous Figures 


Problem 28. Consider the equation ¢ = —w?z. 


a) Find all solutions. 
b) Suppose the solution has the form Asinwt + Bcoswt = Csin(wt + y). 
Prove that A? + B? = C?. 


In what follows we are considering the equation 


# = —VU(z). (B.1) 


(x? + x3), find all solutions 


N| = 


Problem 29. Given the potential U (z1, £2) = 
of the equation (B.1) in the form shown in part b) of Problem 28. 
Problem 30. Given the potential U (£1, £2) = £? + 4£ız2 + 423, find all 


the solutions of the equation (B.1). How do the solutions behave in the 
(£1, £2)-plane? 


Problem 31. Prove that the total energy of the conservative system (B.1) is 
a first integral of this system. 


In what follows we consider the potential U(x) = > + 5 2. With a 


suitable choice of origin for the time t the general solution of (B.1) can be 
written as 


Tı = A sint, 


T2 = Agsin(wt + y). 


Problem 32. Let w = 1, y = 5, 7, and 0. How will the Lissajous figures 
look in the rectangle |zı| < Aj, |v2| < A2? What will happen if y varies from 
O to r? Find the coefficients of the similarity transformation of the ellipses 
U(z1, z2) = E and the Lissajous figure when y = = 


Homework 


Problem 33. Given the potential U (x,y) = 5x? +5y? — 7zy. Solve the equa- 
tion (B.1). 


2 

. , 1 
Problem 34. Given the potential U (z1, £2) = + + Ti we. Consider the 
E Ge. ies 


ellipse U (x1, £2) < E. Let Ey = +i +—, Fy = 2 and E = Ei + Eo. 


Then zı and z> lie in the strips oF p |z| < V2E2. Prove that the 
rectangle |x] < V261 = Aj, lzəl| < V2E> = Ap is inscribed in the ellipse 
U(z1, £2) < E. 
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Problem 35. Prove that if 
a) w = %7, the Lissajous figure is a closed curve, 


b) w ¢ Q, it is not closed. 


Problem 36. Prove that if w = n, there exists a phase y such that the 
Lissajous figure is the graph of a polynomial of degree n (the Chebyshev 
polynomial p(x) = cos(n arccos x)). 


2 3% 
Problem 37. Prove that the ellipse n + x = 1 and the arc of the sinusoidal 


curve y = 3cos(xz/4), 0 < x < 87, have the same length. 


Harmonic Functions 


Problem 38. Prove that the angle subtended by a line segment from a point 
in its plane is a harmonic function on the covering of the plane with the 
endpoints of the interval removed. 


Problem 39. Construct a function f harmonic in the disk and assuming 
given values Cı and C% on the arcs Sı and S2, where Sı U Sg = S!. 


Problem 40. Solve the analogous problem for a partition of the circle into n 
arcs S1,..., Sn with values C),...,Ch. 


Problem 41. Write Newton’s equation of motion for a free particle in polar 


coordinates. 


Homework 


Problem 42. Prove that the Lagrangians 


) eg te : +4, } 
0,7 a7 


minimize the action on a fixed curve (the shortest curve joining two points). 
Problem 43. Compute A in polar and spherical coordinates. 
Problem 44. Prove that the magnitude of the solid angle subtended from a 


variable point of space by a fixed closed contour in R® is a harmonic function 
of the variable point on the manifold covering the complement of the curve. 
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B.2. Written Examination Problems 


1995' 
1 (1). Is the surface 


2 
2 
z=rt+e * 
; () 
in Euclidean space smooth? (Here r? = x? + y? + 2°.) Is it convex? Find its 
curvature. 


2 (2). Find the value at the origin of a double-layer potential of density 1 
distributed on the surface whose equation is (*). 


3 (1,2,3,3,6). Compute the average values over the surface (x) for the 
following functions: a) z; b) 1/r; c) z/r?; d) r?; e) 1/r°. 


4 (5). Solve the interior Dirichlet problem with boundary condition u = 
1/r? on the surface (x) for Laplace’s equation Au = 0. 


5 (5). Consider the single-layer potential of density z distributed on the 
surface (*). Find the average value of this potential over the surface of the 
sphere r? = 1. 


6 (2,4). Find the upper and lower bounds of the Dirichlet integral 


fJ [tera dz dy dz 


over the region bounded by the surface (*) on the set of smooth functions in 
the closure of this region taking the value r? on its boundary. 


7 (6). Find the value of the solution f of the equation 


0 
2 + div(f gradu) = f? 
ot 
for small |t| at the point (x = y = 0, z = 1/2) with initial condition (f = 1 
when t = 0), where u is the solution of Problem 4 in this examination. 


HINT. It is not necessary to know the answer to Problem 4 in order to do 
this. Problem 7 can be solved independently. 


! The number of points awarded for a correct solution of a problem or part of a 
problem is indicated in parentheses next to the problem number. The recom- 
mended criteria for grading in a three-hour examination are: average — 12; good — 
16; excellent - 26 (the maximum possible is 40). 
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1996 


1. Given the vector field v(x, y) = yz + re and the function ul = 
f(y), under what condition on f in a neighborhood of the point (1,0) does 
there exist a solution of the Cauchy problem for the equation Lyu = 0? Is it 
unique? 


2. Consider the equation 
£=-VU(z), (xx) 
where z = (£1, £2) and U(x) = 2? + x3 + ax, 29. For which (real) values of a 
are all the solutions of the equation (**) periodic? 


3. Electric charge is distributed on two lines in R°: on the line z = 1, 
y = x with density 1 and the line z = —1, y = —z with density —1. Find the 
equipotential surfaces of the field created by these charges. 


4. The function v is defined on the sphere S? : r? + y? + 2? = 1 and is 
harmonic except at the point N = (0,0,1). Let R? be the plane z = 0 and 
p : S? \ N — R? the stereographic projection. Let u(x, y) = u(p7'(z,y)). 


20 
Prove that f u,(1,y) dy = 0. 
0 


Universitext 


Aksoy, A.; Khamsi, M. A.: Methods in Fixed 
Point Theory 


Alevras, D.; Padberg M. W.: Linear Opti- 
mization and Extensions 


Andersson, M.: Topics in Complex Analysis 
Aoki, M.: State Space Modeling of Time Se- 
ries 

Arnold, V.I.: Lectures on Partial Differential 
Equations 

Audin, M.: Geometry 

Aupetit, B.: A Primer on Spectral Theory 


Bachem, A.; Kern, W.: Linear Programming 
Duality 

Bachmann, G.; Narici, L.; Beckenstein, E.: 
Fourier and Wavelet Analysis 

Badescu, L.: Algebraic Surfaces 


Balakrishnan, R.; Ranganathan, K.: A Text- 
book of Graph Theory 


Balser, W.: Formal Power Series and Linear 
Systems of Meromorphic Ordinary Differen- 
tial Equations 

Bapat, R.B.: Linear Algebra and Linear Mod- 
els 


Benedetti, R.; Petronio, C.: Lectures on Hy- 
perbolic Geometry 

Berberian, S. K.: Fundamentals of Real Anal- 
ysis 

Berger, M.: Geometry I, and II 

Bliedtner, J.; Hansen, W.: Potential Theory 
Blowey, J.F; Coleman, J.P; Craig, A. W. 


(Eds.): Theory and Numerics of Differential 
Equations 


Borger, E.; Grädel, E.; Gurevich, Y.: The Clas- 
sical Decision Problem 


Böttcher, A; Silbermann, B.: Introduction to 
Large Truncated Toeplitz Matrices 


Boltyanski, V; Martini, H.; Soltan, P. S.: Ex- 
cursions into Combinatorial Geometry 


Boltyanskii, V.G.; Efremovich, V.A.: Intu- 
itive Combinatorial Topology 


Booss, B.; Bleecker, D.D.: Topology and 
Analysis 

Borkar, V. S.: Probability Theory 

Carleson, L.; Gamelin, T. W.: Complex Dy- 
namics 

Cecil, T. E.: Lie Sphere Geometry: With Ap- 
plications of Submanifolds 

Chae, S. B.: Lebesgue Integration 
Chandrasekharan, K.: 
Transform 

Charlap, L. S.: Bieberbach Groups and Flat 
Manifolds 

Chern, S.: Complex Manifolds without Po- 
tential Theory 

Chorin, A.J.; Marsden, J.E.: Mathematical 
Introduction to Fluid Mechanics 


Classical Fourier 


Cohn, H.: A Classical Invitation to Algebraic 
Numbers and Class Fields 

Curtis, M. L.: Abstract Linear Algebra 
Curtis, M. L.: Matrix Groups 

Cyganowski, S.; Kloeden, P; Ombach, J.: 
From Elementary Probability to Stochastic 
Differential Equations with MAPLE 

Dalen, D. van: Logic and Structure 

Das, A.: The Special Theory of Relativity: A 
Mathematical Exposition 

Debarre, O.: Higher-Dimensional Algebraic 
Geometry 

Deitmar, A.: A First Course in Harmonic 
Analysis 

Demazure, M.: Bifurcations and Catastro- 
phes 

Devlin, K.J.: Fundamentals of Contempo- 
rary Set Theory 
DiBenedetto, E.: 
Equations 


Diener, F; Diener, M.(Eds.): Nonstandard 
Analysis in Practice 


Degenerate Parabolic 


Dimca, A.: Singularities and Topology of 
Hypersurfaces 

DoCarmo, M. P: Differential Forms and Ap- 
plications 

Duistermaat, J. J.; Kolk, J. A. C.: Lie Groups 


Edwards, R.E.: A Formal Background to 
Higher Mathematics Ia, and Ib 

Edwards, R.E.: A Formal Background to 
Higher Mathematics IIa, and IIb 

Emery, M.: Stochastic Calculus in Manifolds 
Endler, O.: Valuation Theory 

Erez, B.: Galois Modules in Arithmetic 
Everest, G.; Ward, T.: Heights of Polynomials 
and Entropy in Algebraic Dynamics 
Farenick, D. R.: Algebras of Linear Transfor- 
mations 

Foulds, L. R.: Graph Theory Applications 
Frauenthal, J. C.: Mathematical Modeling in 
Epidemiology 

Friedman, R.: Algebraic Surfaces and Holo- 
morphic Vector Bundles 

Fuks, D. B.; Rokhlin, V. A.: Beginner’s Course 
in Topology 

Fuhrmann, P. A.: A Polynomial Approach to 
Linear Algebra 

Gallot, S.; Hulin, D.; Lafontaine, J.: Rieman- 
nian Geometry 

Gardiner, C. F: A First Course in Group The- 
ory 

Garding, L.; Tambour, T.: Algebra for Com- 
puter Science 

Godbillon, C.: Dynamical Systems on Sur- 
faces 

Goldblatt, R.: Orthogonality and Spacetime 
Geometry 

Gouvêa, F. Q.: p-Adic Numbers 

Gustafson, K. E.; Rao, D. K.M.: Numerical 


Range. The Field of Values of Linear Opera- 
tors and Matrices 


Gustafson, S.J.; Sigal, I.M.: Mathematical 
Concepts of Quantum Mechanics 

Hahn, A. J.: Quadratic Algebras, Clifford Al- 
gebras, and Arithmetic Witt Groups 

Hajek, P; Havránek, T: Mechanizing Hy- 
pothesis Formation 

Heinonen, J.: Lectures on Analysis on Metric 
Spaces 


Hlawka, E.; Schoifengeier, J.; Taschner, R.: 
Geometric and Analytic Number Theory 


Holmgren, R. A.: A First Course in Discrete 
Dynamical Systems 


Howe, R., Tan, E.Ch.: Non-Abelian Har- 
monic Analysis 

Howes, N.R.: Modern Analysis and Topol- 
ogy 

Hsieh, P-F; Sibuya, Y. (Eds.): Basic Theory 
of Ordinary Differential Equations 

Humi, M., Miller, W.: Second Course in Or- 
dinary Differential Equations for Scientists 
and Engineers 


Hurwitz, A.; Kritikos, N.: Lectures on Num- 
ber Theory 


Iversen, B.: Cohomology of Sheaves 

Jacod, J.; Protter, P: Probability Essentials 
Jennings, G. A.: Modern Geometry with Ap- 
plications 

Jones, A.; Morris, S. A.; Pearson, K.R.: Ab- 
stract Algebra and Famous Inpossibilities 
Jost, J.: Compact Riemann Surfaces 

Jost, J.: Postmodern Analysis 

Jost, J.: Riemannian Geometry and Geomet- 
ric Analysis 

Kac, V.; Cheung, P: Quantum Calculus 


Kannan, R.; Krueger, C. K.: Advanced Anal- 
ysis on the Real Line 


Kelly, P; Matthews, G.: The Non-Euclidean 
Hyperbolic Plane 


Kempf, G.: Complex Abelian Varieties and 
Theta Functions 


Kitchens, B. P: Symbolic Dynamics 


Kloeden, P; Ombach, J.; Cyganowski, S.: 
From Elementary Probability to Stochastic 
Differential Equations with MAPLE 


Kloeden, P E.; Platen; E.; Schurz, H.: Nu- 
merical Solution of SDE Through Computer 
Experiments 

Kostrikin, A. I.: Introduction to Algebra 
Krasnoselskii, M. A.; Pokrovskii, A. V.: Sys- 
tems with Hysteresis 


Luecking, D. H., Rubel, L. A.: Complex Anal- 
ysis. A Functional Analysis Approach 


Ma, Zhi-Ming; Roeckner, M.: Introduction 
to the Theory of (non-symmetric) Dirichlet 
Forms 


Mac Lane, S.; Moerdijk, I.: Sheaves in Geom- 
etry and Logic 

Marcus, D. A.: Number Fields 

Martinez, A.: An Introduction to Semiclas- 
sical and Microlocal Analysis 

Matoušek, J.: Using the Borsuk-Ulam Theo- 
rem 

Matsuki, K.: Introduction to the Mori Pro- 
gram 

Mc Carthy, P.J.: Introduction to Arithmeti- 
cal Functions 

Meyer, R. M.: Essential Mathematics for Ap- 
plied Field 

Meyer-Nieberg, P: Banach Lattices 
Mikosch, T.: Non-Life Insurance Mathemat- 
ics 

Mines, R.; Richman, F; Ruitenburg, W.: A 
Course in Constructive Algebra 

Moise, E. E.: Introductory Problem Courses 
in Analysis and Topology 
Montesinos-Amilibia, J. M.: Classical Tessel- 
lations and Three Manifolds 

Morris, P.: Introduction to Game Theory 
Nikulin, V. V.; Shafarevich, I. R.: Geometries 
and Groups 

Oden, J. J.; Reddy, J. N.: Variational Methods 
in Theoretical Mechanics 

Oksendal, B.: Stochastic Differential Equa- 
tions 

Poizat, B.: A Course in Model Theory 
Polster, B.: A Geometrical Picture Book 
Porter, J.R.; Woods, R.G.: Extensions and 
Absolutes of Hausdorff Spaces 

Radjavi, H.; Rosenthal, P.: Simultaneous Tri- 
angularization 


Ramsay, A.; Richtmeyer, R. D.: Introduction 
to Hyperbolic Geometry 


Rees, E. G.: Notes on Geometry 


Reisel, R. B.: Elementary Theory of Metric 
Spaces 


Rey, W.J.J.: Introduction to Robust and 
Quasi-Robust Statistical Methods 


Ribenboim, P: Classical Theory of Algebraic 
Numbers 


Rickart, C. E.: Natural Function Algebras 
Rotman, J. J.: Galois Theory 


Rubel, L. A.: Entire and Meromorphic Func- 
tions 


Rybakowski, K. P: The Homotopy Index and 
Partial Differential Equations 


Sagan, H.: Space-Filling Curves 

Samelson, H.: Notes on Lie Algebras 

Schiff, J. L.: Normal Families 

Sengupta, J.K.: Optimal Decisions under 
Uncertainty 

Séroul, R.: Programming for Mathemati- 
cians 

Seydel, R.: Tools for Computational Finance 
Shafarevich, I. R.: Discourses on Algebra 
Shapiro, J. H.: Composition Operators and 
Classical Function Theory 

Simonnet, M.: Measures and Probabilities 
Smith, K.E.; Kahanpää, L.; Kekäläinen, P; 
Traves, W.: An Invitation to Algebraic Ge- 
ometry 

Smith, K. T.: Power Series from a Computa- 
tional Point of View 

Smorynski, C.: Logical Number Theory I. An 
Introduction 

Stichtenoth, H.: Algebraic Function Fields 
and Codes 

Stillwell, J.: Geometry of Surfaces 

Stroock, D. W.: An Introduction to the The- 
ory of Large Deviations 

Sunder, V. S.: An Invitation to von Neumann 
Algebras 

Tamme, G.: Introduction to Etale Cohomol- 
ogy 


Tondeur, P: Foliations on Riemannian Man- 
ifolds 


Verhulst, F: Nonlinear Differential Equa- 
tions and Dynamical Systems 


Wong, M. W.: Weyl Transforms 


Xamb6-Descamps, S.: Block Error-Cor- 
recting Codes 


Zaanen, A.C.: Continuity, Integration and 
Fourier Theory 


Zhang, F: Matrix Theory 


Zong, C.: Sphere Packings 


Zong, C.: Strange Phenomena in Convex and 
Discrete Geometry 


Zorich, V.A.: Mathematical Analysis I 
Zorich, V.A.: Mathematical Analysis II 


